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This  manual  provides  a detailed  technical  description  of  the  SNAP  II  computer 
program.  The  manual  is  designed  such  that  a qualified  programmer  cannot  only  gain 
a full  working  knowledge  of  the  programming  logic,  but  also  alter  or  add  to  it  by 
following  specific  guidelines. 

Included  in  the  manual  arel-llf  a description  of  the  main  program,  subprograms, 
interrelations  between  programs,  and  each  of  the  variables  (code  words)  used  in  the 
programs;  cross-reference  indicating  the  program  in  which  each  variable  is  used; 
and  Sf^flow  charts  depicting  each  logical  sequence  in  the  overall  program. 
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1.  MAIN  PROGRAM  AND  INTERRELATIONSHIP 
WITH  SUBROUTINES 


SNAP  II  performs  five  main  functions  in  solving  a circuit,  as  discussed 
individually  in  Sections  1 . 1 through  1.5. 


1 . 1 READING  IN  CIRCUIT  ELEMENT  VALUES 

Circuit-element  values  are  read  into  the  main  program  after  data  statements 
have  been  entered  and  dimensioning  has  been  completed.  The  values  are  read  in 
one  at  a time,  and  each  is  checked  for  logical  and  keypunching  errors  by  subroutine 
CHECK. 

The  node  with  the  greatest  numerical  value  in  the  circuit  is  stored  as  cir- 
cuit elements  are  added.  An  incorrect  circuit  element  value  will  not  stop  the  pro- 
gram, although  errors  in  card  format  will  activate  a diagnostic  routine.  Circuit 
Element  Cards  with  errors  are  printed  in  CHECK  as  they  are  found. 

The  last  Circuit  Element  Card  contains  the  word  NOMOR,  which  causes  the 
program  to  either  1)  print  a reference  list  of  error  codes  if  Circuit  Element  Cards 
had  errors,  or  2)  rearrange  the  circuit  elements  in  numerically  increasing  order  if 
all  the  Circuit  Element  Cards  seemed  to  be  correctly  prepared.  This  rearranging 
is  done  in  subroutine  EXCH,  which  is  entered  for  each  circuit  element.  Circuit 
elements  which  were  supplied  with  numbers  are  put  in  the  position  in  the  circuit 
element  array  which  they  would  have  had  if  their  numbers  were  used  as  their  index 
in  the  array.  Circuit  elements  which  do  not  have  numbers  are  assigned  numbers 
that  are  unfilled  by  other  elements.  These  numbers  reflect  the  order  of  the  circuit 
elements  in  the  data  deck.  The  program  will  stop  if  two  circuit  elements  have  been 
assigned  the  same  number. 

The  circuit  elements  are  now  printed  in  their  new  order  by  subroutine  INOUT. 
We  do  not  use  this  portion  of  the  main  program  or  INOUT,  EXCH,  and  CHECK  again, 
unless  a NEXT  card  rather  than  a FINI  card  is  used  at  the  end  of  the  output  requests 
for  thi.s  circuit. 


1.2  FORMING  THE  EQUIVALENT  CIRCUIT 

An  output  request  card,  called  a Type  Card  in  the  User's  Manual,  is  read 
in  to  request  either  an  AC  or  DC  solution.  Subroutine  ACEQ  will  be  called  for  AC, 
subroutine  DCEQ  for  DC.  Each  of  these  subroutines  shorts  or  opens  certain  types 
of  circuit  elements;  tests  the  resulting  circuit  for  connectivity  (using  subroutine 
CONECT);  and  condenses  the  resulting  circuit  so  that  the  names  of  the  nodes  are 
sequential,  beginning  with  zero  and  progressing  through  the  maximum  node.  No  node 
numbers  are  skipped. 


1.3  ASSIGNING  VALUES  TO  CIRCUIT  ELEMENTS 

In  addition  to  denoting  the  type  of  circuit,  the  Type  Card  also  specifies  the 
type  of  solution  required.  The  choices  are:  Nominal,  Special,  Sensitivity,  Monte 
Carlo,  or  Frequency  Plot.  The  Plot  and  Sensitivity  output  have  two  choices  each, 
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and  the  Nominal  and  Special  Solutions  could  also  have  the  nonlinear  option.  Thus, 
some  programming  is  required  to  assign  proper  values  to  the  circuit  elements  for 
these  different  solutions. 

The  main  program  first  tests  to  see  if  the  solution  requires  nonlinear  circuit 
elements.  If  so,  each  unknown  node  requires  a starting  value.  These  values, 
together  with  the  upper  and  lower  limits  on  solutions  and  their  tolerances,  are  read 
in.  Subroutine  NOLIN  is  called  to  initialize  the  solutions.  Next  the  program 
branches,  depending  on  the  type  of  solution  required.  We  will  consider  each  branch 
individually. 

a.  If  the  solution  requested  is  Nominal,  the  program  reads  in  the  frequencies 
required  (in  the  AC  case)  by  calling  subroutine  READFQ.  Next,  the 
program  reads  any  nodal,  branch-current,  and  special-equation  output 
requests.  Finally,  nominal  values  are  inserted  in  the  circuit  elements 

of  the  equivalent  circuit  and  the  program  transfers  to  statement  635,  the 
beginning  of  the  solution  loop. 

b.  A Special  solution  request  causes  the  program  to  read  in  the  frequencies 
required  (in  the  AC  case)  by  calling  READFQ.  It  then  starts  the  loop  on 
the  number  of  special  solutions  this  Type  Card  specified  (columns  7-8) 
and  calls  subroutine  SPECIN,  where  node  output  requests  and  special 
values  for  the  circuit  elements  are  read  in.  The  appropriate  special  value 
is  assigned  to  each  circuit  element  and  their  values  are  printed  out.  The 
program  returns  to  MAIN,  sets  an  indicator  (NO)  to  5 to  show  a Special 
solution  request,  and  goes  to  statement  635. 

c.  A Sensitivity  solution  request  causes  the  program  to  read  in  a card  speci- 
fying the  mode  of  sensitivity  analysis  (see  below);  whether  phase, 
magnitude,  or  both  types  of  output  are  required  (in  the  AC  case);  and  the 
number  of  frequencies,  up  to  five.  Thus  for  DC  Sensitivity  requests,  this 
card  only  contains  the  mode  of  the  output.  The  main  program  then  sets 
an  indicator  (NO)  to  one  to  indicate  a Nominal  request,  and  transfers  to 
the  Nominal  part  of  the  main  program  where  nodal  and  current  output 
requests  are  read  in.  The  difference  in  the  two  modes  does  not  become 
apparent  until  output  has  been  accumulated,  as  described  in  Section  1.  6. 

d.  A Monte  Carlo  solution  request  reads  in  a card  giving  the  number  of  solu- 
tions  desired  for  each  frequency  and  the  starting  value  for  the  random 
number  generator.  If  the  number  of  solutions  requested  exceeds  999,  the 
next  solution  request  card  is  read  and  the  Monte  Carlo  solution  is  not 
done.  As  in  the  Sensitivity  solution,  the  number  of  frequencies  (up  to  five) 
and  their  values  are  read  in  on  this  card.  Up  to  five  Special  Function  out- 
put requests  are  read  in  next.  Since  only  five  output  values  are  allowed, 
the  usual  sequence  of  node  outputs,  branch  currents,  and  special  function 
output  selection  cards  does  not  apply.  Instead  the  output  requests  are 
entered  as  special  functions  and  are  requested  by  entering  equations  in 
EQUOUT. 

The  cards  specifying  the  special  functions  desired  for  output  also  specify 
upper  and  lower  bounds  on  the  Histogram  plot  routine.  Further,  an 
indicator  is  set  for  each  output  so  that  either  phase  or  magnitude  may  be 
requested.  The  program  tests  that  the  lower  limit  is  lower  than  the  upper 
limit.  If  these  cards  were  improperly  prepared,  the  program  searches 
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for  the  next  Type  Card  and  skips  this  Monte  Carlo  request.  If  no  hi^  or 
low  limits  were  set,  the  program  uses  ±20  percent  of  the  nominal  solution. 
Normally  the  nominal  solution  is  not  needed  for  a Monte  Carlo  request,  but 
in  this  case  it  is,  so  we  set  NO  = 1 and  go  to  statement  635  to  compute  the 
nominal  solution  first. 


If  the  upper  and  lower  limits  were  supplied,  then  we  begin  supplying  cir- 
cuit element  values  in  the  following  way:  The  starting  value  of  the  random 
number  generator  is  the  value  read  in  on  the  card  following  the  Type  Cards. 
The  first  sample  will  use  as  many  separately  generated  random  variables 
as  there  are  circuit  elements  in  the  equivalent  circuit.  (Admittance  and 
impedance  will  require  two  random  variables  for  their  values.)  Each  cir- 
cuit element  is  examined  to  discover  its  distribution.  If  the  distribution 
was  read  in  as  zero  on  the  Circuit  Element  Card,  the  program  provides 
the  nominal  value  to  that  element  and  goes  on  to  the  next  one.  If  the  dis- 
tribution code  was  non-zero  (i.e. , from  1 to  4),  a random  variable  is 
generated  by  subroutine  RANDU,  and  the  appropriate  subroutine  is  called 
to  transform  the  uniformly  distributed  random  variable  into  a normal 
(NORM),  lognormal  (LOGNOR),  rectangular  (RECTAN),  or  special 
(SPEC)  distribution.  The  value  produced  is  assigned  to  the  circuit  ele- 
ment. Admittance  and  impedance  elements  receive  two  random  vari- 
ables and  enter  the  appropriate  subroutine  twice.  When  all  the  circuit 
elements  have  been  assigned  values,  NO  is  set  to  6 and  the  program  trans- 
fers to  statement  635  and  begins  the  solution. 


e.  A Frequency  Plot  request  is  used  for  an  AC  equivalent  circuit  only.  It 
collects  solutions  for  various  frequencies  and  produces  a plot  of  solutions 
versus  frequency  using  a Stromberg-Carlson  CRT  display  unit. 

The  card  after  the  Type  Card  gives  the  mode  of  the  plot  desired;  whether 
a phase  or  magnitude  plot,  or  both,  is  required;  and  whether  the  frequency 
axis  is  to  be  logarithmic  or  linear.  Next  the  frequencies  are  read  in  by 
subroutine  READFQ.  The  functional  output  requests  are  read  in,  each 
containing  the  number  of  grid  lines  desired  and  the  maximum  and  mini- 
mum values  of  the  output.  If  we  are  in  mode  1,  the  PLOTF  subroutine 
computes  the  number  of  grid  lines  and  the  maximum  and  minimum  values 
of  the  output,  so  the  Functional  Output  Request  Card  will  not  contain  the 
mode  2 information.  The  PLOTF  can  only  be  used  with  nominal  values, 
so  at  this  point,  the  first  frequency  is  set  and  the  control  goes  to  the  part 
of  the  program  that  supplies  nominal  values  to  each  circuit  element. 

f.  A Nonlinear  solution  for  either  NOML  or  SPCL  output  requests  requires 
additional  cards.  Thus,  immediately  after  the  Type  Card  specifying  a 
nonlinear  circuit  is  processed,  cards  are  read  in  giving  1)  the  n^e  number 
of  a node  voltage  referenced  by  an  equation  in  EQUIN,  the  initial  value  of 
this  voltage;  and  2)  if  available,  upper  and  lower  bounds  on  the  value  of 
the  node  voltage.  The  tolerance  can  also  be  supplied.  If  upper  and  lower 
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bounds  are  not  supplied,  subroutine  NOLIN  sets  these  bounds  to  ±50  percent 
of  the  first  guess.  If  the  tolerance  is  not  supplied,  it  is  set  at  0.  2 percent 
of  the  guess,  and  if  the  tolerance  supplied  is  less  than  0.01  percent  of  the 
guess,  it  is  set  up  to  0. 2 percent.  The  program  then  proceeds  with  a 
solution  as  if  it  were  an  ordinary  Special  or  Nominal  request. 


1 . 4 OBTAINING  SOLUTIONS  TO  EQUIVALENT  CIRCUIT 

This  section  of  the  program  is  common  to  all  the  output  selections,  and  is 
broken  up  into  five  subroutines  to  reflect  the  five  steps  necessary  in  obtaining  a 
solution. 

a.  The  main  program  calls  subroutine  ASSIGN  to  complete  the  assignment  of 
values  to  and  compute  impedance  for  circuit  elements.  Here  the  value  of 
NO  directs  the  subroutine  when  it  must  find  the  imaginary  parts  of  admit- 
tance and  impedance  elements.  The  ASSIGN  subroutine  also  indicates  cur- 
rent flow  by  setting  the  sign  for  active  elements  (V,  E,  D,  I,  B,  H,  G, 
and  F types),  and  computes  the  values  for  circuit  elements  dependent  on 
other  elements  by  computing  their  value  in  an  equation  in  EQUIN. 

b.  On  return  from  ASSIGN,  the  main  program  immediately  calls  subroutine 
SOLUT.  This  subroutine  prints  the  equivalent  circuit,  if  NDEBUG  on  the 
Type  Card  of  this  output  request  was  1,  and  calls  subroutine  CURE  NT 

for  V,  E,  H and  G circuit  elements.  These  elements  are  voltage  sources, 
and  since  we  are  writing  node  equations  when  we  form  our  simultaneous 
linear  equations,  our  active  circuit  elements  must  be  current.  Subroutine 
CURENT  performs  the  conversion  from  voltage  to  current  by  putting  a 
resistor  in  parallel  with  the  voltage  in  the  circuit  and  dividing  the  voltage 
by  the  resistance.  The  value  of  the  resistance  is  1 ohm  unless  supplied 
by  the  engineer  on  the  voltage  circuit  element  card.  The  resistor  is 
entered  as  a circuit  element  in  the  equivalent  circuit,  and  this  circuit  is 
printed  out  again  if  NDEBUG  = 1. 

c.  On  return  to  the  main  program  from  SOLUT,  the  main  program  immedi- 
ately calls  subroutine  COMPUT,  where  the  node  equations  are  formed.  At 
the  end  of  COMPUT,  if  NDEBUG  = 1,  the  matrix  formed  is  printed  out. 

The  values  of  the  coefficients  are  printed  out,  not  an  algebraic  representa- 
tion of  the  equations, 

d.  Finally,  COMPUT  calls  SOLVE,  which  solves  the  simultaneous  linear 
equations  formed  by  COMPUT.  The  technique  employed  is  a modified 
Gauss  elimination  method  in  double  precision.  Real  matrices  are  solved 
in  a different  portion  of  the  subroutine  from  complex  matrices.  Because 
the  hardware  for  double  precision  varies  between  machines,  consider- 
able effort  was  made  to  avoid  unnecessary  computation  in  SOLVE.  This 
is  discussed  in  Section  3 under  subroutine  SOLVE. 

On  return  to  COMPUT  from  SOLVE,  the  solutions  are  put  in  the  appropriate 
position  in  the  OUTPUT  array  so  the  engineer's,  rather  than  the  program's, 
numbering  of  nodes  prevails.  Remember  that  on  conversion  to  the  equiva- 
lent circviit,  the  nodes  were  renumbered  so  that  no  nodes  were  skipped. 
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1 . f)  FURTHER  PROCESSING  OF  VOLTAGE  OUTPUT 
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On  return  to  the  main  program  from  COMPUT,  the  program  first  tests  to  see 
if  we  had  a singular  (or  ill-conditioned)  matrix.  This  would  indicate  that  no  solution 
was  obtainable  with  the  values  the  engineer  supplied.  (It  sometimes  indicates  that  he 
has  not  supplied  values  in  the  "high"  and  "low"  columns  — second  and  third  data  fields  — 
of  the  Circuit  Element  Card,  while  still  specifying  a sensitivity  or  special  solution 
using  those  elements. ) The  program  will  allow  two  such  "no  solution"  attempts 
before  terminating  the  run. 

a.  Nonlinear  Solution  - Next,  if  the  solution  request  was  for  a nonlinear 
circuit,  subroutine  NOLIN  is  called.  NOUN  tests  to  see  if  this  solution 
is  adequate,  or  if  new  values  for  the  unknown  outputs  are  needed.  In 
the  first  case,  INDIC  = 2 and  we  print  out  nominal  or  special  solu- 
tions. If  INDIC  = 1,  we  must  iterate  the  solution  again,  and  we 
branch  to  the  Nominal  or  Special  solution  sections  of  the  program,  which- 
ever is  appropriate.  If  the  nominal  values  are  needed,  they  are  assigned 
to  the  circuit  elements  and  the  program  returns  to  statement  635  which 
calls  subroutine  ASSIGN.  If  special  values  for  the  circuit  elements  are 
needed,  the  main  program  calls  subroutine  SPECIN  and  then  goes  to 
statement  635. 

If  INDIC  = 3,  the  iterations  have  been  used  up  without  obtaining  a solution, 
so  we  print  a diagnostic  and  search  for  the  FINI  or  NEXT  card. 

If  INDIC  = 2 or  if  we  have  a linear  circuit,  we  branch  on  the  kind  of  solu- 
tion obtained  as  in  the  linear  case. 

b.  Nominal  or  Special  Soution  — We  take  the  same  branch  and  write  the  title, 
type  of  solution,  and  solution  subtitle;  write  the  frequency  (if  this  is  an 
AC  circuit);  and  call  POLAR.  This  subroutine  calls  subroutine  OUT, 
which  computes  magnitude  and  phase.  Then,  on  return  to  POLAR,  the 
outputs  requested  by  the  engineer  are  computed  and  printed.  An  expansion 
of  the  logic  employed  in  POLAR  may  be  found  in  Section  3. 

On  return  to  the  main  program,  we  branch  if  the  solution  is  a special  one. 

For  the  Nominal  solution  of  a DC  equivalent  circuit,  we  return  to  state- 
ment 252  and  read  the  next  Type  Card.  If  this  is  an  AC  equivalent  circuit, 
the  frequency  index  test  (LF)  is  incremented  to  see  if  we  have  just  com- 
puted the  last  frequency  (if  so  we  go  to  statement  252),  and  return  to  the 
part  of  the  program  where  nominal  values  are  supplied  to  the  circuit 
elements. 

For  a Special  solution  of  a DC  equivalent  circuit,  we  1)  print  the  title  and 
call  POLAR,  2)  increment  the  counter,  ISP,  and  3)  test  to  see  if  it  exceeds 
NSP,  the  limit  supplied  by  the  engineer  on  the  number  of  Special  solutions 
he  wanted.  If  there  are  still  Special  solutions  to  be  read  in  and  performed, 
subroutine  SPECIN  is  called  to  read  the  cards  giving  the  special  values 
to  be  used.  On  return  to  the  main  program,  we  go  to  statement  635  and 
complete  the  solution.  If  we  have  read  all  the  Special  solution  cards  (when 
ISP  is  greater  than  NSP),  we  go  to  statement  252  to  read  the  next  Type 
Card. 
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For  an  AC  Special  solution  request  card,  we  1)  print  the  title  and  frequency, 

2)  call  POLAR,  3)  increment  the  frequency  index,  LF,  and  4)  test  to  see  if 
the  Special  solution  has  been  performed  for  all  the  frequencies  requested 
(is  LF  greater  than  NLF?).  If  so,  we  increment  the  Special  solution 
counter,  ISP,  and  the  program  behaves  as  it  did  in  the  DC  case,  above. 

c.  Sensitivity  — For  this  branch,  we  first  test  to  see  if  we  have  . .st  completed 
the  Nominal  solution  (needed  as  a base  for  Sensitivity  analysis)  or  if  we  are 
already  varying  the  circuit  elements.  This  is  reflected  by  the  value  of  NO. 

If  NO  = 1,  the  Nominal  solution  has  been  computed  and  is  saved  in  COMP.  t 

The  Sensitivity  variation  is  started  by  initializing  ISN,  which  controls  the 
variation  of  the  circuit  element  values.  If  the  first  element  is  not  to  be 
varied  in  a Sensitivity  analysis,  NSENSE  (1)  will  be  zero.  In  this  case  we 
increment  ISN  by  one,  test  to  see  that  we  have  more  circuit  elements  in 
the  input  circuit  (engineer  supplied,  not  the  equivalent  circuit),  and 
continue  testing  the  NSENSE  array  until  a circuit  element  is  found  that 
should  be  varied.  When  such  an  element  is  found,  we  set  all  other  values 
to  nominal  and  set  the  one  circuit  element  we  are  varying  to  its  low  value 
(columns  41-50  of  the  element's  card);  set  NO  to  2;  and  go  to  statement  635 
to  continue  the  solution. 

If  the  tested  value  of  NO  is  2 (this  is  on  return  from  finding  a low  value 
solution),  we  store  the  low  solution's  node  voltages  in  the  upper  half  of  the 
COMP  array  (the  lower  half  contains  the  nominal  voltages)  and  put  the 
high  values  in  the  ISN  circuit  element  being  varied,  with  the  remiaining 
elements  being  reset  to  their  nominal  value.  NO  is  set  to  3 and  we  go  to 
statement  635  to  obtain  a solution. 

If  the  tested  value  of  NO  is  3,  we  transfer  to  a test  to  determine  the  mode 
of  sensitivity  the  engineer  requested  for  output.  If  he  requested  mode  1 
(the  mode  which  simply  prints  the  outputs  requested  and  doesn't  order 
them  or  do  a worst  case  analysis),  we  print  the  title  if  this  is  the  first 
circuit  element  varied,  print  the  frequency  if  this  is  an  AC  circuit,  and 
then  call  subroutine  SENSPR,  which  computes  and  prints  the  outputs 
requested.  On  return  from  SENSPR,  ISN  is  incremented  and  tested  to  see 
that  all  circuit  elements  have  been  varied.  If  there  are  still  circuit 
elements  to  be  varied,  we  transfer  to  the  part  of  the  program  where  NO 
is  set  to  2 and  put  the  low  values  in  the  next  element  that  the  engineer 
wished  to  be  varied.  If  all  circuit  elements  have  been  varied  (ISN  is 
greater  than  JCOMP  after  being  incremented),  NO  is  set  to  1,  and  we  test 
to  see  if  we  are  on  mode  1 or  mode  2 Sensitivity  analysis;  since  we  are 
on  mode  1,  we  test  to  see  whether  an  AC  or  a DC  circuit  is  being  solved. 

If  we  have  just  finished  varying  a DC  circuit,  we  go  to  statement  252  and 
read  the  next  Type  Card.  For  an  AC  equivalent  circuit,  the  frequency 
index  is  incremented.  If  the  circuit  has  not  been  solved  for  all  requested 
frequencies,  we  return  to  statement  611  and  begin  a Nominal  solution  for 
the  new  frequency. 

For  a mode  2 Sensitivity  solution  (after  solving  the  circuit  for  a high  value), 
subroutine  SENSPl  is  called.  This  subroutine  stores  outputs  until  all 
circuit  elements  have  been  varied.  On  return  from  SENSPl,  ISN  is 
incremented  and  tested  to  see  if  all  circuit  elements  have  been  varied.  If 
not,  we  transfer  to  the  part  of  the  program  where  NO  is  set  to  2 and  begin 
on  the  solution  of  the  circuit  with  a new  circuit  element  set  to  its  low  value 
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j If  the  circuit  has  been  solved  for  all  low  and  high  value  combinations 

requested,  we  then  test  to  see  if  we  are  on  mode  1 or  mode  2.  Since  a 
I mode  2 branch  is  now  being  examined,  we  test  to  see  if  a worst  case 

analysis  (low  or  high)  has  been  solved.  This  is  indicated  when 

* NEXIT  = 1 or  2.  If  NEXIT  = 0 we  have  not  stored  the  worst-case  low 
values  of  the  circuit  elements  so  we  must  obtain  a worst  case  solution. 

We  also  must  order  and  print  the  preceding  variations.  For  this  purpose, 

» SENSPl  is  called  again.  On  return  from  SENSPl,  if  NEXIT  = 1,  ISN  is 

still  greater  than  JCOMP,  the  mode  is  still  2,  and  NEXIT  is  neither  0 
■ nor  5,  then  we  set  NO  = 4 (to  help  ASSIGN  find  the  correct  admittance 

j and  impedance  values)  and  transfer  to  statement  635,  where  we  begin  a 

solution  of  the  worst-case  low  variation. 

I On  return  from  COMPUT,  we  take  the  Sensitivity  solution  branch,  test  that 

* NO  is  greater  than  2 and  that  the  mode  is  2,  and  call  SENSPl  again. 

SENSPl  stores  the  low  worst-case  solution,  computes  the  high  worst-case 
values  for  the  circuit  elements,  sets  NEXIT  to  2,  and  returns  to  the  main 
program.  ISN  is  still  greater  than  JCOMP,  and  NEXIT  is  now  2,  so  we 
eliminate  the  previous  tests  and  go  to  statement  635  where  the  solution 

1 begins.  On  return  to  SENSPl,  the  worst  case  high  and  low  solutions  are 

printed  out  and  SENSPl  tests  to  see  if  there  are  more  output  requests  to 
vary  in  a worst  case  analysis.  (Remember  that  worst  case  combinations 
of  circuit  elements  are  dependent  on  the  sign  of  the  output  when  we  sub- 
j tract  high  solutions  from  low  solutions.  Thus  we  must  compute  separate 

' worst  cases  for  each  output  requested. ) If  all  outputs  have  not  been  con- 

sidered, the  next  output  values  are  ordered  and  printed  out,  worst-case 

!Iow  values  are  computed  for  the  components,  NEXIT  is  set  to  1,  and  we 

return  to  the  main  program.  If  all  the  outputs  have  been  computed, 

NEXIT  is  set  to  5 and  on  return  to  main  we  either  read  the  next  Type  Card 
j (DC  case)  or  vary  the  frequency  (AC  case). 

‘ d.  Frequency  Plot  Solution  — On  return  to  the  main  program  from  COMPUT, 

we  branch  to  a part  of  the  program  which  calls  the  PLOTF  subroutine. 

I PLOTF  stores  the  outputs  obtained,  but  doesn't  do  any  plotting  until 

LF  = NLF,  where  LF  is  the  index  of  the  present  frequency  and  NLF  is  the 
index  of  the  last  frequency  to  be  computed.  On  return  to  the  main  program 
j from  PLOTF,  the  program  increments  LF  by  one,  and  tests  to  see  if  all 

I NLF  frequencies  have  been  computed.  If  not,  we  go  to  statement  611 

where  the  nominal  values  of  the  circuit  elements  in  the  equivalent  circuit 
are  supplied.  If  all  NLF  frequencies  are  computed,  then  the  last  time 
PLOTF  was  called  a plot  of  the  outputs  was  produced,  so  we  are  through 
' with  this  output  selection.  The  program  transfers  to  statement  252  and 

reads  a new  Type  Card. 

I e.  Monte  Carlo  Solution  — On  return  to  the  main  program  from  COMPUT,  we 

branch  to  a part  of  the  program  which  tests  if  NO  = 1.  This  indicator  is 

I set  to  1 when  we  have  a Nominal  solution,  and  the  only  time  a Nominal 

solution  would  be  needed  in  a Monte  Carlo  analysis  would  be  when  we  didn't 
supply  high  and  low  limits  on  the  histogram  requested  as  the  output  for  the 
Monte  Carlo  analysis.  So  when  NO  = 1 we  must  transfer  back  to  the  part 
I of  the  program  which  initializes  a Monte  Carlo  solution.  Here  (state- 

' ment  483),  we  compute  the  nominal  values  of  those  outputs  which  do  not 

have  high  and  low  limits  and  then  set  the  high  and  low  limits  to  ±20  percent 

I 
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r 

} 

of  the  nominal  value.  Then  we  begin  the  Monte  Carlo  initialization  as 
described  in  the  beginning  of  this  section. 

If  NO  is  not  1,  we  call  subroutine  STAT,  which  stores  the  solutions  for  1 

later  printing  in  HIST,  the  histogram  plot  routine.  On  return  from  STAT, 
we  transfer  to  an  earlier  portion  of  the  program  where  IR  is  incremented. 

IR  is  the  coimter  that  keeps  track  of  the  number  of  samples  taken  so  far. 

Each  IR  generates  a new  solution  of  the  circuit  with  new  values  supplied  to 

the  circuit  elements  which  the  engineer  wanted  to  vary  randomly.  After 

incrementing  IR,  the  program  tests  to  see  if  we  have  as  many  samples  as 

the  engineer  requested  (NORV).  If  IR  is  not  greater  than  NORV,  we  1 

transfer  to  the  part  of  the  program  which  supplies  new  values  to  the  circuit 

elements.  j 

If  we  have  completed  NORV  samples,  we  go  to  statement  252  and  read  a 
new  T3rpe  Card. 

I 
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2.  CROSS-REFERENCE  OF  ARRAYS 


. Variable 

A 

AC 

ADMIT 

ALL 

ASTI 

B 

BEG 

BEND 

BIN 

BLANK 

BLANKl 

BLANK4 

BVAL 

C 

CARD 

CARLO 

CHECK 

COL8L 

COL8R 

COLMID 

COMP 

CUR 

CURR 

D 

i 

DC 

I 


Found  in 

CHECK,  ASSIGN,  SPCE,  COMPUT, 
SOLVE,  EXCH,  RECTAN,  INOUT, 
MAIN 

SPECIN,  POLAR,  MAIN 
CHECK,  ASSIGN,  MAIN 
SPECIN,  MAIN 
HIST 

CHECK,  ASSIGN,  EXCH,  RECTAN, 
PLOTF,  SOLUT,  CURENT, 
COMPUT,  MAIN 

STAT  ■’ 

PLOTF 

STAT 

READFQ 

CHECK,  MAIN 

CHECK,  MAIN 

COMPUT 

CHECK,  ASSIGN,  DCEQ,  MAIN 

EQUOUT,  CARD,  EQUIN 

CHECK,  INOUT,  MAIN 

READFQ 

HIST 

HIST 

HIST 

SPECIN,  SENSPl,  BRANCH, 
POLAR,  MAIN 

SENSPl 

CURR,  EQUOUT,  EQUIN 

CHECK,  ASSIGN,  DCEQ, 

COMPUT,  ACEQ,  MAIN 

SPECIN,  SENSPl,  MAIN 
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Variable 

Found  in 

1 

,l 

DEL 

PLOTF,  NOLIN 

- 

BELT 

NOLIN 

.1 

DEN 

ASSIGN 

DIAG 

DCEQ,  CURRENT,  ACEQ,  MAIN 

1 

DIF 

SPCE,  STAT,  NOLIN,  NORM 

DIST 

SPCE 

r 

DIV 

COMPUT,  BRANCH,  SOLVE 

DIVI 

SOLVE 

* 

1 

DIVR 

SOLVE 

1 

DSAVE 

NOLIN 

1 

DUM 

ASSIGN,  SOLUT,  CONECT,  EQUIN 

f 

DUMB 

READFQ 

DX 

PLOTF 

1 

E 

CHECK,  MAIN,  ASSIGN,  EXCH, 

SOLUT,  DCEQ,  CURENT,  COMPUT 

1 

ER 

PLOTF,  EQUOUT 

ERR 

EQUIN 

1 

ERROR 

ASSIGN 

1 

F 

CHECK,  MAIN,  ASSIGN,  EXCH, 

PLOTF 

1 . 

1 

F13 

STAT 

1 

F25 

HIST 

FINI 

MAIN 

! 

FINISH 

STAT 

FREQ 

i 

1 

MAIN,  ASSIGN,  PLOTF,  STAT, 

SENSPR,  POLAR,  EQUIN, 

READFQ,  SENSPl,  EQUOUT 

1 

! FS 

! 

STAT 

.1 

G 

CHECK,  MAIN,  ASSIGN,  EXCH, 

CURENT 

1 

GI 

MAIN 

GSAVE 

NOLIN 

1 

GUESHI 

MAIN,  NOLIN 

* 

! 

i. 
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ill 

I 

Found  in 

MAIN,  NOLIN 
NOLIN 


Variable 

GUESLO 

GUESS 


L 

f * 

i. 

r 

i: 

i; 


r 

i: 

r 

i 


I’ 

I 

I 


H CHECK,  MAIN,  ASSIGN,  EXCH, 

SOLUT,  CURENT,  COMPUT 

HI  MAIN 

HIHIST  MAIN,  STAT 


I 

II 
12 

lAD 

IC 

ICI 

ICOM 

ICOMP 


ICR 

IE 

II 

IK 

IL 

IM 

IN 

INDIC 

INL 

INODE 

INPUT 

lOUT 

lOVER 


OUT,  EXCH,  PLOTF,  CURENT, 
READFQ 

SOLUT,  DCEQ,  SENSPR,  SOLVE, 
ACEQ 

SOLVE 

CHECK,  MAIN,  ASSIGN,  EXCH, 
INOUT 

CHECK,  MAIN,  EXCH,  INOUT 

CURENT 

SENSPl 

CARD,  MAIN,  ASSIGN,  SPECIN, 
SOLUT,  DCEQ,  CURENT, 
COMPUT,  BRANCH,  POLAR, 
SENSPl,  CURR,  ACEQ 

POLAR 

CHECK,  MAIN,  ASSIGN,  EXCH, 
INOUT 

COMPUT,  SOLVE 

STAT 

MAIN 

STAT,  COMPUT 
MAIN 

MAIN,  NOLIN 
MAIN 

MAIN,  ASSIGN,  SOLUT,  CONECT, 
DCEQ,  CURENT,  COMPUT, 
POLAR,  SENSPl,  CURR,  ACEQ 

EQUIN 

DCEQ 

SPCE 
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Variable 

IP 

IPOW 

IPR 

IPRl 

IPR2 

IR 

IS 

ISN 

ISP 

1ST 

ISTAT 

IT 

IT123 

ITI 

ITYPE 

IV 

I VALUE 

IW 

IWHICH 

J 

J2 

J3 

JB 

JBOX 

JC 

JCOMP 

JI 

JNODE 


-• 

t 

I 

.t 

Found  in  | 

CHECK,  MAIN,  EXCH,  SOLUT, 

CURE NT  j 

DCEQ 

SPECIN,  SENSPR,  POLAR  ' 

SENSPR 

SENSPR 

MAIN,  STAT,  COMPUT,  SOLVE 

CHECK,  MAIN,  ASSIGN,  EXCH, 

INOUT,  SPECIN,  SENSPl 

MAIN,  SENSPR,  SENSPl 

MAIN,  SPECIN 

NORM 

MAIN,  STAT 

CHECK,  MAIN,  EXCH,  INOUT, 

STAT,  COMPUT 

READFQ  j 

CURENT 

MAIN,  ASSIGN,  SOLUT,  DCEQ,  ) 

CURENT,  COMPUT,  ACEQ  \ 

CURENT 

MAIN,  SPECIN,  SOLUT,  DCEQ,  j 

CURENT,  BRANCH,  SENSPl,  ^ 

CURR,  ACEQ 

MAIN,  STAT  I 

MAIN,  PLOTF,  SENSPR,  SENSPl 

CHECK,  EXCH,  BRANCH,  ^ 

READFQ,  SENSPl,  HIST 

PLOTF,  SOLVE  j 

STAT 

COMPUT  I 

STAT 

MAIN,  DCEQ,  CURENT,  SENSPl,  | 

ACEQ,  CURR  1 

MAIN,  INOUT,  SPECIN,  SOLUT, 

CURENT,  POLAR,  SENSPl  | 

SOLVE 

MAIN,  SOLUT,  DCEQ,  CURENT,  ; 

COMPUT,  ACEQ 
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i: 

1 

Variable 

Found  in 

JR 

SOLVE 

i 

JS 

SPECIN,  PLOTF,  STAT,  SENSPR, 

POLAR,  SENSPl 

( 

JSl 

MAIN,  SENSPR 

i 

JSC 

CURENT 

1 

JT 

STAT 

I 

K 

CHECK,  MAIN,  ASSIGN,  PLOTF, 

1 

STAT,  DCEQ,  NOLIN,  SENSPR, 

1 

COMPUT,  SOLVE,  POLAR, 

READFQ,  SENSPl,  HIST,  ACEQ 

1- 

K1 

SPECIN,  DCEQ,  CURENT,  SOLVE 

1 

K2 

SPECIN,  DCEQ,  CURENT,  SOLVE, 

HIST 

1. 

K22 

CURENT 

KB 

STAT,  COMPUT 

1 

KBR 

COMPUT 

KCOMP 

CURENT 

(■ 

KCT 

HIST 

1 , 

KE 

ASSIGN 

r 

KEl 

ASSIGN 

1 

KI 

COMPUT,  SOLVE 

KM 

STAT 

KR 

ASSIGN,  COMPUT,  SOLVE 

KS 

MAIN,  STAT,  SOLUT,  CURENT 

t 

KSENl 

MAIN 

KSEN2 

MAIN 

1 ■ 

j 

KV 

MAIN,  SPECIN 

KVl 

CURENT 

li 

KV2 

CURENT 

L 

ASSIGN,  CHECK 

i 

LI 

HIST 

i 

LABEL 

HIST,  STAT 

! l! 

LB4 

HIST 

^ in 

LEQ 

SENSPl,  PLOTF 

- ! '' 
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Variable 

Found  in 

.1 

LEQU 

MAIN,  HIST,  SENSPl,  POLAR, 

SENSPR,  STAT,  PLOTF, 

SPECIN 

1 

LF 

MAIN,  SENSPl,  POLAR,  SOLUT, 

SPECIN 

• 

LG 

MAIN 

LIG3 

STAT 

i 

J 

LIM 

NOLIN 

LIMIT 

NOLIN 

j 

1 

LL 

READFQ 

i 

LPRNT 

HIST 

• 

LTYPE 

READFQ 

M 

M2 

MAX 

MB 

MBR 

MC 

MD 

MER 

MIN4 

MINMID 

MODE 

MONT 

MS 

N 

N1 

N2 

N200 

NAC 

NAD 

NALL 


SPECIN,  PLOTF,  STAT,  DCEQ, 
ACEQ,  SENSPl,  ASSIGN 

STAT 

HIST 

COMPUT 

COM  PUT 

MAIN 

SPCE 

CHECK 

HIST 

HIST 

PLOTF,  MAIN 
MAIN 

SOLUT,  MAIN 

EQUIN,  EQUOUT,  CONECT.  CURR, 
ACEQ,  SENSPl,  SOLVE,  SPCE, 
PLOTF,  COMPUT,  DCEQ,  ASSIGN 

SOLVE 

SOLVE,  EXCH 
MAIN 

MAIN,  NOLIN 

SOLVE 

MAIN 


i 
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Variable 


Found  in 


NAM  BAD 

CHECK 

NAME 

SENSPR,  SPECIN,  POLAR,  MAIN, 
SOLUT,  SENSPl,  INOUT,  EXCH, 
CHECK 

NAMCP 

SENSPl 

NAMO 

SENSPl 

NARLO 

MAIN,  EXCH,  CHECK 

NB 

CURENT,  SPECIN,  CONECT,  STAT, 
SOLUT,  PLOTF,  COMPUT,  DCEQ 

NBETWH 

STAT 

NBETWL 

STAT 

NBOX 

STAT 

NBR 

SPECIN 

NBRAN 

MAIN 

NC 

SPECIN,  POLAR,  STAT,  MAIN, 
SENSPl,  SOLUT,  NOLIN,  COMPUT 

NCLOS 

NOLIN 

NCOMP 

CURENT,  CARD,  CURR,  MAIN, 
INOUT,  COMPUT,  EXCH,  DCEQ, 
CHECK 

NCON 

CURR,  MAIN 

NCR 

MAIN 

NCUR 

SPECIN,  POLAR,  SENSPl 

ND 

CURENT,  CONECT,  ACEQ,  MAIN, 
SOLUT,  NOLIN,  COMPUT,  DCEQ 

NDl 

CURENT,  CONECT,  ACEQ,  SOLUT, 
COMPUT,  DCEQ 

ND2 

COMPUT 

NDC 

MAIN,  COMPUT,  ASSIGN 

NDEBUG 

MAIN,  SOLUT,  NOLIN,  COMPUT 

NDIF 

SOLVE,  READFQ 

NDIM 

SOLVE 

NDIST 

MAIN,  INOUT,  EXCH,  CHECK 

NDMIT 

MAIN,  INOUT,  EXCH,  ASSIGN, 
CHECK 

NDR 

CONECT,  ACEQ,  DCEQ 

NEl 

BRANCH 

NE2 

BRANCH 
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Variable 

NEG 

NEQU 


NEQUAT 

NER 

NERROR 

NEW 

NEXIT 

NEXT 

NF 

NFl 

NFG 

NFGRID 

NFIRST 

NGRID 

NHOLD 

NIM 

NINT 

NL 

NLE 

NLF 

NLF2 

NLF3 

NLFSP 

NMl 

NM2 

NMC 

NMCl 

NND 

NO 

NOD 


Found  in 

NORM 

SENSPR,  CURENT,  BRANCH, 
POLAR,  CURR,  ACEQ,  MAIN, 
SENSPl,  SOLUT,  NOLIN, 

COMPUT,  DCEQ,  ASSIGN 

MAIN,  INOUT,  EXCH,  ASSIGN, 
CHECK 

MAIN,  COMPUT,  CHECK 
MAIN,  SOLUT,  CHECK 
HIST 

MAIN,  SENSPl 
MAIN 

CURENT,  CONECT,  STAT,  SOLUT, 
PLOTF,  COMPUT,  DCEQ 

STAT 

PLOTF 

MAIN,  PLOTF 

CONECT 

MAIN 

SPECIN,  MAIN 

MAIN 

PLOTF 

SPECIN,  MAIN 
PLOTF 

READFQ,  MAIN 

READFQ 

READFQ 

READFQ 

STAT,  SOLVE 

STAT 

MAIN,  INOUT,  EXCH,  CHECK 

CHECK 

ACEQ,  DCEQ 

CONECT,  ACEQ,  MAIN,  SENSPl, 
ASSIGN 

SENSPR,  MAIN,  SENSPl 


.1 

.1 


! 

I 


! 


I 


I 
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Variable 


NODE 


NODMAX 

NODMXl 

NODNAM 

NOME 

NOLIN 

NORV 

NOTCON 

NOUT 

NOUND 

NOUNDHI 

NOUNDLO 

NOW 

NPEC 

NPOW 

NPW 

NR 

NRDC 

NRDS 

NRM 

NRT 

NS 

NSl 

NSAV 

NSAVE 

NSCR 

NSENSE 

NSP 


Found  in 


SENSPR,  BRANCH,  POLAR,  CURR, 
ACEQ,  MAIN,  SENSPl,  SOLUT, 
INOUT,  COMPUT,  EXCH,  DCEQ, 
ASSIGN,  CHECK 

SENSPR,  BRANCH,  SPECIN, 
POLAR,  ACEQ,  MAIN,  SENSPl, 
SOLUT,  DCEQ 

ACEQ,  MAIN,  SOLUT,  COMPUT, 
DCEQ,  ASSIGN 

MAIN,  NOLIN 

MAIN 

MAIN,  NOLIN 

STAT,  MAIN 

CONECT,  ACEQ,  DCEQ 

SENSPR,  SPECIN,  POLAR,  MAIN, 
SENSPl 

STAT 

STAT 

STAT 

CONECT,  STAT 

SPECIN,  MAIN,  EXCH,  ASSIGN, 
CHECK 

CURENT,  MAIN,  SOLUT,  EXCH, 
CHECK 

MAIN,  CHECK 

STAT,  CHECK 

CHECK 

CHECK 

STAT 

STAT 

CONECT,  ACEQ 

CONECT 

READFQ 

CURENT,  CONECT,  ACEQ,  DCEQ 

CURENT,  CURR,  ACEQ,  MAIN, 
SOLUT,  DCEQ 

MAIN,  INOUT,  EXCH,  CHECK 
MAIN 
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Variable 


NSPEC 

NST 

NSUM 

NSUT 

NT 

NTl 

NTOP 

NTYPE 

NUM 

NV 

NW 

NWC 

NWORST 

NX 

NX2 

NXT 

NY 

NY2 

NY  GRID 
NZ 

OU 

OUT 

OUTPUT 


PI 

PIl 

PI2 

PLOT 

PLTLE 

PNOM 


Found  in 

SPECIN,  MAIN,  INOUT,  EXCH, 
CHECK 

SPECIN,  MAIN 
STAT 

SPECIN,  MAIN 
STAT,  EXCH 
STAT 
NORM 

SENSPR,  SPECIN,  POLAR,  ACEQ, 
MAIN,  SENSPl,  SOLUT,  INOUT, 
EXCH,  DCEQ,  CHECK 

SENSPR,  MAIN,  SENSPl,  CHECK 

SENSPR,  SPECIN,  MAIN 

SENSPl 

SENSPl 

SENSPl 

PLOTF 

PLOTF 

HIST 

PLOTF 

PLOTF 

PLOTF 

MAIN 

MAIN,  CURR,  BRANCH,  POLAR, 
SENSPR,  PLOTF 

COMPUT 

MAIN,  ASSIGN,  CURR,  BRANCH, 
POLAR,  SENSPl,  SENSPR,  OUT, 
EQUOUT,  PLOTF,  STAT,  NOLIN 

SOLVE 

SOLVE 

SOLVE 

MAIN 

MAIN 

SENSPl 
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, Variable 

Found  in 

1 

SOLVE 

1 PRl 

SOLVE 

PR2 

1 

SOLVE 

1 

R 

1 

SOLUT,  COMPUT,  CURENT, 
ASSIGN,  CHECK,  MAIN 

1 RANGE 

STAT,  PLOTF 

RCT 

HIST 

1 RE 

MAIN,  OUTPUT,  COMPUT 

RI 

SPCE 

[ RN 

STAT 

RNl 

STAT 

( RR 

SPCE 

‘ RV 

SPCE,  NORM,  RECTAN,  MAIN 

J RV2 

MAIN 

1 RVAL 

COMPUT 

1 SENS 

MAIN 

SIG 

HIST,  NORM,  STAT 

1 SIGNO 

STAT 

SPCL 

MAIN 

1 SPEC 

MAIN,  CHECK,  ASSIGN,  INOUT, 
SPECIN 

SQUAR 

STAT 

1 STAN 

STAT 

START 

MAIN,  READFQ 

j STEP 

READFQ 

STORE 

MAIN,  SENSPl,  PLOTF 

1 STR 

STAT 

' SUBTLE 

MAIN,  SENSPl,  SPECIN,  STAT, 
SOLUT 

[ SUM 

SOLVE 

SUMDIF 

NOLIN 

[ SUMI 

SOLVE 

SUMR 

SOLVE 
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Variable 

Found  in 

] 

( 

SUMTOL 

NOLIN 

1 

SV 

MAIN,  SENSPR,  SPECIN 

TENPC 

PLOTF 

THIS 

STAT 

TITLE 

MAIN,  STAT,  SPECIN,  SENSPl 

TOL 

MAIN,  NOLIN 

1 

TRANS 

MAIN,  EXCH,  INOUT,  COMPUT, 

CHECK 

, 1 

TRY 

NOLIN 

TTLPT 

PLOTF 

( 

TYPE 

MAIN,  SOLUT  STAT,  SPECIN, 

SENSPl 

V 

CHECK,  ASSIGN,  ACEQ,  CURENT, 

DCEQ,  COMPUT,  EXCH,  SOLUT, 

MAIN 

VAL 

RECTAN,  NORM,  CONECT,  SPCE, 

MAIN 

VAL2 

MAIN 

VALUE 

ASSIGN,  COMPUT 

VARI 

STAT 

VI 

SPCE 

VR 

SPCE 

VOLT 

EQUIN,  CURENT,  EQUOUT 

i 

VOLT2 

CURENT 

1 

; 

W 

ASSIGN 

■ 

WCMAXM 

SENSPl 

WCMAXP 

SENSPl 

WCN 

SENSPl 

WCSAVIN 

SENSPl 

j 

WCSAVX 

SENSPl 

1 

wcx 

SENSPl 

1 

2-12 

■1 

i 

1 

1 
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Variable 


Found  in 


X 

NOLIN 

XF 

PLOTF 

XHIGH 

CHECK,  ASSIGN,  SENSPl, 

SENSPR,  INOUT,  EXCH,  SPECIN 

XIM 

OUT 

XINT 

STAT 

XL  HIST 

MAIN,  STAT 

XLOW 

CHECK,  MAIN,  ASSIGN,  SENSPl, 
SENSPR,  INOUT,  EXCH,  SPECIN 

XM 

NORM,  SPCE 

XMAG 

MAIN,  POLAR,  SENSPl, 

PLOTF,  OUT 

XMAX 

SOLVE,  PLOTF,  MAIN 

XMCHI 

CHECK,  MAIN,  ASSIGN,  INOUT, 
EXCH,  SPECIN 

XMCLO 

CHECK,  MAIN,  ASSIGN,  INOUT, 
EXCH,  SPECIN 

XMEAN 

NORM,  STAT 

XMED 

STAT 

XMIN 

SOLVE,  PLOTF 

XN 

EXCH 

XNINT 

PLOTF 

XNOM 

CHECK,  MAIN,  ASSIGN,  SENSPl, 
CURENT,  INOUT,  EXCH,  SPECIN, 
SOLUT 

XNORM 

NORM 

XP 

STAT 

XPHS 

MAIN,  POLAR,  SENSPl,  PLOTF, 
OUT 

XR 

CURENT 

Y 

NOLIN,  COMPUT,  DCEQ,  ASSIGN, 
CHECK 

YD 

PLOTF 

YL 

HIST 

Z 

INOUT,  SPCE,  EXCH,  ASSIGN, 
CHECK 

I 
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ACKQ  (NODK,  NTYPE,  ITYPE,  INODE,  ND,  V,  D,  NODMAX,  NEQU,  DIAG,  JNODE 


ACEQ  is  called  by  the  main  program  whenever  an  AC  circuit  solution  is 
requested  and  either  1)  this  is  the  first  solution  request  after  the  Circuit  Element 
Cards  are  read  in,  or  2)  the  previous  solution  request  was  for  a DC  circuit.  ACEQ 
converts  the  read-in  circuit  to  an  AC  equivalent  circuit  by  opening  DC  current 
sources  and  shorting  DC  voltage  sources.  In  the  case  of  voltage,  ACEQ  accomplishes 
this  by  bringing  the  two  nodes  on  each  side  of  the  voltage  source  "together"  so  ftat 
a node  is  lost.  The  dropped  node  is  renumbered  to  correspond  with  the  retained 
node,  and  all  references  to  it  in  the  read-in  circuit  are  also  renumbered.  Any 
circuit  elements  in  parallel  with  the  shorted  voltage  source  necessarily  become  shorted. 

Current  sources  are  opened  by  dropping  any  circuit  elements  of  type  D from 
the  equivalent  circuit  being  formed.  Because  the  network  has  been  opened,  it  may 
now  be  unconnected,  so  subroutine  CONECT  is  called  to  test  this.  If  the  circuit  is 
now  without  power  in  all  its  subcircuits,  it  is  an  incorrect  equivalent  circuit.  A flag 
is  set  to  indicate  this  and  an  error  return  is  taken  to  the  main  program. 

Otherwise  (if  the  network  is  powered),  the  equivalent  circuit  is  constructed 
from  the  scratch  array,  NSCR,  in  which  the  nodes  in  the  equivalent  circuit  are  stored. 
Forming  the  equivalent  circuit  consists  of  going  through  the  scratch  array  and  finding 
which  nodes  are  connected  to  each  other.  ACEQ  also,  at  this  time  in  the  program, 
eliminates  any  node  numbers  that  do  not  have  links  coming  from  them.  This 
compresses  the  network  and  makes  the  "table  of  contents"  array,  NEQU,  necessary. 

Finally,  the  program  tests  that  at  least  one  node  is  named  "0"  for  ground,  i.  e.  , 
that  "0"  was  not  an  eliminated  node;  and  then  returns  to  the  main  program. 


1 

DC  circuit  element. 

DIAG 

f 

0,  usually.  If  subroutine  CONECT  finds  the  circuit  separate 
when  it  shouldn't  be,  DIAG  = 1. 

1 11 

I + 1,  necessary  when  examining  arrays  concerning  nodes, 
since  ground  is  zero. 

1 ICOMP  (400) 

An  array  of  circuit  element  numbers  which  serve  as  an  index 
between  the  read-in  circuit  and  the  equivalent  circuit. 

INODE  (1) 

One  of  the  nodes  between  which  a circuit  element  in  the 
equivalent  circuit  lies. 

ITYPE  (1) 

The  circuit  element  type  in  the  equivalent  circuit,  either  R,  L, 
C,  V,  E,  D,  I,  A,  Z,  B,  H,  F,  or  G. 

1 IVALUE  (2,400) 

A place-saver  in  the  COMMON  listing. 

JC 

The  number  of  circuit  elements  times  two  in  the  equivalent 
circuit. 

JNODE  (2,  1) 

An  array  that  serves  as  a table  of  contents  for  the  array 

describing  the  equivalent  circuit.  See  appendix  on  Network 
Storage. 


ACEQ  (Continued) 
K 

M 

N 

ND 

NDl 

NDR 

NEQU  (1) 

NND 

NO 

NODE  (2,  1) 

NODMAX 
NODMXl 
NOTCON  (51) 

NS 

NSAVE  (51) 

NSCR  (2,  201) 


I 


Temporary  storage  for  a particular  value  of  ICOMP.  This  is 
necessary  since  the  entries  of  ICOMP  are  themselves  used  as 
indices. 

Temporary  storage  of  a node  in  the  read-in  circuit. 

Temporary  storage  of  a node  in  the  read-in  circuit. 

The  maximiun  node  in  the  equivalent  circuit. 

The  number  of  nodes  in  the  equivalent  circuit  plus  one  (to 
accoimt  for  ground  being  zero). 

A counter  of  the  number  of  nodes  with  links  in  the  equivalent 
circuit.  Since  some  nodes  may  be  skipped  at  this  point  in  the 
program,  NDR  may  be  less  than  ND,  the  maximum  node  in  the  j 

equivalent  circuit. 

An  array  of  node  name  equivalences  between  the  read-in 
circuit  (contained  in  the  index  1 of  NEQU)  and  the  node 
number  of  the  equivalent  circuit.  NEQU  (5)  = 3 means  that 
node  4 of  the  read-in  circuit  is  numbered  node  3 of  the 
equivalent  circuit. 

Temporary  storage  for  ND  + 1 when  the  compressed  equivalent 
circuit  is  being  constructed. 

The  number  of  nodes  in  the  NOTCON  array. 

The  node  numbers  between  which  this  component  lies . 

NODE  (1,  J)  is  the  primary  node,  NODE  (2,  J)  is  the  secondary 
node. 

The  maximum  node  in  the  read-in  circuit. 

The  maximum  value  of  the  nodes  plus  one  in  the  read-in  circuit.  : 

An  array  of  nodes  which,  after  return  from  subroutine  CONECT, 
the  program  has  found  unconnected  to  the  nodes  in  NSAVE. 

.M 

The  number  of  nodes  in  the  NSAVE  array. 

-f 

An  array  of  connected  nodes,  as  distinct  from  NOTCON,  in  | 

the  newly  generated  equivalent  circuit. 

An  array  identical  with  the  NODE  array  (of  the  read-in  circuit)  I 

at  the  beginning  of  this  subroutine.  Gradually  NCSR  is  trans- 
formed into  an  array  containing  the  node  linkages  of  the  new 
equivalent  circuit.  ' | 


i 
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ACEQ  (Continued) 
NTYPE  (1) 


The  circuit  element  type,  either  R,  C,  L,  A,  Z,  E,  V,  D,  I, 
B,  H,  F or  G. 


1 

! 


V A DC  voltage  circuit  element. 

ACEQ  calls  subroutine  CONECT. 


I 


1 


i: 

r 

i: 

r 

I' 

r 

I 
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ASSIGN  (REQU AT,  IE,  lAD,  ADMIT,  ITYPE,  A,  Z,  V,  Y,  E,  D,  INODE,  NEQU, 

XNOM,  NODE,  NDC,  R,  C,  L.  B,  FREQ,  ERROR,  NODMXl,  NO,  NDMIT, 
IS,  SPEC,  NPEC,  OUTPUT,  H,  F,  G,  XLOW,  XHIGH,  XMCHI,  NMCLO) 

ASSIGN  is  called  by  the  main  program  for  any  solution  requested.  On  entry  to 
ASSIGN,  the  values  in  the  IVALUE  array  are  the  v^ues  supplied  by  the  main  program 
from  the  Circuit  Element  Cards  (or  from  random  number  calculation  in  the  case  of  a 
Monte  Carlo  solution),  but  are  not  necessarily  in  the  form  of  impedance  as  they  will 
be  on  exit  from  ASSIGN.  After  testing  to  see  if  the  circuit  element's  value  is  supplied 
by  an  equation  (in  which  case  the  impedance  for  that  element  is  not  calculated  until 
elements  not  part  of  the  equivalent  circuit  have  been  computed),  ASSIGN  tests  to  see  if 
the  element  is  an  A or  Z type.  If  it  is,  ASSIGN  must  supply  the  imaginary  part  of  the 
value  from  the  ADMIT  array.  If  it  is  not  a nominal  solution,  ASSIGN  must  find  which 
value  in  the  ADMIT  array  to  use  for  the  imaginary  part.  It  does  this  by  a combination 
of  being  directed  by  NO  and  by  the  value  already  in  the  real  part  of  VALUE. 

Next,  ASSIGN  tests  to  see  if  we  have  a power-source  circuit  element.  If  so, 
the  sign  of  the  value  must  reflect  current  flow  through  the  element.  If  the  circuit 
is  an  AC  one,  then  C and  L circuit  elements  must  have  a frequency  multiplier  for 
computation  of  their  impedance. 

Finally,  if  the  value  of  the  circuit  element  depends  on  an  equation  in  EQUIN, 
ASSIGN  computes  its  value;  and,  if  it  is  a power  source,  sets  its  sign  correctly. 

A Admittance- type  circuit  element. 

ADMIT  (1,  1)  Array  containing  the  imaginary  parts  of  the  A and  Z circuit 


B 

C 

D 

DEN 

DUM 

E 

ERROR 

F 

FREQ 

G 

H 

lAD 


elements.  ADMIT  (1,  I)  = circuit  element  number;  ADMIT 
(2,  I)  = nominal  imaginary  value,  etc. 

Voltage-dependent  current  source. 

Capacitor 

DC  current  source. 

Reciprocal  of  the  magnitude  of  a complex  number. 

Place  holder  in  the  COMMON  array. 

AC  voltage  source. 

A dummy  not  used  by  ASSIGN. 

Current-depencent  current  source. 

Frequency. 

Current- dependent  voltage  source. 

Voltage-dependent,  voltage  source. 

Index  controlling  ADMIT. 
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ASSIGN  (Continued) 


I COMP 

IE 

INODE  (1) 


Array  of  circuit  element  numbers  for  those  elements  in  the 
equivalent  circuit. 

Index  controlling  NEQUAT. 

Node  in  the  equivalent  circuit  that  the  Ith  circuit  element  lies 
between. 


IS 

ITYPE  (1) 
K 


Index  used  by  the  NPEC  and  SPEC  array. 

T)q3e  of  circuit  element  in  the  equivalent  circuit. 

Value  of  ICOMP  for  this  circuit  element,  i.  e.  , the  circuit 
element  number. 


KE 

KEl 

KR 

L 

M 

N 


Second  node  number  plus  1 on  the  read-in  circuit. 

Actual  KE  node  number. 

Second  node  number  plus  1 in  the  read- in  circuit. 
Inductor-type  circuit  element. 

Second  index  in  the  ADMIT  array. 

Value  of  NO  used  to  find  the  right  value  in  the  ADMIT  array. 


NDC 


0 for  an  AC  circuit,  1 for  a DC  circuit. 


NDMIT  (6,  1)  Array  containing  the  imaginary  parts  of  the  A and  Z circuit 

elements.  ADMIT  (1,  I)  = circuit  element  number,  ADMIT 
(2,  I)  = nominal  imaginary  value,  etc. 

NEQU  (1)  Array  of  node  equivalences  between  nodes  in  the  read- in  circuit 

and  nodes  in  the  equivalent  circuit. 


NEQUAT  (2,  1)  Array  containing  circuit  element  numbers  in  NEQUAT  (1,  I)  and 
equation  numbers  in  EQUIN  in  NEQUAT  (2,  I). 

NO  A variable  directing  ASSIGN  as  to  the  type  of  solution  requested. 

1 means  nominal  values  are  to  be  put  in  value  of  the  circuit 
element.  2 means  low  values  for  a sensitivity  variation  on 
ISN.  3 means  high  values  for  a sensitivity  variation  on  ISN. 

4 means  worst  case  analysis  and  we  must  test  real  part  to  get 
high.  5 means  special  analysis.  Test  real  part  to  obtain 
correct  value.  6 means  Monte  Carlo  analysis.  Imaginary 
part  already  provided. 

NODE  (2,  1)  The  primary  (NODE  (1,  1))  and  secondary  (NODE  (2,  I))  between 

which  the  1th  circuit  element  in  the  read-in  circuit  lies. 


NODMXl 


Number  of  nodes  plus  1 in  the  read-in  circuit. 
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ASSIGN  (Continued) 

NPEC  (6,  1)  SPEC  array  in  fixed  format,  used  to  hold  the  circuit  element's 

number  in  NPEC  (1,  I). 


OUTPUT  (2, 

R 

SPEC  (6,  1) 

V 

VALUE  (2,  1 
VV 

XHIGH  (1) 
XLOW  (1) 
XMCHI  (1) 
XMCLO  (1) 

XNOM  (1) 

Y 
Z 

ASSIGN  calls 


1)  Node  voltage  output  required  by  subroutine  EQUIN  for  nonlinear 
solutions . 

Resistor. 

NPEC  array  holding  special  values  for  some  of  the  circuit 
elements . 

T 

A DC  voltage  circuit  element.  ^ 

I Impedance  of  the  Ith  circuit  element.  ^ 

2 7r  times  frequency. 

High  value  of  the  circuit  element  for  sensitivity  analysis  or  a 
special  value. 

Low  value  of  the  circuit  element  for  sensitivity  analysis  or  a 
special  value. 

Standard  deviation  of  the  Monte  Carlo  distribution  (if  a normal 
or  lognormal  distribution  is  required),  or  a special  value. 

Mean  of  the  Monte  Carlo  distribution  (if  a normal  or  lognormal 
distribution  is  required),  or  a special  value. 

Nominal  values  of  the  circuit  elements.  * 

AC  current  source. 

Impedance-type  circuit  element, 
subroutine  EQUIN. 

1 
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BRANCH  (J,  NEQU,  OUTPUT,  OU,  NODE) 


BRANCH  is  called  by  subprograms  POLAR,  SENPR,  SENSPl,  and  CURR.  It 
computes  the  branch  current  throi^h  circuit  element  J,  the  first  argument.  BRANCH 
first  finds  the  nodes  between  which  the  circuit  element  lies,  makes  sure  they  have 
not  been  dropped  on  conversion  to  the  equivalent  circuit,  and  then  computes  the 
voltage  between  them.  Next,  BRANCH  finds  the  impedance  of  the  circuit  element, 
divides  it  into  the  voltage  to  obtain  the  current,  and  returns  to  the  calling  program. 


COMP  (201) 
DIV 

ICOMP  (400) 


Place-saver  for  BRANCH'S  COMMON  statement. 

Magnitude  squared  of  the  impedance.  Later,  DIV  is  the 
reciprocal  of  the  impedance  squared. 

The  array  of  circuit  element  numbers  in  the  equivalent 
circuit. 


IVALUE  (2,400) 


J 

NEl 


Impedance  of  the  circuit  elements  in  the  equivalent  circuit. 
IVALUE  (1,  I)  contains  the  real  part,  IVALUE  (2,  I)  contains 
the  imaginary  part. 

Circuit  element  nvimber  of  the  elements  whose  current  has 
been  computed. 

The  circuit  element  lies  between  two  nodes,  with  NEl  being 
the  first  node  on  the  Circuit  Element  Card. 


NE2 

NEQU  (1) 

NODE  (2,  1) 

NODMAX 
OU  (1) 

OUTPUT  (2,  1) 


The  second  node  on  the  Circuit  Element  Card. 

An  array  of  node  equivalences.  If  node  numbers  were  changed 
on  conversion  to  the  equivalent  circuit,  say  node  7 was 
changed  to  node  4,  then  NEQU  (8)  = 4. 

The  primary  (NODE  (1,  I))  and  secondary  (NODE  (2,  I))  nodes 
that  a circuit  element  lies  between  in  the  read-in  circuit. 

The  maximum  number  of  nodes  in  the  read-in  circuit. 

Temporary  storage  of  the  voltage,  later  converted  to  current 
before  return  to  the  calling  program. 

Contains  the  node  voltages.  OUTPUT  (1,  3)  is  the  real  part 
and  OUTPUT  (2,  3)  the  imaginary  part  of  node  2 output  voltage. 


BRANCH  calls  no  subroutines. 
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CARD  (M) 

CARD  is  called  by  EQUOUT  or  EQUIN  when  the  impedance  of  a particular 
circuit  element  is  needed.  Subroutine  CARD  finds  the  entry  in  the  equivalent  circuit 
and  returns  to  the  calling  program 

CARD  The  value  of  circuit  element  M. 

ICOMP  (400)  An  array  of  circuit  element  numbers  in  the  equivalent  circuit. 

NCOMP  (201)  A place  holder  in  this  subroutine. 

CARD  calls  no  subroutines. 


I 


! 


I' 

i 


i 


T" 


I 


r 

j 


r* 

1 


r 

r 
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CHECK  (J,  NODE,  NTYPE,  NUM,  IE,  NEQUAT,  IT,  TRANS,  XLOW,  XHIGH, 
XMCLO,  XMCHI,  B,  K,  C,  Z,  A,  L,  V,  Y,  E,  D,  NERROR,  NAME, 
NAMBAD,  NSENSE,  NDIST,  XNOM,  NSPEC,  NMC,  BLANK4,  lAD,  ADMIT, 
BLANKl,  IC,  NARLO,  CARLO,  IS,  NPEC,  R,  NPW,  NPOW,  IP,  H,  F,  G) 


CHECK  is  called  by  the  main  program  once  for  each  component  card  read  in. 
This  subroutine  makes  eight  checks  on  the  accuracy  of  the  circuit  element's  entries 
and  reads  in  any  additional  cards  required  by  this  element.  The  eight  checks  are: 

1 . Are  the  node  numbers  greater  than  50  ? 

2.  If  the  circuit  element  is  a B,  H,  F,  or  G type,  does  it  have  a circuit 
element  reference? 

3.  Is  the  circuit  element  type  a recognized  one? 

4.  If  this  is  an  A or  Z circuit  element,  is  there  a follow-on  card  giving 
the  imaginary  parts  of  its  values? 

5.  Is  the  Monte  Carlo  distribution  specified  as  greater  than  4? 

6.  If  the  Monte  Carlo  distribution  is  1,  2,  or  3,  is  the  mean  and  standard 

deviation  (or  if  3,  the  upper  and  lower  bounds)  supplied? 

7.  If  the  distribution  is  a 4,  are  follow-on  cards  in  the  data  deck? 

8.  If  the  distribution  is  a 4,  is  the  distribution  provided  a cumulative  one 

and  is  the  ordinate  scale  normalized  between  0 and  1 ? 


A Admittance-t3q)e  circuit  element. 

ADMIT  (6,  50)  Array  used  to  hold  the  Imaginary  parts  of  A and  Z circuit 

elements.  ADMIT  (1,  I)  is  called  NDMIT  and  holds  the 
element's  number. 


B 

BLANK  1 

BLANK  4 

C 

CARLO  (41,30) 
D 


Voltage-dependent  current  source. 

A word  containing  one  blank  used  for  comparison  when  we 
are  checking  the  follow-on  card  of  an  A or  Z circuit  element. 

A word  containing  four  blanks  used  for  comparison  on  the  A 
and  Z follow-on  cards,  NAME  array. 

Capacitor. 

Same  as  NARLO,  but  with  a floating  point  name. 

DC  current  circuit  element. 


E AC  voltage  circuit  element. 

F Current- dependent  current  source  circuit  element. 
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CHECK  ( Continued) 

G 

H 

lAD 

IC 

IE 

IP 

IS 

IT 

J 

K 

L 

MER  (10) 

NAME  (201) 

NAMBAD  (100,  2) 
NARLO  (41,  30) 

NCOMP  (201) 

NDIST  (201) 

NDMIT  (6,  50) 
NEQUAT  (40,  30) 


Current- dependent  voltage- source  circuit  element. 

Voltage-dependent  voltage  source  circuit  element. 

Index  of  the  ADMIT  array,  used  for  imaginary  parts  of  A 
and  Z circuit  elements. 

Index  of  the  CARLO  (or  NARLO  array). 

Index  of  the  NEQUAT  array  which  holds  the  circuit  element 
number  and  its  equation  number. 

Index  used  by  the  NPOW  array. 

Index  used  by  the  NPEC  and  SPEC  array. 

Index  of  the  TRANS  array  which  holds  the  circuit  element 
number  of  a B,  H,  F,  or  G element  and  its  associated 
reference  number. 

The  number  of  a circuit  element  in  the  data  deck.  If  this  is 
the  third  element  read  in,  J = 3.  Comes  from  the  main 
program . 

Available  dummy  used  in  CHECK. 

Inductor-type  circuit  element. 

Array  used  to  flag  particular  types  of  circuit  element  errors. 

Up  to  4 alphanumeric  characters  used  to  describe  the  circuit 
element,  besides  its  type. 

An  error  array  no  longer  used  by  CHECK. 

Array  holding  pairs  of  points  describing  a special  distribution 
required  by  the  circuit  element  whose  number  is  in 
NARLO  (1,  I). 

Circuit  element's  unique  number  assigned  by  the  engineer, 
or,  if  blank  on  the  Circuit  Element  Card,  assigned  by  CHECK 
to  be  200  + J. 

Array  of  either  0,  1,  2,  3 or  4 giving  the  codes  for  the 
Monte  Carlo  distributions  of  the  circuit  elements. 

Fixed  point  array  of  ADMIT,  not  used  in  CHECK. 

Array  of  circuit  element  numbers  (in  NEQUAT  (1,  I))  and  their 
associated  equation  numbers  (in  NEQUAT  (2,  I))  if  their  value 
must  be  computed.  NEQUAT  (3,  I)  is  0,  1 or  2 depending  on 
the  number  of  equation-labeling  cards  following  the  Circuit 
Element  Card.  Positions  NEQUAT  (4,  I)  through  NEQUAT  (38,  I) 
hold  the  equation  image  for  printout  in  INOUT. 
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CHECK  (Continued) 
NER  (10) 

NERROR 

NMC  (201) 

NMCl 

NODE  (2,  201) 

NPEC  (6,30) 

NPOW  (2,  1) 


Array  used  for  printing  out  the  MER  array  when  errors  are 
found. 

Error  indicator,  set  to  1 by  CHECK  if  any  errors  were  found 
on  the  Circuit  Element  Cards. 

Array  which  tells  the  program,  if  a circuit  element  has 
special  distribution  for  Monte  Carlo  analysis,  how  many 
points  are  in  the  distribution. 

Index  computed  from  NMC.  The  number  of  values  read  into 
the  NARLO  array  is  twice  the  number  of  points  in  the 
special  distribution. 

Node  numbers  between  which  this  circuit  element  lies. 

NODE  (1,  J)  is  the  primary  node,  and  NODE  (2,  J)  the 
secondary  node. 

SPEC  array  in  fixed  format,  used  to  hold  the  circuit 
element's  number  in  NPEC  (1,  I). 

Array-holding  resistor  codes  for  converting  voltage  to 
current.  NPOW  (1,  I)  holds  the  circuit  element  number. 


NPW 


( 

! 

I 

i 

i 


NR 

NRDC 

NRDS 

NSENSE  (201) 
NSPEC  (201) 

NTYPE  (201) 
NUM  (2) 

R 


Read-in  code  for  a value  of  the  resistor  to  be  used  to  convert 
voltage  sources  to  current.  This  is  stored  in  the  NPOW 
array. 

Number  of  errors  found  on  a Circuit  Element  Card. 

0 or  1 , depending  on  whether  we  have  read  in  the  imaginary 
values  of  a special  distribution  for  an  A or  Z circuit  element. 

0 or  1 , depending  on  whether  we  have  read  in  the  imaginary 
values  of  special  values  for  an  A or  Z circuit  element. 

Array  containing  ones  and  zeros,  depending  on  whether  the 
J circuit  element  is  to  be  varied  in  a sensitivity  solution. 

Array  which  tells  the  program  if  this  circuit  element  has 
special  values  in  addition  to  those  in  columns  41-80  of  the 
Circuit  Element  Card. 

Circuit  element  type,  either  R,  C,  L.  A,  Z,  E,  V,  D,  I,  B, 
H,  F,  or  G. 

Circuit  element  reference  number  for  a B,  H,  F or  G 
element.  Number  NUM  (2)  is  used  for  equation  input.  If 
NUM  (2)  =0,  no  follow-on  cards  are  needed. 

Resistor  circuit  element. 


I 
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CHECK  (Continued) 


NPEC  array  holding  special  values  for  some  of  the  circuit 
elements . 


SPEC  (6,  30) 
TRANS  (2,  20) 

V 

XHIGH  (201) 
XLOW  (201) 
XMCHI  (201) 

XMCLO  (201) 

XNOM  (201) 

Y 
Z 

CHECK  calls  no 


Array  which  holds  the  circuit  element  numbers  of  B,  H,  F 
or  G elements  and  the  circuit  element  numbers  of  their 
reference  elements. 

DC  voltage  circuit  element 

High  value  of  the  circuit  element  for  sensitivity  analysis, 
or  a special  value. 

Low  value  of  the  circuit  element  for  sensitivity  analysis, 
or  a special  value. 

Standard  deviation  of  the  Monte  Carlo  distribution  (if  a 
normal  or  lognormal  distribution  is  required),  or  a special 
value. 

Mean  of  the  Monte  Carlo  distribution  (if  a normal  or  lognormal 
distribution  is  required),  or  a special  value. 

Nominal  value  of  circuit  elements. 

AC  current  circuit  element. 

Impedance-type  circuit  element. 

subroutines. 
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C'OMPUT  (ITYPE,  INODE,  JNODE,  IT,  B,  TRANS,  V,  E,  D,  Y,  NDl,  NODE, 

OUT,  NDC,  NODMXl,  NEQU,  NDEBUG,  NER,  R,  H) 

COMPUT  is  called  by  the  main  program  after  control  has  returned  from 
subroutine  SOLUT.  COMPUT,  called  once  for  each  circuit  solution  required  by 
SNAP  II,  constructs  the  admittance  matrix  from  the  equivalent  circuit.  The  simul- 
taneous equations  are  node  equations,  and  the  right-hand  sides  of  these  equations 
are  current.  The  first  equation  formed  represents  the  current  going  into  and  out  of 
the  first  node  of  the  equivalent  circuit.  The  right-hand  side  of  the  second  equation, 
for  example,  would  equal  only  active-current  sources  attached  to  node  2 of  the 
equivalent  circuit.  (Realize  that  node  2 of  the  equivalent  circuit  might  equal  node  6, 
say,  of  the  read-in  circuit,  in  which  case  NEQU  (7)  = 2. ) 

The  coefficients  of  the  second  equation  are  formed  as  follows;  The  coefficient 
for  the  first  position  (a2i)  is  the  negative  sum  of  the  admittances  of  all  circuit 
elements  which  lie  between  nodes  one  and  two.  The  coefficient  for  the  second 
position  (a22)  is  the  sum  of  all  the  admittances  of  circuit  elements  attached  to 
node  two.  The  coefficient  for  the  third  position  is  minus  the  sum  of  all  the  admit- 
tances of  circuit  elements  which  lie  between  nodes  two  and  three.  Thus  COMPUT 
must  calculate  the  admittance  from  the  impedance  and  put  the  sums  in  the  correct 
element  of  the  matrix  A.  Since  A is  a complex,  double  precision  matrix,  and  the 
complex  indexing  and  computation  must  be  done  by  hand  in  the  program  rather  than 
by  the  computer. 

The  dependent  sources  (B,  H,  G and  F types)  require  more  programming.  If 
a dependent  source  lies  between  nodes  MB  and  KB,  and  if  its  value  depends  on  the 
current  through  a circuit  element  lying  between  nodes  JB  and  IB,  then  in  equation  MB 
and  KB  the  terms  for  nodes  IB  and  JB  will  have  a sum  added  to  and  subtracted  from  them, 
respectively.  This  sum  will  depend  on  the  type  of  dependent  source  with  which  we 
are  dealing. 

COMPUT  calls  subroutine  SOLVE  to  compute  the  node  voltages.  Before 
returning  to  the  main  program,  COMPUT  puts  these  node  voltages  in  the  array 
OUTPUT  so  they  are  in  the  right  element  for  the  read  in  node  names. 

A Admittance-type  circuit  element. 

B Voltage-dependent  current  source. 

BVAL  Value  of  a dependent  power  source. 

D DC  current  circuit  element. 

DIV  Double-precision  reciprocal  of  the  magnitude  of  a complex 

impedance. 

E AC  voltage  circuit  element. 

H Voltage-dependent  voltage  source  circuit  element. 

IB  One  of  the  nodes  between  which  the  reference  circuit  element 

of  a dependent  power  source  lies.  JB  is  the  other  node. 

ICOMP  (400)  Circuit  element  numbers  in  the  equivalent  circuit. 
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COM  PUT  (Continued) 


II 

IM 

INODE  (400) 

IR 

IT 


ITYPE  (400) 

JB 

JNODE  (2,  51) 


K 

KB 

KBR 

KI 

KR 

MB 


2*1,  indexes  the  imaginary  part  of  an  array. 

Temporary  storage  for  the  imaginary  part  of  a sum  in  double 
precision. 

One  of  the  nodes  between  which  a circuit  element  in  the 
equivalent  circuit  lies. 

II-l,  indexes  the  real  part  of  an  array. 

Index  of  the  TRANS  array  which  holds  the  circuit  element 
number  of  a B,  H,  F or  G element  and  its  associated  reference 
element  number. 

The  circuit  element  types  in  the  equivalent  circuit.  Either 
R,  L,  C,  V,  E,  D,  I,  A,  Z,  B,  H,  F,  or  G. 

One  of  the  nodes  between  which  the  reference  circuit  element 
of  a dependent  power  source  lies.  IB  is  the  other  node. 

An  array  which  serves  as  a table  of  contents  for  the  array 
describing  the  equivalent  circuit.  See  appendix  on  Network 
Storage. 

Temporary  value  of  a node  number. 

One  of  the  nodes  between  which  a dependent  power  source  lies. 
MB  is  the  other  node. 

Index  addressing  real  part  of  the  KB  element. 

Imaginary  index  of  the  Kth  element. 

Real  index  of  the  Kth  element. 

One  of  the  nodes  between  which  a dependent  power  source  lies. 
KB  is  the  other  node. 


MBR  Index  addressing  the  real  part  of  the  MB  element. 

N Reference  circuit  element  number  for  a dependent  power  source. 

NB  The  beginning  range  of  a DO  loop  which  searches  through  one 

node's  entries  in  the  equivalent  circuit's  arrays. 

NC  0 or  1 depending  on  whether  NDC  = 1 or  0 . 

NCOMP  (201)  A place  saver  in  COMPUT. 

ND  Number  of  nodes  in  the  equivalent  circuit. 


t 


I 
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COM  PUT  (Continued) 

NDl  Number  of  nodes  in  the  equivalent  circuit  plus  one. 

ND2  2+ND 

NDC  0 for  DC  circuit,  1 for  AC  circuit. 

NDEBUG  1 or  0 depending,  respectively,  on  whether  or  not  we  want 

intermediate  debugging  prints. 

NEQU  (1)  An  array  of  node  name  equivalences  between  the  read  in  circuit 

(contained  in  the  index  + 1 of  NEQU)  and  the  node  number  of  the 
equivalent  circuit.  Thus  NEQU  (5)  = 3 means  that  node  4 of  the 
read  in  circuit  is  numbered  node  3 of  the  equivalent  circuit. 

NER  0 or  1 on  return  from  subroutine  SOLVE,  depending  on  whether 

or  not  the  matrix  was  solvable. 

NF  The  end  of  a range  of  a DO  loop.  See  NB. 

NODE  (2,  1)  The  node  numbers  between  which  the  circuit  element  lies. 

NODE  (1,  J)  is  the  primary  node,  NODE  (2,  J)  the  secondary 
node. 

NODMXl  The  number  of  nodes  plus  one  in  the  read-in  circuit. 

OUT  (2,  1)  The  output  array  of  node  voltages.  OUT  (1,  7)  = the  node 

voltage,  real  part,  off  of  node  6.  This  is  called  OUTPUT  in 
other  subroutines. 

R Resistor-type  circuit  element. 

RE  Temporary  storage  for  the  real  part  of  a sum  in  double  precision. 

rVAL  The  value  of  the  reference  circuit  element  of  a dependent  power 

source. 

TRANS  (2,  20)  The  array  which  holds  the  circuit  element  numbers  of  B,  H,  F, 

or  G elements  and  the  circuit  element  numbers  of  their 
reference  elements. 

V DC  voltage  circuit  element. 

VALUE  (2,400)  The  impedance  of  the  circuit  element  in  the  equivalent  circuit. 

IVALUE  (1,  I)  = the  real  part,  IVALUE  (2,  I)  = the  imaginary 
part. 

Y AC  current  circuit  element. 

COMPUT  calls  subroutine  SOLVE. 
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CONECT  (INODE,  JNODE,  ND,  NDR,  NO,  NS,  N) 


CONECT  is  called  once  by  ACEQ  and  DCEQ  when  a new  equivalent  eircuit  is 
being  formed.  CONECT  tests  the  conneetivity  of  the  circuit  network  formed  in  ACEQ 
or  nCEO.  The  argument  N controls  whether  this  is  the  first  time  CONECT  has  been 
entered  for  this  equivalent  circuit,  in  which  case  we  find  the  first  node  in  the  JNODE 
array  and  break  the  circuit  into  two  parts:  the  portion  of  the  circuit  connected  to  the 
first  node  by  a series  of  links  (stored  in  NSAVE),  and  those  not  connected  to  the  first 
node  (stored  in  NOTCON). 

If  this  is  not  the  first  time  CONECT  has  been  entered  for  this  circuit,  then  N 
will  equal  the  node  we  are  to  begin  with.  CONECT  was  constructed  this  way  because 
if  the  calling  program  finds  that  the  NSAVE  array  doesn't  have  a power  source,  then 
the  part  of  the  circuit  with  a power  source  is  in  NOTCON,  which  may  consist  of 
several  unconnected  networks.  Thus  N is  a node  in  the  NOTCON  array  constructed 
in  the  previous  call  to  CONECT. 

The  method  used  in  CONECT  is  a simple  algorithm  which  traces  up  and  down 
the  tree  formed  by  the  network,  starting  with  the  first  node.  The  nodes  connected  to 
the  first  node  are  located,  the  first  node  and  its  connected  nodes  are  put  in  NSAVE, 
and  the  particular  node  being  traced  is  kept  track  of.  We  then  pick  a new  node  from 
those  in  NSAVE,  check  the  nodes  connected  to  it  to  find  if  they  are  already  in  NSAVE. 
If  not,  they  are  added.  We  continue  until  we  have  tested  all  the  nodes  in  NSAVE  and 
don't  find  any  new  ones. 

NOTCON  is  then  constructed  by  filling  it  with  the  nodes  in  the  equivalent  circuit 
which  were  not  in  NSAVE.  If  there  were  no  such  nodes,  NOTCON  is  empty  on  return. 

DUM  (601)  Place  holder  in  the  COMMON  statement. 

INODE  (1)  Node  in  the  equivalent  circuit.  See  appendix  on  Network  Storage. 

JNODE  (2,  1)  Table  of  contents  for  the  equivalent  circuit  arrays. 

N Number  of  the  node  to  be  used  as  the  first  node. 

NB  Starting  value  for  a DO  loop  searching  through  the  links  attached 

to  one  node. 

ND  Maximum  index  of  JNODE  in  the  equivalent  circuit. 

NDl  ND  + 1. 

NDR  Actual  number  of  nodes  in  the  equivalent  circuit. 

NF  Final  value  for  a DO  loop.  See  NB. 

NFIRST  First  or  source- attached  node  in  the  equivalent  circuit. 

NO  Index  of  NOTCON.  NO  = 0 return  to  the  calling  program  if  the 

circuit  is  connected. 

NOTCON  (50)  Array  of  nodes  not  connected  to  the  nodes  of  NSAVE.  They  are 
not  necessarily  connected  to  each  other. 
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CONECT  (Continued) 

NOW  Node  being  considered.  The  links  attached  to  this  node  are 

examined  to  determine  if  other  nodes  also  common  with  these 
links  should  be  added  to  NSAVE . 

NS  Index  of  NSAVE.  NS  = NDR  if  the  circuit  is  connected. 

NSl  Used  to  find  the  next  node. 

NSAVE  (50)  The  array  of  nodes  connected  to  each  other. 

VAL  (1202)  Place  holder  in  the  COMMON  statement. 

CONECT  calls  no  subroutines. 


I 


CURP:NT  (I,  ITYPE,  INODE,  JNODE,  HEQU,  ND,  R,  B,  KS,  DIAG,  IP,  NPOW, 
XNOM,  JCOMP,  V,  E,  H,  G)  

CURE  NT  is  called  by  subroutine  SOLUT  once  for  each  voltage  source  (whether 
a V,  E,  H,  or  G type  circuit  element),  for  each  solution  required.  If  this  is  a new 
equivalent  circuit,  CURENT  finds  the  voltage  source's  entry  in  the  NPOW  array  and 
uses  the  number  it  finds  there  as  a code  for  the  resistor  to  be  put  in  parallel  with  the 
voltage  source  in  the  equivalent  circuit.  If  this  is  not  the  first  time  the  equivalent 
circuit  has  been  used  in  this  run,  then  the  resistor  is  already  part  of  the  equivalent 
circuit  and  we  only  need  to  look  up  its  value  in  the  NSAVE  array,  which  was  con- 
structed during  the  first  solution  of  this  equivalent  circuit.  CURENT  also  sets  an 
entry  in  the  NSCR  array  to  flag  SOLUT  that  this  voltage  source  has  been  converted  to 
current  so  that  when  SOLUT  sees  the  reflection  of  this  circuit  element,  it  will  not 
call  CURENT  again. 


B 

DIAG 

E 

G 

H 

I 

I Cl 

ICOMP  (400) 
INODE  (1) 

IP 

ITI 

ITYPE  (1) 


IV 

IVALUL  (2,  1) 

JC 

JCOMP 


Voltage-dependent  current  source.  Not  used  in  this  version  of 
CURENT. 

A dummy  variable  previously  used  for  diagnostics. 

AC  voltage  source. 

Cur  rent- dependent  voltage  source. 

Voltage-dependent  volt^e  source. 

Index  of  the  voltage  source  in  the  equivalent  circuit;  there  are 
two  for  each  circuit  element,  one  of  which  is  I. 

The  voltage-source  circuit  element  number. 

The  circuit  element  numbers  referencing  the  read-in  circuit. 

One  of  the  nodes  between  which  a circuit  element  in  the 
equivalent  circuit  lies. 

The  number  of  entries  in  the  NPOW  array  and  equal  to  the 
number  of  voltage  sources  in  the  read-in  circuit. 

Type  of  voltage-source  circuit  element,  either  V,  E,  F,  or  G. 

Circuit  element  type  in  the  equivalent  circuit,  either  a R,  L,  C, 
H,  Z,  E,  V,  I,  D,  B,  H,  F,  or  G. 

Value  of  the  resistor  to  be  put  in  parallel  with  the  voltage  source 
when  we  are  computing  it  from  the  NPOW  array. 

Value  of  a circuit  element  for  this  solution  of  the  equivalent 
circuit. 

Two  times  the  number  of  circuit  elements  in  the  equivalent 
circuit. 

Number  of  circuit  elements  in  the  read-in  circuit. 
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CURENT  (Continued) 

JNODE  (2,  1)  An  array  which  directs  the  program  in  the  INODE,  ITYPE, 

IVALUE,  ICOMP  arrays.  See  appendix  on  Network  Storage. 


I 


I 

I 


I 

1 


JSC 


The  number  of  circuit  elements  in  the  read-in  circuit  plus  KS, 
the  number  of  resistors  added. 


K1 

K2 

K22 

KCOMP 

KS 

KVl 

KV2 

NB 

NCOMP  (201) 

ND 

NDl 

NEQU  (1) 


NF 

NPOW  (2,  1) 


One  of  the  nodes  (K2  is  the  other)  between  which  the  voltage 
source  lies.  (This  node  is  in  the  equivalent  circuit  and  may  not 
have  the  same  number  as  the  node  in  the  read-in  circuit. ) 

One  of  the  nodes  (K1  is  the  other)  between  which  the  voltage 
source  lies  in  the  equivalent  circuit. 

Aji  index  in  the  JNODE  array  for  use  in  inserting  the  resistors 
into  the  equivalent  circuit  storage. 

The  read-in  index  or  circuit  element  number  of  the  voltage 
source  we  are  converting  to  current. 

The  number  of  entries  in  the  NSAVE  array.  If  this  is  the  first 
solution  for  this  equivalent  circuit,  KS  = 0. 

Index  in  the  equivalent  circuit  arrays  of  the  first  use  of  the 
voltage-source  circuit  element  we  are  converting. 

Index  in  equivalent  circuit  arrays  of  the  second  use  of  the 
voltage-source  circuit  element  we  are  converting. 

The  beginning  range  of  a DO  loop  used  to  search  through  one 
node's  links  for  a particular  circuit  element. 

A diunmy  place  holder  for  COMMON  storage. 

Number  of  nodes  in  the  equivalent  circuit. 

ND,  the  number  of  nodes  in  the  equivalent  circuit. 

An  array  of  node  equivalences.  When  an  equivalent  circuit 
has  been  established,  we  renumber  the  nodes.  The  index  of 
NEQU  minus  one  is  the  old  read-in  circuit  node  number,  while 
the  value  of  NEQU  (I)  is  the  equivalent  circuit  renumbering  of 
that  node. 

The  end  of  a DO  loop  used  to  search  through  all  the  links  from 
a particular  node. 

An  array  of  voltage-source  circuit  element  numbers 
(NPOW  (1,  1))  and  the  negative  exponent  of  their  associated 
resistors  to  be  put  in  parallel  with  them. 
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CURENT  (Continued) 

NSAVE  (2,  25)  An  array  of  voltage-source  circuit  element  numbers  (NSAVE  (1,  I)) 
and  the  negative  exponent  power  of  their  parallel  resistors 
(NSAVE  (2,  I)).  NSAVE  is  constructed  when  we  first  add  the 
resistors  to  an  equivalent  circuit. 

An  array  used  to  flag  voltage  sources  in  the  equivalent  circuit 
which  have  already  been  converted  to  current  sources . This 
keeps  SOLUT  from  calling  CURENT  a second  time  for  the  same 
circuit  element  in  single  solution. 

Resistor. 

DC  voltage  source. 

The  real  value  of  the  voltage  source  for  this  solution.  VOLT 
will  change  sign  for  the  two  entries  in  IVALUE  to  reflect  the 
direction  of  current  flow . 

0 imless  the  value  for  the  voltage  source  was  supplied  by  EQUIN. 

Nominal  value  of  the  circuit  elements  in  the  read-in  circuit. 

The  value  of  the  resistor  to  be  put  in  parallel  with  the  voltage 
source  to  make  it  a current  source.  XR  is  computed  from  an 
entry  in  NPOW  (2,  I)  or  from  NSAVE  (2,  I). 

CURENT  calls  no  subroutines. 


i 


NSCR  (402) 

R 

V 

VOLT 

VOLT2 
XNOM  (1) 
XR 


CURR  (N) 


CURR  computes  branch  current  through  a circuit  element,  N,  whose  circuit 
element  number  is  CURR's  only  argument.  CURR  is  only  called  by  subroutine 
EQUOUT  when  current  is  a value  needed  by  an  equation. 


CURR 

ICOMP  (400) 


INODE  (400) 
IVALUE  (2,  400) 
N 

NCOMP  (201) 
NCON  (100) 
NEQU  (51) 

NODE  (2,  201) 
NSCR  (2,  201) 
OU  (2) 

OUTPUT  (2,  51) 


The  current  we  computed  and  are  returning  to  subroutine 
EQUOUT. 

An  array  of  the  circuit  element  numbers  of  the  elements  in  the 
equivalent  circuit.  The  engineer  may  have  requested  branch 
current  through  a circuit  element  dropped  on  conversion  to  the 
equivalent  circuit. 

A place  saver  in  the  COMMON  statement. 

A place  saver  in  the  COMMON  statement. 

The  circuit  element  number  of  the  element  through  which  we 
must  find  the  current. 

A place  saver  in  the  COMMON  statement. 

A place  saver  in  the  COMMON  statement. 

An  argument  needed  by  subroutine  BRANCH. 

An  argument  needed  by  subroutine  BRANCH. 

A place  saver  in  the  COMMON  statement. 

The  result  of  calling  BRANCH.  The  real  part  of  the  current 
is  in  OU  (1),  the  im^nary  part  is  in  OU  (2). 

An  argument  needed  by  subroutine  BRANCH. 


CURR  calls  subroutine  BRANCH. 


T 


DCEQ  (NODE,  NTYPE,  ITYPE,  INODE,  ND,  NODMAX,  E,  Y,  L,  DIAG,  C,  V,  D, 
NEQU,  JNODE) 


DCEQ  is  called  by  the  main  program  whenever  we  have  a DC  circuit  solution 
request  and  this  is  either  the  first  solution  request  after  reading  in  the  component 
cards,  or  the  previous  solution  request  was  for  an  AC  circuit.  DCEQ  converts  the 
read-in  circuit  to  a DC  equivalent  circuit.  This  is  accomplished  by  shorting  AC 
voltage  sources  and  inductors,  and  opening  AC  current  sources  and  capacitors.  In 
the  first  case,  we  lose  a node  and  those  components  in  parallel  with  the  shorted 
component.  In  the  case  of  capacitors,  we  just  "disconnect"  the  component. 

Since  opening  circuits  can  cause  the  network  to  be  unconnected,  we  call  sub- 
routine CONECT  to  test  the  connectivity  of  the  network.  If  the  network  is  now  in  two 
or  more  pieces,  some  without  a power  source,  the  pieces  without  power  are  dropped 
from  the  equivalent  circuit.  DCEQ  may  have  to  call  CONECT  several  times  to 
eliminate  all  such  pieces.  When  all  such  unpowered,  imconnected  sections  of  the 
network  are  eliminated,  we  construct  the  equivalent  circuit  from  the  scratch  array, 
NSCR,  in  which  we  have  been  flagging  the  nodes  we  wish  to  eliminate.  Forming  the 
equivalent  circuit  consists  of  going  through  the  scratch  array  and  finding  which  nodes 
are  connected  to  each  other.  DCEQ  also,  at  this  time  in  the  program,  eliminates 
any  node  numbers  which  do  not  have  links  coming  from  them.  This  compresses  the 
network  and  makes  the  "table  of  contents"  array,  NEQU,  necessary. 


Finally,  the  program  tests  that  at  least  one  node  is  designated  "0"  for  ground — 
that  is,  "0"  was  not  an  eliminated  node — and  then  returns  to  the  main  program. 

C Capacitor-type  circuit  element. 

D DC  current  circuit  element. 


DIAG 

E 

II 

ICOMP  (400) 
INODE  (1) 

lOUT 

IPOW 
ITYPE  (1) 


0 usually.  If  subroutine  CONECT  finds  the  circuit  without 
power,  DIAG  = 1 . 

AC  voltage  circuit  element. 

1 + 1,  necessary  when  examining  arrays  concerning  nodes, 
since  ground  is  zero. 

An  array  of  circuit  element  numbers  which  serve  as  an  index 
between  the  read-in  circuit  and  the  equivalent  circuit. 

One  of  the  nodes  between  which  lies  a circuit  element  in  the 
equivalent  circuit. 

0 if  no  danglers  (links  without  attachment  to  other  links  on 
one  end)  are  found  after  all  the  appropriate  links  were 
dropped;  1 if  danglers  are  found. 

Controls  whether  or  not  the  part  of  the  network  in  NSAVE  is 
powered. 

Type  of  circuit  elements  in  the  equivalent  circuit,  either  R, 
L,  C,  V,  E,  D,  1,  A,  Z,  B,  H,  F,  or  G. 
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DCEQ  (Continued) 


A place  saver  in  the  COMMON  array. 

Twice  the  number  of  circuit  elements  in  the  equivalent  circuit. 


IVALUE  (2,  400) 


JC 

JNODE  (2,  1) 

K 

K1 

K2 

M 

N ■ 

NB 

NCOMP  (201) 

ND 

NDl 

NDR 

NEQU  (1) 

NF 

NND 

NODE  (2,  1) 

NODMAX 


An  array  which  serves  as  a table  of  contents  for  the  array 
describing  the  equivalent  circuit.  See  appendix  on  Network 
Storage. 

Temporary  storage  for  a particular  value  of  ICOMP.  This 
is  necessary  since  the  entries  of  ICOMP  are  themselves 
used  as  indices. 

Temporary  storage  for  an  element  of  the  JNODE  array. 

Temporary  storage  for  an  element  of  the  JNODE  array. 

Temporary  storage  of  a node  in  the  read-in  circuit. 

Temporary  storage  of  a node  in  the  read-in  circuit. 

Beginning  range  of  a DO  loop  which  searches  through  one 
node's  entries  in  the  equivalent  circuit's  arrays. 

An  array  of  circuit  element  numbers  in  read-in  circuit  order. 
In  fact,  NCOMP  (N)  = N after  processing  by  EXCH. 

Maximum  node  in  the  equivalent  circuit. 

The  number  of  nodes  in  the  equivalent  circuit,  plus  one  (to 
account  for  ground  being  zero). 

A counter  of  the  number  of  nodes  with  links  in  the  equivalent 
circuit.  Since  some  nodes  may  be  skipped  at  this  point  in  the 
program,  NDR  may  be  less  than  ND,  the  maximum  node  in 
the  equivalent  circuit. 

An  array  of  node  name  equivalences  between  the  read-in  cir- 
cuit (contained  in  the  index  + 1 of  NEQU)  and  the  node 
number  of  the  equivalent  circuit.  Thus  NEQU  (5)  = 3 means 
that  node  4 of  the  read-in  circuit  is  numbered  node  3 of  the 
equivalent  circuit. 

The  end  of  a rai^e  of  a DO  loop.  See  NB. 

Temporary  storage  for  ND  + 1 when  we  are  constructing  the 
compressed  equivalent  circuit. 

The  node  numbers  this  component  lies  between.  NODE  (1,  J) 
is  the  primary  node,  NODE  (2,  J)  the  secondary  node. 

The  maximum  node  in  the  read-in  circuit. 
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DCEQ  (Continued) 
NODMXl 

NOTCON  (50) 

NSAVE  (50) 

NSCR  (2,  201) 

NTYPE  (1) 

V 

Y 

DCEQ  calls  subroutine 


The  maximum  value  of  the  nodes  plus  one  in  the  read-in 
circuit. 

An  array  of  nodes  which,  after  return  from  subroutine 
CONECT,  the  program  has  found  unconnected  to  the  nodes  in 
NSAVE . 

An  array  of  connected  nodes,  as  distinct  from  NOTCON,  in 
the  newly  generated  equivalent  circuit. 

An  array  identical  with  the  NODE  array  (of  the  read-in 
circuit)  at  the  beginning  of  this  subroutine.  Gradually  NCSR 
is  transformed  into  an  array  containing  the  node  linkages  of  > 

the  new  equivalent  circuit. 

The  circuit  element  type,  either  R,  C,  L,  A,  Z,  E,  V,  D,  I, 

B,  H,  F,  or  G. 

DC  voltage  circuit  element.  ^ 

AC  current  circuit  element. 


CONECT. 


EQUIN  (N,  INPUT,  DUM,  VOLT,  FREQ,  ERl^) 


EQUIN  is  called  by  subroutine  ASSIGN  when  a circuit  element's  value  is 
expressed  as  a functional  relationship  rather  than  read  in  as  a constant  from  the 
Circuit  Element  Card.  The  engineer  prepares  a GO  TO  card  and  the  equations  in 
FORTRAN  IV.  The  integer  N refers  to  the  equation  number  and  is  also  a statement 
number  of  the  FORTRAN  arithmetic  statement  prepared  by  the  engineer.  On  return 
to  ASSIGN,  the  impedance  of  the  circuit  element  is  in  INPUT. 


CARD  (M) 
CURR  (N) 

DUM 

ERR 

FREQ 

INPUT 

N 

VOLT 


Impedance  of  the  Mth  circuit  element. 

Branch  current  through  the  Nth  circuit  element,  used  only  for 
nonlinear  solutions. 

A space  saver  so  that  the  value  INPUT  has  real  and  imaginary  values. 

0 if  an  equation  exists  for  this  value  of  N;  1 otherwise. 

Frequency  at  which  equivalent  circuit  is  being  solved. 

Impedance  of  the  circuit  element  on  return  to  ASSIGN. 

The  equation  number  to  be  used  in  computing  this  circuit  element's 
value . 

The  voltage  off  the  Nth  node,  used  for  nonlinear  solutions. 


EQUIN  c Us  no  subroutines. 


EQUOUT  (N,  VOLT,  OUTPUT,  FREQ,  ER] 


EQUOUT  is  called  by  subroutines  POLAR,  SENSPR,  SENSPl,  PLOTF,  and 
ST  AT.  It  computes  functional  outputs  from  node  voltages,  frequency,  circuit  element 
values , and  branch  currents . The  engineer  prepares  the  equations , each  of  which 
must  have  a unique  sequential  statement  number. 


CARD  (M) 
CURR  (M) 
ER 

FREQ 

N 


Impedance  of  circuit  element  M. 

Branch  current  through  circuit  element  M. 

0 if  there  is  an  equation  whose  nvunber  is  in  EQUOUT.  If  ER  = 1,  this 
signals  subroutine  POLAR  that  there  is  no  equationwith  that  number  in  it. 

The  frequency  at  which  the  node  voltage  in  VOLT  was  obtained. 

The  unique  equation  number,  which  is  identical  with  one  of  the 
statement  numbers  in  EQUOUT. 


OUTPUT 

VOLT  (M) 


Result  of  the  calculation  performed  in  EQUOUT.  OUTPUT  (1)  is  the 
real  part,  OUTPUT  (2)  is  the  imaginary  part. 

Node  voltage  of  node  M. 


EQUOUT  calls  no  subroutines. 


7 
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EXCII  (1,  J,  NTYPE,  NAME,  NODE,  NSENSE,  NDIST,  XNOM,  XHIGH,  XLOW, 

XMCLO,  XMCHI,  NSPEC,  NMC,  N2,  IS,  IC,  lAD,  IE,  IT,  B,  A,  Z,  NPEC, 
NARLO,  TRANS,  NDMIT,  NEQUAT,  IP,  NPOW,  V,  E,  H,  G,  F) 

EXCH  is  called  by  the  main  program  once  for  each  circuit  element  in  the  data 
deck.  The  engineer  is  not  required  to  number  (in  columns  1-3)  all  circuit  elements, 
but  only  those  referenced  by  others.  However,  to  save  computer  storage  and 
simplify  programming,  it  is  necessary  to  number  all  the  circuit  elements  after  they 
have  been  read  in.  The  main  program  does  this,  assigning  numbers  to  the  circuit 
elements  not  numbered.  Then  subroutine  EXCH  puts  this  new  circuit  element 
number  in  the  NPOW,  TRANS,  NEQUAT,  ADMIT,  and  SPEC  arrays  where  appro- 
priate, and  puts  the  element's  entries  in  the  other  arrays  in  the  position  it  should 
have  with  its  new  number.  EXCH  is  not  called  again  once  the  circuit  has  been 
rearranged. 

A Admittance-type  circuit  element  which  requires  an  entry  in  the 

NDMIT  array  for  imaginary  values. 

B Voltage-dependent  current  source.  This  requires  a reference  ^ 

circuit  element  number  in  the  TRANS  array. 

E AC  voltage  source. 


F 

G 

H 

I 

lAD 

IC 

IE 

IP 

IS 

IT 

J 


Current-dependent  current  source. 

Current- dependent  voltage  source. 

Voltage-dependent  voltage  source. 

Position  in  the  circuit  element  arrays  that  this  element 
currently  occupies. 

Number  of  circuit  elements  with  entries  in  the  NDMIT  or  ADMIT 
array. 

Number  of  circuit  elements  with  entries  in  the  NARLO  or 
CARLO  array. 

Number  of  circuit  elements  whose  value  depends  on  an  equation 
in  subroutine  EQUIN.  Thus  there  are  IE  entries  in  NEQUAT. 

Number  of  voltage  circuit  elements  (i.  e. , the  number  of  entries 
in  NPOW). 

Number  of  circuit  elements  with  entries  in  the  SPEC  or  NPEC 
array. 

Number  of  circuit  elements  with  entries  in  the  TRANS  array. 

Circuit  element  number  and  position  in  the  element  arrays  that 
a circuit  element  will  occupy  when  EXCH  returns  to  the  main 
program . 
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EXCH  (Continued) 
N2 

NAME  (1) 

NARLO  (41,  1) 

NCOMP  (201) 
NDIST  (1) 

NDMIT  (6,  1) 
NEQUAT  (40,  1) 

NMC  (1) 

NODE  (2,  1) 
NPEC  (6,  1) 

NPOW  (2,  1) 

NSENSE  (1) 
NSPEC  (1) 

NT 

NTYPE  (1) 
TRANS  (2,  20) 


Number  over  200  used  as  the  circuit  element  number  until 
EXCH  was  called  for  this  circuit  element. 

Up  to  4 alphanumeric  characters  labeling  the  circuit  element. 

Array  of  special  distributions  to  be  used  when  varying  this 
circuit  element  in  a Monte  Carlo  solution.  NARLO  (1,  I)  con- 
tains the  circuit  element  number. 

Array  of  circuit  element  numbers.  When  EXCH  has  been  called 
for  all  the  circuit  elements,  NCOMP  (I)  = I. 

0 if  the  circuit  element  is  going  to  take  on  its  nominal  value  in 
a Monte  Carlo  solution.  If  NDIST  (1)  = 4,  there  are  special 
distribution  entries  in  array  NARLO. 

Array  of  imaginary  parts  for  any  A or  Z circuit  element. 

NDMIT  (1,  I)  contains  the  circuit  element  number. 

Array  of  circuit  element  numbers  (in  NEQUAT  (1,  I))  and  equa- 
tion numbers  for  subroutine  EQUIN  (in  NEQUAT  (2,  I))  for  those 
circuit  elements  whose  value  depends  on  functional  relationships. 
NEQUAT  (3,  I)  is  zero,  1,  or  2 depending  on  whether  0,  1 or  2 
follow-on  cards  followed  the  Circuit  Element  Card. 

Number  of  points  in  a circuit  element's  special  distribution 
(if  NDIST  (I)  = 4). 

Two  nodes  between  which  the  circuit  element  lies. 

Array  of  special  values  for  the  circuit  element.  NPEC 

(1,  I)  = the  circuit  element's  number,  and  NPEC  (2,  I)  through 

NPEC  (6,  I)  the  special  values. 

Array  of  resistor  negative  exponents  to  use  in  converting 
voltage  (V,  E,  G,  H circuit  element  types)  to  current. 

NPOW  (1,  I)  contains  the  circuit  element  number. 

0 if  the  circuit  element  is  not  going  to  be  varied  in  a sensitivity 
output;  otherwise,  NSENSE  = 1. 

An  array  of  special  values  which  the  circuit  element  can  use  in 
a special  solution. 

Temporary  storage  used  by  EXCH. 

The  letter  A,  Z,  R,  L,  C,  V,  E,  I,  D,  B,  H,  G,  or  F 
describing  the  type  of  circuit  element. 

Array  of  reference  circuit  element  numbers  for  use  with  B,  H, 
G,  and  F types.  TRANS  (1,  I)  contains  the  B,  H,  G,  or  F 
circuit  element  number. 
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EXCII  (Continued) 


DC  voltage  source. 

XHIGH  (1)  Sensitivity  high  value  of  the  circuit  element  (or  a special  value). 

XLOW  (1)  Sensitivity  low  value  of  the  circuit  element  (or  a special  value). 

XMCHI  (1)  Monte  Carlo  standard  deviation  or  a special  value. 

XTVICLO  (1)  Monte  Carlo  distribution  mean  value,  or  a special  value. 

XI<I  Temporary  storage  used  by  EXCH. 

XNOM  (1)  Nominal  value  of  circuit  element. 

Z Impedance-type  circuit  element  requiring  an  entry  in  the 

NDMIT  array. 

EXCH  calls  no  subroutines. 

1 

I 

l 

I 

I 

I 
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HIST  (NXT,  RCT,  LABEL,  LEQU,  SIG) 


HIST  is  called  by  subroutine  STAT  when  all  the  Monte  Carlo  variations  have 
been  performed  by  the  program  and  the  results  are  ready  to  be  plotted  in  a histogram. 
HIST  is  called  once  for  each  output. 

The  array  of  points  to  be  plotted  must  be  returned  to  STAT  without  changes, 
since  STAT  prints  out  a listing  of  their  range  after  calling  HIST.  Therefore  HIST 
stores  the  array  of  ordinates  in  a temporary  array  which  it  then  scales  so  the  histo- 
gram will  fit  on  a page.  This  scaling  is  accomplished  by  dividing  the  ordinates  by  2 
until  their  maxim.um  is  less  than  50.  The  abscissa  is  divided  up  two  different  ways, 
depending  on  the  number  of  samples  used  in  the  Monte  Carlo  analysis.  If  NOUV, 
the  number  of  samples  read  in  after  the  type  card,  is  less  than  300,  then  the  interval 
is  broken  up  into  12  intervals  with  a bar  at  l/2-sigma  intervals  about  the  mean.  If 
NORV  is  greater  than  300,  the  histogram  is  broken  up  into  24  intervals  with  a bar 
at  l/4-sigma  intervals.  The  labeling  is  at  1/2-signal  intervals,  however. 

HIST  starts  at  the  top  of  the  page  (the  title  has  already  been  printed  in  STAT); 
puts  blanks  in  the  printing  array,  LPRNT;  and  then  uses  one  of  two  parts  of  coding, 
depending  on  whether  we  are  printing  12  bars  or  24.  In  each  case,  the  ordinate 
array  is  tested  to  see  if  there  are  any  entries  that  would  be  full  enough  to  reach  the 
top  of  the  page.  If  so,  it  is  necessary  to  enter  asterisks  or  other  symbols  into  that 
bar  representing  that  range  of  ordinates.  When  we  have  tested  all  the  ordinates  for 
the  top  of  the  page  printout,  we  print  the  line  of  print,  go  to  the  end  of  the  print  loop, 
and  start  over  (first  blanking  the  line)  on  the  next  to  top  line.  This  procedure  is 
continued  until  the  histogram  has  been  printed. 

Finally,  the  bottom  line  is  printed  after  testing  to  see  if  there  are  12  bars  or  24. 
The  labeling  of  this  bottom  line  is  independent  of  the  munber  of  bars  printed  and 
consists  of  letters  ±SG1,  ±SG2,  etc.  , and  the  numeric  value  for  these  points. 

HIST  then  returns  control  to  STAT. 


ASTI  used  to  form  the  left  hand  margin. 

COLMID  Four  dots  used  for  the  inside  of  a double  bar  when  NXT  = 12. 


COL8L 

COL8R 

F13  (19) 
F25  (19) 


J 


K 


3 blanks  and  a 1,  thus  " 1"  used  to  form  the  bar  when 

NXT  = 24  and  the  preceding  KCT  was  zero. 

The  characters  ".  . . 1"  used  to  form  the  bar  when  NXT  = 24  and 
the  preceding  KCT  was  not  0. 

Bottom  line  of  the  histogram  if  NXT  = 12. 

Bottom  line  of  the  histogram  if  NXT  = 24. 

Index  which  keeps  track  of  the  LPRNT  index.  Every  word  of 
LPRNT  controls  four  characters  across  the  line  being  printed  so 
if  J = 5,  we  are  looking  at  the  17-20  characters  of  a line  of  printing. 

A reverse  index.  The  printing  DO  loop  goes  from  1 to  50,  but  we 
start  at  the  top  of  the  page,  so  we  must  test  the  KCT  array  to  see 
if  it  has  elements  equ^  to  50  (or  to  K).  As  the  DO  loop  proceeds 
and  IP  increases,  K decreases. 
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HIST  (Continued) 


1 
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i 

I 

1 

1 


K2 


KCT  (25) 

LI 

LABEL  (25) 

LB4 

LEQU  (2) 

LPRNT  (28) 


MAX 

MIN4 

MINMID 

NEW 

NXT 

RCT  (1) 

SIG  (7) 

YL 


Scaling  number  used  to  keep  track  of  how  many  times  we  had  to 
divide  KCT  elements  by  2 before  getting  the  maximum  KCT  element 
below  50. 

RCT  array  scaled  so  it  will  fit  in  the  50  lines  of  a computer  page. 

Characters  "*1"  used  for  the  left  hand  margin  if  a bar  must  be 
printed  there . 

Array  used  for  labeling  the  abscissa  axis.  Label  contains  blanks 
and  words  like  "^SGl",  "-SG1",  etc. 

Four  blank  characters,  used  to  initialize. 

The  label  for  this  output.  Subroutine  ST  AT  transmits  the  LEQU 
(2,  I)  and  LEQU  (3,  I)  entries  of  main  to  HIST  so  they  look  like 
LEQU  (1)  and  LEQU  (2). 

A print  array  containing  the  characters  to  be  printed  for  one  line  of 
the  histogram.  LPRNT  is  formed  by  testing  KCT  to  see  if  the  KCT 
entries  would  cause  any  histogram  printout  at  that  position  of  the 
histogram. 

Maximum  KCT  element. 

Characters  " — 1"  used  to  form  the  top  of  the  bar  when  NXT  = 24. 

Characters  " " used  to  form  the  first  part  of  the  top  of  the  bar 

when  NXT  = 12.  The  right  half  of  the  top  of  the  bar  is  formed 
using  MIN4. 

0 if  the  adjacent  bar  was  not  printed,  1 if  we  don't  need  to  print  a 
beginning  bar. 

Number  of  entries  in  RCT,  the  ordinate  array.  NXT  is  12  or  24, 
depending  on  whether  NORV  (in  the  main  program)  was  less  than 
or  greater  than  300,  respectively. 

Array  of  integers  used  to  obtain  the  ordinate  values  of  the  histo- 
gram. In  ST  AT,  when  an  output  is  within  a certain  range,  we  add 
one  to  that  range's  RCT  entry. 

Numeric  values  of  the  iSc  points,  SIG  (7)  and  SIG  (1);  ±20'  points, 
SIG  (6)  and  SIG  (2);  ±1<7  points,  SIG  (5)  and  SIG  (3);  and  the  mean 
(SIG  4)  used  for  labeling. 

Ordinate  axis  labels  giving  the  value  of  KCT.  YL  is  blank  for  four 
lines,  and  then  prints  the  value  of  the  fifth  entry  of  KCT  on  every 
fifth  line. 
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INOUT  (JCOMP,  IE,  NEQUAT,  NTYPE,  NAME,  NODE,  A,  Z,  NSENSE,  NDIST, 
XNOM,  XLOW,  XHIGH,  XMCLO,  XMCHI,  NSEC,  NMC,  lAD,  ADMIT,  IS, 
SPEC,  CARLO,  IT,  TRANS,  IC) 


INPUT  is  called  by  the  main  program  once  for  each  original  circuit  after  EXCH 
has  been  called  for  all  circuit  elements . INOUT  prints  each  circuit  element  and  the 
values  and  information  which  appeared  on  its  Circuit  Element  Card.  INOUT  checks 
to  see  if  the  circuit  element's  value  is  determined  by  an  equation  rather  than  numeric 
input,  and  checks  to  see  if  the  element  is  an  admittance  or  impedance  type.  (If  so, 
it  must  print  real  and  imaginaiy  parts  tor  it. ) Finally,  INOUT  prints  special  arrays 
formed  in  CHECK  for  special  Monte  Carlo  distributions. 


A Admittance-type  circuit  element. 

ADMIT  (G,  1)  Array  used  to  hold  the  imaginary  parts  of  A and  Z circuit 

elements.  ADMIT  (1,  I)  is  called  NDMIT  and  holds  the 
circuit  element's  number. 


CARLO  (21,  30) 
lAD 


IC 

IE 


IS 

IT 

JCOMP 

NAME 

NCOMP  (201) 
NDIST  (201) 


Same  as  NARLO,  but  with  a floating  point  name. 

Index  of  the  ADMIT  array,  used  for  imaginary  parts  of  A 
and  Z circuit  elements. 

Index  of  the  CARLO  (or  NARLO  array). 

Index  of  NEQUAT  array,  which  holds  the  circuit  element 
number  and  its  equation  number. 

Index  used  by  the  NPEC  and  SPEC  array. 

Index  of  the  TRANS  array  which  holds  the  circuit  element 
number  of  a B,  H,  F,  or  G element  and  its  associated 
reference  element's  number. 

The  number  of  circuit  elements  in  the  read-in  circuit. 

Up  to  4 alphanumeric  characters  used  to  describe  the  circuit 
element,  besides  its  type. 

The  vmique  circuit  element  number  assigned  by  the  engineer 
or  by  CHECK. 

Array  of  either  0,  1,  2,  3 or  4 giving  the  codes  for  the  Monte 
Carlo  distributions  of  the  circuit  elements. 


NDMIT  (6,  50)  The  fixed  point  array  of  ADMIT. 

NEQUAT  (40,  30)  An  array  of  circuit  element  numbers  (in  NEQUAT  (1,  I))  and 

their  associated  equation  numbers  (in  NEQUAT  (2,  I))  if  their 
value  must  be  computed.  NEQUAT  (3,  I)  to  NEQUAT  (39,  I)  is 
used  for  follow-on  cards. 


I 


I 

I 
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INPUT  (Continued) 
NMC 


1 


i 

1 

I 

I 

I 

1 

I 

I 

I 

I 


An  array  which  tells  the  program,  if  this  circuit  element  has 
special  distribution  for  Monte  Carlo  analysis,  how  many  points 
are  in  the  distribution. 

NODE  The  node  numbers  between  which  this  circuit  element  lies. 

NODE  (1,  J)  is  the  primary  node,  NODE  (2,  J)  is  the  secondary 
node. 

NSENSE  (201)  An  array  containing  ones  and  zeros  depending  on  whether  the 

J circuit  element  is  to  be  varied  in  a sensitivity  solution. 

NSPEC  (201)  An  array  which  tells  the  program  if  this  circuit  element  has 

special  values  in  addition  to  those  in  columns  41-80  of  the 
Circuit  Element  Card. 

NTYPE  (201)  The  circuit  element  type,  either  R,  C,  L,  A,  Z,  E,  V,  D,  I, 

B,  H,  F,  or  G. 

SPEC  (6,  30)  The  NPEC  array  holding  special  values  for  some  of  the  circuit 

elements. 

TRANS  (2,  20)  The  array  which  holds  the  circuit  element  numbers  of  B,  H, 

F,  or  G elements,  and  the  circuit  element  numbers  of  their 
reference  elements. 

XHIGH  (201)  High  value  of  the  circuit  element  for  sensitivity  analysis,  or 

a special  value. 

XLOW  (201)  Low  value  of  the  circuit  element  for  sensitivity  analysis,  or  a 

special  value. 

XMCHI  (201)  Standard  deviation  of  the  Monte  Carlo  distribution  (if  a normal 

or  lognormal  distribution  is  required),  or  a special  value. 

XMCLO  (201)  Mean  of  the  Monte  Carlo  distribution  (if  a normal  or  lognormal 

distribution  is  required),  or  a special  value. 

XNOM  (201)  Nominal  values  of  the  circuit  elements. 

Z Impedance-type  circuit  element. 

INOUT  calls  no  subroutines. 
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NOLIN  (NODNAM,  GUESS,  TOL,  NEQU,  OUTPUT,  ND,  INDIC,  NDEBUG,  NOUN, 

NAC,  GUESLO,  GUESHI) 

NOUN  is  called  by  the  main  program  once  at  the  beginning  of  a solution  to 
initialize  the  first  guesses  into  the  output  array,  and  is  subsequently  called  each 
time  the  main  program  has  a new  solution.  NOUN  is  only  called  when  a nonlinear 
set  of  equations  is  input  through  subroutine  EQUIN. 

On  initialization,  if  the  engineer  did  not  supply  upper  and  lower  bounds,  they 
are  set  to  ±50  percent  of  the  first  approximation  supplied  by  the  engineer.  If  no 
tolerance  limit  was  set,  the  tolerance  is  set  to  0.  2 percent  of  the  first  approxima- 
tion. If  the  tolerance  provided  was  smaller  than  0.  01  percent  of  the  first 
approximation,  it  is  set  at  0.  2 percent  of  the  first  approximation. 

After  the  first  solution,  using  the  first  approximation,  has  been  found,  the 
program  finds  the  "area"  that  the  solution  falls  into  and  decides  on  a delta  to  vary 
the  next  approximation.  The  areas  are  defined  as  follows: 

Area  1:  The  difference  between  the  approximation  used  and  the  output  obtained  ^ 

from  this  approximation  is  less  than  ±20  times  the  tolerance. 

I 

Area  2:  Between  area  1 and  the  high  or  low  limits.  ■ 

Area  3;  Above  or  below  the  high  or  low  limits. 

The  logic  used  to  obtain  a solution  differs  considerably  for  these  three  areas, 
and  depends  on  the  preceding  delta,  the  differences,  the  preceding  approximations, 
and  the  resulting  outputs . 

When  the  square  of  the  sum  of  the  differences  (between  the  output  obtained  and 
the  approximation  used)  of  the  nonlinear  outputs  is  less  than  the  square  of  the  sum 
of  the  tolerances,  we  have  a solution  and  return  to  the  main  program  with  INIC  set 
to  2 to  indicate  this. 

If  we  have  iterated  more  than  20  times  the  number  of  nonlinear  outputs,  we 
return  to  the  main  program  with  the  INDIC  set  to  3 to  indicate  that  no  solution  was 
obtained  and  that  we  should  stop  iterating. 

DEL  Temporary  storage  for  DELTA  (I). 

DELT  (9)  Increment  added  to  TRY  (1,  I,  2)  to  obtain  the  next  approxima- 

tion. DELT  (J)  is  the  delta  used  for  the  node  in  NODNAM  (J). 

DIF  (9,  2)  Array  containing  the  differences  between  the  output  obtained  and  j 

the  approximation  used.  DIF  (3,  1)  is  this  difference  for  the  j 

node  found  in  NODNAM  (3) . The  difference  between  the 
approximation  used  to  obtain  the  present  output  and  the  previous 
approximation  (in  TRY  (1,  3,  3))  is  found  in  DIF  (3,  2).  On  j 

return  to  the  main  program,  when  we  have  decided  on  a new 
approximation,  DIF  (3,  1)  contains  the  difference  between 
approximations,  not  the  approximations  and  outputs  themselves. 

DSAVE  (9)  Not  used  in  this  version  of  NOUN. 
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NOLIN  (Continued) 


GSAVE  (2,  9)  Previous  output  voltage  obtained  for  a particular  node. 

GSAVE  (1,  3)  is  the  previous  output  of  the  node  in  NODNAM  (3). 


GUESHl  (2,  1) 
GUESLO  (2,  1) 


GUESS  (2,  1) 

INDIC 


K 


High  limits  of  the  nodes  in  NODNAM.  See  GUESLO. 

Low  limits  of  the  nodes  in  NODNAM.  The  first  unknown  node 
is  in  GUESLO  (1,  1).  To  find  which  node  it  is  in  the  equivalent 
circuit,  we  must  look  in  NODNAM  (1). 

First  approximations  of  the  values  of  the  nodes  in  NODNAM. 
Always  supplied  by  the  engineer. 

Indicator  whose  value  flags  the  main  program  as  to  progress  in 
NOLIN,  INDIC  = 0 means  this  is  the  first  time  we've  entered 
NOLIN  and  initialization  is  necessary.  INDIC  = 1 means  we  are 
still  trying  to  achieve  a solution.  INDIC  = 2 means  we  have  a 
solution,  and  INDIC  = 3 means  we  have  not  achieved  a solution 
and  are  giving  up. 

Node  number  of  the  node  voltage  we  are  presently  approximating. 


LIM  A counter  which  counts  the  number  of  iterations  performed. 

When  LIM  is  greater  than  LIMIT,  INDIC  is  set  to  3 and  we 
return  to  main. 


LIMIT 

NAC 

NC 

NCLOS  (9) 


ND 

NDEBUG 

NEQU 

NODNAM  (1) 
NONLIN 


Number  of  iterations  possible  before  setting  INDIC  = 3 and 
giving  up  on  a solution.  At  the  present  time,  LIMIT  is  set  to 
20  times  the  number  of  nonlinearities. 

0 if  this  is  a DC  circuit.  If  NAC  = 1 , this  is  an  AC  equivalent 
circuit,  which  can't  be  in  NOLIN. 

Temporary  storage  for  NCLOS  (I). 

Tag  which  is  1,  2,  or  3 (plus  or  minus)  depending  on  what  area 
the  solution  was  from  the  previous  iteration.  Thus  if 
NCLOS  (2)  = 1,  it  means  the  output  from  the  node  in  NODNAM  (2) 
is  less  than  20  times  the  tolerance  from  the  approximation  used 
to  obtain  this  output  (i.  e. , the  voltage  area  is  area  1). 

Number  of  nodes  in  the  equivalent  circuit. 

Read- in  variable  from  the  main  program.  NDEBUG  = 0 when 
intermediate  results  are  not  wanted.  NDEBUG  = 1 when  the 
initial  values  of  NOLIN  are  not  to  be  printed  out. 

Array  not  needed  in  this  version  of  NOLIN. 

Node  numbers  of  the  nodes  whose  output  is  dependent  on  other 
nodes,  and  which  thus  must  be  solved  for  iteratively  in  NOLIN. 

Number  of  dependent  nodes  in  this  equivalent  circuit. 
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NOUN  (Continued) 

OUTPUT  (2,  1)  Node  voltages  of  the  nodes  in  the  read-in  circuit.  The  output 

voltage  of  node  6 is  found  in  OUTPUT  (1,  7).  Only  real  voltage 
is  considered  in  NOLIN,  so  OUTPUT  (2,  7)  is  not  used. 

SUMDIF  (2)  Sum  of  the  squares  of  the  differences  when  they  are  still  the 

differences  between  the  outputs  obtained  and  the  approximations 
used  to  obtain  them . 

SUMTOL  Sum  of  the  squares  of  the  tolerance,  used  to  compare  with 

sumdif  to  see  if  we  have  a solution. 

TOL  (1)  Value  for  each  unknown  nonlinear  node  which  tells  us  how  close 

to  a solution  we  must  come  before  we  can  call  it  a solution. 

TRY  (2,  9,  3)  Array  of  approximations  and  outputs  obtained.  Three  previous 

values  of  TRY  have  been  kept,  but  the  third  one  is  no  longer 
used  in  this  version  of  NONLIN.  While  room  for  imaginary 
parts  are  available  in  TRY  (2,  I,  J),  it  is  not  used.  As  an 
example,  TRY  (1,  3,  2)  contains  the  value  used  to  obtain  this 
output  voltage  for  the  node  contained  in  NODNAM  (3).  Thus  it 
contains  the  previous  approximation.  TRY  (1,  3,  1)  will  con- 
tain a new  approximation  on  return  to  main,  but  until  this  is 
set,  it  will  contain  the  voltage  obtained  from  the  preceding 
approximation. 

X +1  or  -1.  X,  like  Y,  is  used  to  find  if  a sign  has  changed. 

Y +1  or  -1.  Y,  like  X,  is  used  to  find  if  a sign  has  changed. 

NOLIN  calls  no  subroutines. 
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NORM  (RV,  SIG,  XMEAN,  VAL,  1ST) 

NORM  is  called  by  the  main  program  when  a Monte  Carlo  solution  is  requested 

and  when  at  least  one  circuit  element  in  the  equivalent  circuit  must  be  varied 

according  to  a normal  distribution.  NORM  is  called  once  for  each  such  circuit 

element,  except  in  the  case  of  A and  Z circuit  element  types,  when  it  is  called  twice. 

DIF  Difference  between  the  closest  value  of  XNORM  and  RV. 

1ST  0 the  first  time  NORM  is  entered  during  a run.  It  causes  NORM  to 

compute  the  abscissas  of  the  normal  distribution. 

NEC  0 if  RV  is  less  than  0.  5;  1 when  RV  is  greater  than  0.5.  When  the 

normalized  random  variable  is  found,  and  if  NEG  = 1,  VAL  is  rescaled. 

NTOP  A counter  which  varies  from  1 to  19  and  causes  NORM  to  vary  the 

normalized  random  variable  in  increments  of  0. 1 when  they  would  be 
greater  than  8.  5 or  less  than  1.5. 

RV  The  uniformly  distributed  random  variable  supplied  to  NORM  and 

returned  to  main. 

SIG  The  standard  deviation  of  the  normal  distribution  desired.  Supplied 

by  main. 

VAL  The  normally  distributed  random  variable  computed  by  NORM  and 

returned  to  main. 

XM  Slope  of  the  interpolated  curve. 

XMEAN  Mean  of  the  normal  distribution  desired.  Supplied  by  main. 

XNORM  (2,  57)  the  array  containing  the  normal  distribution;  the  first  entry 
is  the  abscissa,  the  second  is  the  ordinate. 

NORM  calls  no  subroutines. 


I 

I 

( 

3-39 

1 


f 

\ 


OUT  (XMAG, 

XPHS,  OUTPUT) 

OUT  is  called  by  subroutines  POLAR,  SENSPl,  PLOTF,  ST  AT,  and  SENSP. 

Its  only  function  is  to  convert  rectangular  coordinates  to  polar  coordinates,  which  is 
necessary  for  AC  output.  The  sign  of  the  magnitude  for  DC  output  is  preserved  in 

OUT. 

I 

Always  1,  and  only  necessary  because  of  a previous  version  of  OUT. 

OUTPUT 

The  rectai^ular  coordinates  of  the  voltage  output.  The  real  part  is 
in  OUTPUT  (1),  the  imaginary  part  is  in  OUTPUT  (2). 

J 

RE 

Temporary  storage  of  OUTPUT  (1,  I). 

XIM 

Temporary  storage  of  OUTPUT  (2,  I). 

XMAG 

The  magnitude  of  the  voltage  output  in  the  AC  case;  for  DC  circuits  it 
is  simply  OUTPUT  (1).  This  is  computed  in  OUT. 

- 

XPHS 

The  phase  of  the  voltage  in  the  AC  case.  XPHS  = 0 for  DC  equivalent 
circuits.  XPHS  must  be  adjusted  by  180  degrees  if  it  lies  in 
quandrant  II  or  III. 

OUT  calls  no  subroutines.  * 
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PLOTF  (LF,  NLF,  LEQU,  JS,  OUTPUT,  IWHICH,  FREQ,  MODE,  LG,  XMAX, 
XMIN,  NYGRID,  NFGRID,  TTLPT) 

PLOTF  is  called  by  the  main  program  and  used  when  a frequency  plot  solution 
is  requested.  If  LF  is  less  than  NLF,  we  are  still  computing  outputs  at  different 
frequencies  and  are  not  ready  to  plot.  Instead,  we  compute  the  output  desired  from 
EQUOUT,  convert  it  to  magnitude  and  phase  in  OUT,  and  store  it  in  the  array  STORE 
depending  on  the  value  of  IWHICH. 

When  LF  = NLF,  we  can  plot  the  outputs  saved  in  STORE.  If  MODE  = 1, 
maximum  and  minimum  values  must  be  computed  for  each  output  and  for  the  frequency 
range.  These  values  are  then  adjusted  so  that  1)  the  grid  lines  fall  on  non-fractional 
points,  and  2)  the  plot  is  centered  on  the  space  available,  and  does  not  run  up  to  the 
top  and  bottom  of  the  plot.  If  MODE  = 2,  we  use  the  limits  set  by  the  engineer,  but  we 
do  test  to  see  that  there  are  no  more  than  10  frequency  grid  lines,  since  labeling  is 
clumsy  and  cramped  if  we  do.  PLOTF  calls  subroutines  described  in  the  North 
American  Aviation,  Inc.  , Engineer's  Computing  Manual  for  the  General  Dynamics 
SC-4020  Plotter. 

A separate  section  of  PLOTF  is  used  for  logarithmic  plots. 

PLOTF  returns  to  the  main  program  when  either  a group  of  outputs  for  one 
frequency  have  been  stored  (if  LF  is  less  than  NLF),  or  when  all  the  outputs  have 
been  plotted  (when  LF  equals  NLF). 


B 


BEND 


BSl  ART 


DEL 


DX 


ER 


FREQ  (100) 


A workii^  variable  used  by  PLOTF  to  obtain  the  correct 
maximum  and  minimum  v^ues  of  the  output  when  MODE  = 1 . 

The  final  output  value  plus  a scaling  factor  to  assure  centering 
and  nonfractional  grid  lines.  BEND  = XMAX  if  MODE  = 2. 

Starting  value  of  the  output  axis.  BSTART  = XMIN  when 
MODE  = 2,  otherwise  it  is  calculated. 

The  interval  between  grid  lines  along  the  output  axis  This 
is  computed  in  both  the  MODE  1 and  MODE  2 cases. 

The  floating-point  form  of  the  increment  used  between  grid 
lines  along  the  frequency  axis. 

Error  flag  resulting  from  a bad  equation  in  EQUOUT.  Not 
tested  in  PLOTF. 

The  range  of  frequencies  used. 

The  frequencies  at  which  the  solutions  have  been  obtained. 
(Comes  from  the  main  program. ) 

A variable  which  controls  the  number  of  labeled  frequency 
grid  lines.  When  1 = 1,  every  grid  line  along  the  frequency 
scale  will  be  labeled.  If  I = 10.  every  10th  grid  line  will  be 
labeled. 
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ivvincH 


J2 

JS 


K 


LEQ 

LEQU  (3,  1) 

LF 

LG 


M 

MODE 


N 

NB 

NF 

NFG 

NFGRID 


1,  2,  or  3,  coining  from  the  main  program.  IWHICH  = 1 
indicates  that  all  outputs  are  plotted  in  magnitude;  IWHICH  = 2 
indicates  that  all  outputs  are  plotted  in  phase;  and  IWHICH  = 3 
indicates  that  both  magnitude  and  phase  are  to  be  plotted. 

A variable  which  controls  the  number  of  labeled  output  grid 
lines.  It  acts  as  I does,  but  on  the  output  axis. 

The  number  of  outputs  requested.  If  IWHICH  = 1 or  2,  there 
will  be  JS  plots  produced.  If  IWHICH  = 3,  there  will  be  2*JS 
plots.  Supplied  by  the  main  program. 

The  value  of  J from  I to  JS  (the  number  of  outputs  desired)  if 
IWHICH  = 1 or  2.  If  IWHICH  = 3,  K will  also  take  on  the 
values  J + 5 from  1 to  JS. 

Temporary  storage  for  LEQU  (1,  I),  which  contains  the 
equation  number  of  the  output  desired  (defined  in  PLOT  F) . 

Output  equation  number  (LEQU  (1,  I))  and  title  (LEQU  (2,  I) 
and  LEQU  (3,  I))  for  each  output.  Supplied  by  the  main  program. 

The  frequency  index,  supplied  by  the  main  program. 

0 or  1.  If  LG  = 0,  the  frequency  axis  will  be  laid  out 
logarithmically.  If  LG  = 1,  the  frequency  axis  is  linear. 

LG  comes  from  the  main  program. 

A variable  which  causes  output  grid  lines  to  be  retraced  for 
emphasis.  PLOTF  sets  them  to  the  output  values  of  J2. 

1 or  2,  comes  from  the  main  program.  If  MODE  = 1,  PLOTF 
computes  the  output  maximum  and  minimum  and  number  of 
grid  lines.  If  MODE  = 2,  these  are  supplied  by  the  main 
program. 

A variable  which  causes  the  frequency  grid  lines  to  be  retraced 
for  emphasis.  PLOTF  sets  them  to  the  value  of  I. 

A temporary  storage  variable  used  when  B/DEL  must  be 
truncated  to  obtain  a nonfractional  grid  interval. 

Half  the  number  of  frequeneies.  It  is  used  to  compute  the  grid 
size  of  the  frequency  axis  when  MODE  = 1. 

The  number  of  frequency  grid  lines  adjusted  to  a number 
under  11  if  too  many  grid  lines  were  specified. 

The  number  of  grid  lines  along  the  frequency  axis.  This  is 
computed  in  PLOTF  if  MODE  = 1;  it  is  supplied  by  the  main 
program  if  MODE  = 2,  but  may  be  corrected  by  PLOTF  if  it  is 
too  large.  It  is  independent  of  output. 


I 
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PLOTF  (Continued) 

NINT 

The  number  of  intervals  on  the  output  axis.  NINT  is  used  when 

MODE  = 1 . 

NLE 

The  last  frequency  index,  supplied  by  the  main  program . If 

LF  = NLF,  we  store  the  last  output  and  begin  plotting  STORE. 

NX 

A variable  used  by  NXV  fimction  to  connect  the  points  of  the 
plot. 

NX2 

A variable  similar  to  NX. 

NY 

A variable  used  by  NYV  function  to  connect  the  points  of  the 
plot. 

NY2 

A variable  similar  to  NY. 

NYGRID  (1) 

The  number  of  grid  lines  along  the  output  axis.  This  is 
computed  in  PLOTF  if  MODE  = 1;  otherwise  it  is  supplied  by 
the  main  program. 

OU 

Computed  output  resulting  from  calling  EQUOUT.  OU  (1) 
contains  the  real  part,  OU  (2)  = the  imaginary  part  of  the 
output. 

1 

1 

i 

{ 

OUTPUT  (2,  50) 

The  node  voltage  outputs  (Re  = OUTPUT  (1,  I), 

Im  = OUTPUT  (2,  I))  resulting  from  tl.>s  frequency's  solution 
of  the  circuit.  The  real  part  of  the  voltage  at  node  4 is  found 
in  OUTPUT  (1,  5).  OUTPUT  comes  from  the  main  program. 

RANGE 

The  difference  between  the  maximum  and  minimum  values  for 
a particular  output.  When  MODE  = 1 , RANGE  must  be  com- 
puted to  begin  the  computation  of  NYGRID. 

i 

r 

STORE  (100,  10) 

Array  defined  in  PLOTF  and  used  to  store  the  computed  out- 
puts which  will  be  plotted  when  LF  = NLF.  For  example,  the 
second  output  for  the  fourth  frequency  is  stored  in  STORE  (4,  2). 

! 

TENPC 

A factor  added  and  subtracted  from  the  computed  maximum  and 
minimum  output  for  a MODE  1 plot.  TENPC  centers  the  plot 
on  the  space  assigned  to  it. 

1 

TTLPT 

A title  supplied  by  the  main  program  for  each  plot  if 

MODE  = 2. 

1 

XF 

The  floating  point  form  of  NF  when  NF  has  been  adjusted  to 
the  number  we  think  it  should  be  (used  when  MODE  = 1). 

1 

XMAG 

Magnitude  of  the  output  OU. 

1 

‘ 1 
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PLOTF  (Continued) 


XMAX  (10) 

XTVIIN  (10) 

XNINT 

XPHS 

YD 


The  maximuni  output  for  a mode  2 plot.  If  IWHICH  = 3, 

XMAX  (6)  through  XMAX  (10)  contains  the  phase  maximums  of 
up  to  5 outputs.  XMAX  is  computed  by  PLOTF  if  MODE  = 1; 
otherwise  it  comes  from  the  main  program. 

The  minimum  output  for  a mode  2 plot.  If  IWHICH  = 3, 

XMIN  (6)  through  XMIN  (10)  contains  the  phase  maximums  of 
up  to  5 outputs.  XMIN  is  computed  by  PLOTF  if  MODE  = 1; 
otherwise  it  comes  from  the  main  program. 

The  floating-point  form  of  NINT . 

Phase  of  the  output  OU. 

The  number  of  frequency  grid  lines  when  MODE  = 2.  If  YD 
exceeds  10,  the  labeling  will  not  be  good,  so  in  this  case  the 
number  of  grid  lines  is  decreased. 


PLOTF  calls  subroutines  EQUOUT,  OUT,  CAMRAV, 
APLOTV,  LINEV,  SMXYV,  ENDPLT. 


GRIDIV,  PRINTV,  APRNTV, 


POLAR  (NC,  AC,  FREQ,  LF,  NODMAX,  NOUT,  NEQU,  OUTPUT,  JS,  LEQU,  OU, 
NCUR,  NODE,  JCOMP,  NTYPE,  NAME,  INODE,  B,  H,  G,  F,  E,  V,  D,  Y) 

POLAR  is  called  by  the  main  program  when  a nominal  or  special  solution  has 
been  requested  and  when  the  main  program  has  solved  the  impedance  matrix  and 
obtained  node  voltages.  POLAR  computes  the  node  voltages,  branch  currents,  and 
special  equation  outputs  requested  by  the  engineer.  If  node  voltages  for  nodes  dropped 
on  conversion  to  the  equivalent  circuit  were  requested,  POLAR  prints  a list  of  these 
nodes.  In  the  same  manner,  if  branch  currents  through  circuit  elements  that  were 
dropped  on  conversion  to  the  equivalent  circuit  or  dependent  or  independent  sources 
are  requested  by  the  engineer,  POLAR  prints  a list  of  the  branch  currents  which  it  can't 
compute. 

If  the  equation  requested  by  the  engineer  for  a functional  output  was  not  supplied 
in  EQUOUT  by  the  engineer,  POLAR  prints  a diagnostic. 

AC  AC  equivalent  circuit. 

COMP  (201)  The  number  of  the  circuit  element.  (This  is  not  used  after  we 

have  finished  with  subroutine  EXCH. ) 


FREQ  (1,  1) 

ICOMP 

ICR 

INODE  (400) 

I PR 

JCOMP 


Frequency  at  which  the  circuit  was  solved. 

The  array  of  circuit  element  numbers  in  the  equivalent  circuit. 

0 or  1.  ICR  = 1 if  branch  currents  are  encountered  which  can't 
be  computed  because  the  circuit  element  has  been  dropped  on 
conversion  to  the  equivalent  circuit. 

The  node  of  interest  in  the  equivalent  circuit.  See  appendix  on 
Network  Storage. 

0 or  1.  1 if  we  encoimtered  node  voltages  which  we  couldn't 
compute  because  they  had  been  dropped  on  conversion  to  the 
equivalent  circuit;  0 otherwise. 

Total  number  of  circuit  elements  supplied  by  the  engineer. 


JS 


Number  of  functional  outputs  requested. 


K 0 or  1 depending  on  whether  we  have  any  branch  current  requests 

(K  = 1)  or  not  (K  = 0),  and  whether  or  not  we  have  already 
printed  our  heading. 

LEQU  (3,  1)  The  functional  output  equation  number  (in  LEQU  (1,  I)  and  the 

title  (in  LEQU  (2, 1)  and  LEQU  (3,1)). 


LF 


The  index  of  FREQ . 


NAME  (201)  Up  to  4 alphanumeric  characters  to  title  a circuit  element. 

NC  Either  DC  or  AC  depending  on  the  type  of  circuit. 
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POLAR  (Continued) 

1 

NCUR  (1) 

An  array  of  0 or  1 to  indicate  the  circuit  element  number  through 
which  we  want  the  branch  current.  If  NCUR  (5)  = 1 , then  we 
want  the  branch  current  through  circuit  element  5. 

r 

NEQU  (1) 

An  array  of  node  equivalences.  If  node  numbers  were  changed  on 
conversion  to  the  equivalent  circuit,  say  node  7 was  changed  to 
node  4,  then  NEQU  (8)  = 4. 

NODE  (2,  1) 

The  primary  (NODE  (1,  I))  and  secondary  (NODE  (2,  I))  nodes 
between  which  a circuit  element  lies  in  the  read-in  circuit. 

A 

NOD  MAX 

The  maximum  number  of  nodes  in  the  read-in  circuit. 

NOUT  (1,  1) 

An  array  of  0 or  I's.  If  NOUT  (5)  = 1,  then  we  want  the  output 
voltage  from  node  5. 

NTYPE  (201) 

The  type  of  circuit  element,  either  R,  L,  C,  A,  Z,  E,  V,  D,  I, 

F,  G,  B,  or  H. 

1 

OU  (2) 

Temporary  storage  for  the  real  output,  OU  (1)  and  imaginary 

t 

output,  OU  (2). 

! 

OUTPUT  (2,  1) 

Contain  the  node  voltages.  OUTPUT  (1,  3)  is  the  real  part  of 
node  2 output  voltage.  OUTPUT  (2,  3)  is  the  imaginary  part  of 
node  2 output  voltage. 

i 

XMAG 

Magnitude  of  the  output. 

XPHS 

Phase  of  the  output. 

i 
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READFQ  (NLF,  FREQ.  BLANK) 


I 

1 

1 

I 

i 

I 


i 

I 

i 

I 

1 


READFQ  is  called  by  the  main  program  and  subroutine  SPECIN  (when  a special 
solution  is  requested)  and  reads  in  or  computes  the  frequencies  for  AC  nominal, 
special,  or  frequency-plot  output  requests.  Three  types  of  frequency  constructions 
can  be  requested: 

1)  A list  of  frequencies  supplied  by  the  engineer; 

2)  A linear  construction  of  frequencies  over  a range  with  equal  increments; 

3)  A logarithmic  group  of  frequencies  constructed  between  a range  of  decades 
with  a choice  as  to  the  number  of  points  wanted  in  each  decade. 

A combination  of  types  1 and  2 or  1 and  3 can  be  selected,  and  READFQ  will 
order  the  frequencies  from  least  to  greatest. 

BLANK  Alphanumeric  character  defined  as  a blank  in  the  main  program. 

CHECK  A read-in  character  that  should  be  blank  if  the  data  cards  are  in 

the  correct  order. 

DUMB  (100)  Array  used  for  sorting  the  frequencies  when  both  read-in  and 
constructed  values  are  used. 

Array  of  different  frequencies  to  be  supplied  to  the  calling  pro- 
gram on  return. 

Index  of  the  DUMB  array  when  we  are  ordering  the  frequencies  in 
a combination  type  of  frequency. 

Read-in  variable  controlling  whether  the  frequencies  are  to  be 
read  in  (1),  constructed  (2),  or  a combination  of  both  (3). 

Index  of  the  DUMB  array  when  we  are  orderir^  the  frequencies 
in  a combination  type  of  frequency. 

Index  of  the  FREQ  array  when  we  are  ordering  the  frequencies  in 
a combination  type  of  frequency. 

0 if  this  is  a logarithmic  construction  of  frequencies;  1 if  a linear 
construction. 

0,  1 , 2,  or  3 and  is  read  in  when  LL  = 0.  It  controls  the  number 
of  intervals  minus  1 within  the  decades  we  will  compute. 

LTYPE  = 0 if  LL  = 1;  1 if  decades  only  are  to  be  used;  2 if  one 
intermediate  point  at  3 is  desired;  and  3 if  two  points,  at  2 and  5, 
are  desired. 

Number  of  frequencies  not  constructed  in  a combination  type 
because  the  final  decade  was  reached  before  NLFSP  frequencies 
were  constructed. 


FREQ  (1,  1) 

I 

IT  123 

J 

K 

LL 

LTYPE 

NDIF 
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READFQ  ( C ontinued) 


NLF  Total  number  of  frequencies  to  be  used  in  solving  the  circuit. 

NLF2  Number  of  frequencies  to  be  constructed  if  this  is  a combination  of 

the  two  types  of  frequency. 

NLF3  First  read-in  frequency  index  in  the  DUMB  array  when  a combina- 

tion of  the  two  types  of  frequencies  is  called  for. 

NLFSP  Number  of  frequencies  to  be  read  in  if  this  is  a combination  of  the 

two  types  of  frequency. 

NSAV  A temporary  storage  location  for  NLF  so  we  can  use  the  same  part 

of  the  program  for  logarithmic  and  combination-type  frequency 
construction. 

START  The  starting  frequency  value  for  a constructed  set  of  frequencies. 

STEP  The  step  size  if  LT.  - 1 (a  linear  construction),  the  final  decade 

value  if  LL  = 0 

READFQ  calls  no  subroutines. 
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RKCTAN  (RV,  A,  B.  VAL) 


I 

I 

I 

1 

I 

I 


I 

I 

I 

I 

1 


RECTAN  is  called  when  a Monte  Carlo  solution  is  requested  and  a circuit 
element  must  be  varied  according  to  a rectangular  distribution.  RECTAN  is  entered 
once  for  each  such  circuit  element  unless  the  element  is  an  A or  Z type,  in  which 
case  RECTAN  is  entered  tw'ice. 

A Lower  boundary  of  distribution. 

B Upper  boundary  of  distribution. 

RV  A uniformly  distributed  random  variable  from  the  main  program. 

VAL  The  random  variable  obtained  from  subroutine  RECTAN. 

RECTAN  calls  no  subroutines. 


1 


SKNSPR  (NODMAX,  NOUT,  ISN,  NTYPE,  NAME,  NV,  SV,  IWHICH,  JSl,  JS,  LEQU, 
OUTPUT,  NODE,  NOD,  NEQU,  XLOW,  XfflGH,  NUM,  FREQ,  NCUR, 
JCOMP,  INODE,  XNOM) 

SENSPR  is  called  by  the  main  program  when  a Type  Card  has  specified  a mode  1 
sensitivity  output.  It  is  called  once  for  each  circuit  element  that  is  varied  in  a 
sensitivity  solution  for  each  frequency.  The  main  program  obtains,  with  nominal 
values,  a low  value,  and  a high  value  node  voltages  before  calling  SENSPR.  Thus  any 
entry  to  SENSPR  will  be  with  three  sets  of  solutions  for  the  variation  of  one  circuit 
element  in  the  equivalent  circuit.  SENSPR  tests  to  see  what  type  of  output  is  required 
and  stores  the  computed  output  in  a printing  array  as  it  is  computed.  When  all  the 
outputs  for  this  circuit  element,  or  when  eight  outputs  (to  make  up  a line  of  printing) 
have  been  calculated,  SENSPR  prints  a line  and,  if  the  outputs  have  not  all  been 
printed,  continues  computing,  storing,  and  printing.  No  attempt  is  made  to  sort  or 
rearrange  the  outputs . They  are  printed  in  the  order  of  their  selection  on  the  data 
cards  by  the  engineer,  with  node  voltages,  then  branch  currents,  then  special 
functions . 

The  circuit  elements  are  varied  in  the  order  of  their  appearance  in  the  data 
deck  with  elements  whose  numbers  (col.  1-3)  were  assigned  by  engineers  being  put 
in  that  number's  position  in  the  array. 


SENSPR  returns  to  the  main  program  so  that  the  next  circuit  element  in  the 
equivalent  circuit  can  be  varied. 


FREQ 

Frequency  at  which  the  circuit  was  solved. 

11 

I - 1,  used  in  NEQU  because  zero  indices  are  not  allowed. 

IPR 

First  index  of  the  SV  array;  when  IPR  = 8,  we  are  ready  to 
print  since  we  can  only  get  8 outputs  across  a line  of 

130  characters. 

1 

IPRl 

Index  which  addresses  the  first  part  of  the  labeling  array  NV. 

i 

IPR2 

Index  which  addresses  the  second  part  of  the  labeling  array  NV. 

(NV  is  single  dimensioned,  so  it  has  2*IPR  entries. ) 

1 

ISN 

Circuit  element  number  which  we  have  varied  to  obtain  these 
solutions  to  the  circuit. 

1 

IWHICH 

1,  2,  or  3 depending  on  whether  we  have  a DC  circmt 
(IWHICH  = 1);  an  AC  circuit  where  we  want  only  the  magnitude 
(IWHICH  = 1);  the  phase  (IWHICH  = 2);  or  both  magnitude  and 
phase  (IWHICH  = 3). 

1 

JS 

Number  of  special  output  functions  requested  by  the  engineer 
(JS  - JSl). 

t 

1 

JSl 

Nximber  of  special  output  functions  requested  by  the  engineer. 

K 

Index  used  to  get  the  nominal  and  low  values  from  the  COMP 

} 

array. 
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I-KQU  (.-5.  1) 

NAME  (1) 
NKQU  (1) 

NOD 

NODE  (2,  1) 

NODMAX 
NOUT  (1,1) 

NTYPE  (1) 

NUM  (1) 

NV  (2) 


OU  (2) 

OUTPUT  (2,  1) 
SV  (8,  3) 


The  array  of  equation  numbers  (LEQU  (1,  I))  and  labels  for  the 
special  output  functions  resulting  from  solving  these  equations 
(LEQU  (2,  1))  and  (LEQU  (3,  1)). 

The  four  alphanumeric  characters  used  for  labeling  a component. 

An  array  of  equivalent  node  names.  The  index  minus  1 of 
NEQU  (I)  is  the  read  in  node  number  while  the  value  of  NEQU  (I) 
is  the  equivalent  circuit's  equivalent  node  number. 

Alphanumeric  word  "NODE"  used  to  label  node  voltage  output. 

The  pair  of  nodes  between  which  a circuit  element  lies  between. 
(Needed  for  branch  current  computation. ) 

Maximum  number  of  nodes  in  the  read-in  circuit. 

An  array  of  0 and  1 for  node  voltage  selection.  If  NOUT  (5,  1)  = 1 
then  the  output  voltage  off  of  node  5 is  desired. 

Letter  designating  the  type  of  circuit  element.  Used  for  labeling 
in  SENSPR. 

Array  of  numbers  from  1 to  50  used  to  label  node  voltage. 

Contains  a label  for  the  outputs.  For  node  voltage  it  contains 
"NODE"  in  NV  (1)  and  the  number  of  the  node  in  NV  (2).  For 
branch  current  it  contains  NTYPE  (ISN)  and  NAME  (ISN), 
respectively;  and  for  special  functions  it  contains  LEQU  (2,  I) 
and  LEQU  (3,  I),  respectively. 

Temporary  storage  of  the  real  (OU  (1))  and  imaginary  (OU  (2)) 
parts  of  an  output. 

High  values  of  the  node  voltages . Output  ( 1 , 1)  = the  real  part 
of  the  voltage  of  node  I-l. 

Used  to  save  three  lines  of  output  printing.  The  line  of  low  out- 
puts (outputs  resulting  from  putting  the  low  value  of  the  circuit 
element  in  the  circuit) , the  line  of  high  outputs , and  a bottom 
line  of  the  low  minus  the  high  divided  by  the  nominal  output  make 
up  these  three  lines.  Up  to  eight  outputs  per  line  can  be  printed. 


XHIGH  (1)  The  high  value  of  the  circuit  element  used  to  obtain  the  node 

voltage  in  OUTPUT  (I). 

XLOW  (1)  The  low  value  of  the  circuit  element,  used  to  obtain  the  node 

voltages  in  COMP  (101  + I). 

SENSPR  calls  subroutines  OUT,  EQUOUT,  BRANCH. 
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SENSPl  (ISN,  JCOMP,  NTYPE,  NAME,  NOD,  NUM,  NC,  DC,  LF,  NODMAX, 

NOUT,  NEQU,  LEQU,  IWHICH,  COMP,  OUTPUT,  JS,  ICOMP,  IVALUE, 
XNOM,  XHIGH,  XLOW,  NEXIT,  FREQ,  TITLE,  TYPE,  SUBTLE,  NCUR, 
NODE,  INODE,  STORE) 


SENSPl  is  called  by  the  main  program  when  a node  2 sensitivity  output  is 
requested,  and  after  the  main  program  has  completed  a nominal,  low,  and  high  value 
solution  for  a particular  circuit  element.  SENSPl  stores  the  circuit  element's 
identification  as  well  as  computing  up  to  three  outputs  requested  by  the  engineer; 
computes  the  quotient  of  the  high  value  minus  the  low  value  divided  by  the  nominal 
value  when  all  the  circuit  elements  in  the  equivalent  circuit  which  should  have  been 
varied  have  been  varied;  and  prints  a headii^  describing  the  analysis  it  has  done  and 
orders  the  quotients  in  STORE  for  one  output.  These  are  printed  out  and  high  and 
low  worst  cases  are  prepared  by  supplying  the  low  or  high  value  of  a particular  circuit 
element's  range,  depending  on  whether  the  quotient  was  positive  or  negative.  After 
computing  the  two  worst  cases  (by  returning  to  the  main  program)  and  printing  them, 
the  next  output  requested  and  stored  is  ordered  and  printed.  When  all  outputs  and 
their  associated  worst  cases  have  been  computed  and  printed,  a variable  is  set  to 
flag  this  completion  and  SENSPl  returns  to  the  main  program. 


COMP 


CUR  (2) 

DC 

FREQ  (1) 
ICOM  (200) 


ICOMP  (1) 
INODE  (1) 

IS 


ISN 


Circuit  element's  number  in  the  read-in  circuit,  but  now  used 
to  store  output  voltages.  COMP  (I),  1=1, 100  contains 
nominal  values.  COMP  (I),  I = 101,  200  contains  low  value 
solutions . 

Branch  current  output. 

DC  equivalent  circuit. 

Frequency  at  which  the  equivalent  circuit  is  being  solved. 

Circuit  element's  number  or  index  in  the  read-in  circuit. 

Since  not  all  circuit  elements  are  varied  in  a sensitivity 
analysis,  and  since  their  order  is  disarranged  for  each 
output,  computer  time  is  considerably  lessened  by  saving 
circuit  element  indices  in  ICOM  arranged  in  the  order  of  the 
NWORST  array,  so  that  when  the  circuit  is  being  prepared 
for  a worst  case  analysis,  it  can  be  done  so  with  a minimum 
of  programming  effort. 

Circuit  element  numbers  in  the  equivalent  circuit. 

One  of  the  nodes  in  the  equivalent  circuit.  See  appendix  on 
Network  Storage. 

Index  of  the  circuit  element  variation.  While  ISN  steps 
from  1 to  JCOMP,  some  circuit  elements  will  not  be  varied 
in  a sensitivity  analysis  and  so  ISN  will  skip  values.  IS  does 
not,  and  serves  as  fte  middle  index  of  the  STORE  array  as 
well  as  an  index  in  the  ICOM,  NAME  CP,  and  NWORST  arrays. 

Circuit  element  number  (of  the  read-in  circuit)  of  the  element 
currently  being  varied. 
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SKNSP1  (Continued) 

U'ALUE  (2,  1)  Present  value  of  the  circuit  element  in  the  equivalent  circuit. 

IVALUE  (1,  I)  contains  the  real  part,  IVALUE  (2,  I)  the 
imaginary  part. 

IWHICH  1,  2 or  3 depending  on  whether  we  want  magnitude,  phase,  or 

both,  respectively,  when  computing  and  printing  the  outputs. 

J An  index  used  to  compute  the  index  in  the  COMP  array  from  I 

to  L.  J is  the  middle  index  when  sorting  and  printing  STORE 
array.  It  is  also  the  index  which  controls  the  circuit  element. 

JC  Two  times  the  number  of  circuit  elements  in  the  equivalent 

circuit. 

JCOMP  Total  number  of  circuit  elements  in  the  engineer's  read-in 

circuit. 


( 

f 


JS  Total  functional  outputs  requested  for  this  solution  request. 

K A number  from  1 to  6,  used  to  store  nominal  (K  = 1 or  4), 

low  (K  = 2 or  5),  and  high  (K  = 3 or  6)  output  solutions  in 
STORE.  K = 4,  5 and  6 is  only  used  for  phase  storage  when 
IWHICH  =3.  K is  also  used  as  the  circuit  element's  number 
when  setting  up  the  worst  case.  K is  also  used  to  obtain 
output  in  worst  case  analysis. 

L Used  to  obtain  nominal  and  low  output  values  from  the  COMP 

array.  When  L = 0 we  are  obtaining  nominal  outputs;  when 
L = 100  we  are  getting  low  value  outputs.  L is  also  used  as 
the  IVALUE  index  in  setting  up  worst  cases. 

LEQ  Temporary  storage  for  LEQU  (1,  K),  the  equation  number  of 

a function^  output  request. 

LEQU  (3,  1)  The  functional  output  equation  number  (LEQU  (1,  I))  and  titles 

for  labeling  purposes  (LEQU  (2,  I)  and  LEQU  (3,  I)). 

LF  Index  of  the  frequency  array. 

M 1 through  6,  used  for  the  first  index  in  printing  the  STORE 

array. 

N The  output  index  while  printing. 

NAME  (1)  The  circuit  element's  title  used  with  NTYPE  for  labeling. 

NAME  CP  (2,  200)  The  labels  titling  the  name  of  the  circuit  elements  which  are 

varied.  NAMECP  (1,  I)  = NTYPE  (ISN)  and 
NAME  CP  (2,  I)  = NAME  (ISN). 


I 

1 
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SENSPl  (Continued) 
NAMO  (2,  5) 


NC 

NCUR  (1) 

NEQU  (1) 

NEXIT 


NO 

NOD 

NODE  (2,  1) 

NO  DM  AX 
NOUT  (1,  1) 

NTYPE  (1) 

NUM  (1) 


The  labels  used  for  titling  the  outputs  requested.  If  the  output 
is  node  voltage,  NAMO  (1,  NO)  = the  word  "NODE",  and 
NAMO  (2,  NO)  = a number  from  1 to  50.  If  branch  current  is 
the  output,  the  word  "CURRENT"  is  contained  in  the  two 
positions  of  NAMO.  If  functional  output  is  requested, 

NAMO  (1,  NO)  = LEQU  (2,  I)  and  NAMO  (1,  NO)  = LEQU(3,  I). 

Either  AC  or  DC,  depending  on  the  type  of  equivalent  circuit. 

An  array  of  O's  and  I's  to  indicate  if  a particular  circuit 
element  is  to  have  the  branch  current  going  through  it  com- 
puted as  an  output.  If  NCUR  (4)  = 1,  then  the  current  through 
circuit  element  4 is  requested  as  output. 

An  array  of  equivalent  node  names.  When  we  convert  from 
the  read-in  circuit  to  an  equivalent  circuit,  we  reniunber  the 
nodes.  Thus  if  NEQU  (6)  = 2,  node  2 is  now  in  the  equivalent 
circuit. 

-1  if  this  is  the  first  time  we  have  entered  SENSPl  for  this 
solution  request  card.  We  then  initialize  the  output  labels 
and  set  NEXIT  = 0. 

-0  when  we  are  storing  the  outputs  for  the  varying  compo- 
nents, but  are  not  yet  ready  to  print  the  ordered  outputs. 

1 when  we  have  filled  the  IVALUE  array  with  component 
values  set  to  produce  a worst  case  low  output. 

2 when  we  have  the  circuit  element's  IVALUE 's  filled  with 
worst  case  high  values.  We  now  print  both  high  and  low 
worst  case. 

5 when  the  sensitivity  analysis  is  completed. 

Index  of  the  output  we  are  on,  used  as  the  last  index  in  the 
STORE  array  and  as  an  index  in  the  NWC,  NAMO,  and 
PNOM  arrays. 

The  letters  "NODE"  used  for  labeling  node  voltage  request 
prints . 

The  primary  (NODE  (1,  I))  and  secondary  (NODE  (2,  I)) 
nodes  of  a circuit  element  in  the  read-in  circuit. 

Maximum  number  of  nodes  read  into  the  circuit. 

An  array  of  O's  and  I's.  If  NOUT  (4)  = 1,  then  we  want  to 
compute  the  output  voltage  from  node  4. 

Type  of  circuit  element,  either  A,  Z,  R,  C,  L,  V,  E,  D,  I, 

B,  H,  F,  or  G.  Used  here  for  labeling. 

A number  from  1 to  50,  used  for  labeling  node  voltage  request 
prints . 
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SENSPl  (Continued) 

mv 

NWC  (2,  5) 


NWORST  (200) 


OUTPUT  (2,  1) 

PNOM  (10) 

STORE  (6,  200,  3) 


SUBTLE  (1) 
TITLE  (1) 
TYPE  (1) 


Temporary  storage  for  NWC  (1,  N),  the  array  which  directs 
us  to  the  right  part  of  the  program  for  the  particular  type  of 
output  being  processed. 

Controls  the  type  of  output  being  considered.  Thus  NWC  saves 
program  testing  by  directing  us  to  the  part  of  the  program 
needed  to  compute  a particular  type  of  output.  If  NWC  (1,  3)  = 1 , 
then  the  3rd  output  is  a node  voltage  kind,  and  if  NWC  (1,  3)  =2 
it  is  a branch  current  output;  while  if  NWC  (1,  3)  = 3,  then 
the  third  output  would  be  a functional  output.  The  second  index 
is  a reference  to  which  one  of  that  type.  Thus  if  NWC  (2,  3)  = 6 
for  NWC  (1,  3)  = 1,  we  would  be  requesting  node  voltage  off 
the  6th  node;  NWC  (2,  3)  = 6 for  NWC  (1,  3)  = 2 would  be 
requesting  branch  current  through  the  sixth  circuit  element; 
and  NWC  (2,  3)  = 6 when  NWC  (1,  3)  = 3 would  be  requesting 
a functional  output  employing  the  sixth  equation  of  EQUOUT. 

An  array  of  I's  and  2's  set  to  flag  if  the  quotient  is  negative 
or  positive,  respectively,  for  this  output.  NWORST  is  set 
when  we  are  sorting  the  quotients  before  printing  them.  It  is 
used  to  control  how  the  circuit  elements  are  set  (whether  their 
low  or  high  value  is  to  be  used)  when  we  get  ready  to  compute 
worst  case  solutions. 

The  node  voltage  outputs.  For  example,  OUTPUT  (1,  4)  con- 
tains the  real  part  of  the  voltage  at  node  3,  while 
OUTPUT  (2,  4)  contains  the  imaginary  part  of  the  voltage  at 
node  3.  In  SENSPl,  output  contains  high  value  solutions. 

Contains  the  nominal  solutions  for  the  outputs  requested.  If 
IWHICH  = 3,  the  phase  outputs  are  in  PNOM  (J  + 5)  for 
J = 1,  JS. 

The  array  used  to  store  the  outputs  found  for  up  to  200 
variations  of  the  circuit  elements.  STORE  (1,  I,  N)  contains 
the  nominal  value  of  the  Nth  output  at  first;  then  after  all  the 
computations  have  been  made,  it  contains  the  high-minus-low 
divided  by  nominal  value  for  the  Ith  component  varied. 

STORE  (2,  I,  N)  contains  the  low  value  solutions,  and 
STORE  (3,  I,  N)  the  high  value  solutions.  If  both  magnitude 
and  phase  are  requested  for  all  the  outputs  (if  IWHICH  = 3), 
the  magnitude  solutions  are  stored  in  STORE  (1,  I,  N), 

STORE  (2,  I,  N)  and  STORE  (3,  I,  N);  while  the  phase  solu- 
tions are  stored  in  STORE  (4,  I,  N),  STORE  (5,  I,  N),  and 
STORE  (6,  I,  N). 

Subtitle  on  the  solution  request  card  of  the  main  program. 

Title  of  this  circuit  read  into  the  main  program. 

Type  of  solution  requested — in  this  case,  "SENS". 
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SENSPl  (Continued) 

WCMAXM 

WCMAXP 

WCN 

WCSAVIN 

WCSAVX 

WCX 

XHIGH  (1) 

XLOW  (1) 

XMAG 
XNOM  (1) 

XPHS 


Magnitude  of  the  output  for  a worst  case  high  solution. 

Phase  of  the  output  for  a worst  case  high  solution. 

Worst  case  low  output,  either  phase  or  magnitude,  whichever 
was  requested. 

Worst  case  low  phase  output,  saved  when  IWHICH  = 3 for 
later  printout. 

Worst  case  high  phase  output,  saved  when  IWHICH  = 3 for 
later  printout. 

Worst  case  high  output,  either  phase  or  magnitude,  whichever 
was  requested. 

High  value  of  the  circuit  element  in  the  read-in  circuit. 

Low  value  of  the  circuit  element  in  the  read-in  circuit. 
Magnitude  of  the  output. 

Nominal  value  of  the  circuit  element  in  the  read-in  circuit. 
Phase  of  the  output. 


SENSPl  calls  subroutines  OUT,  BRANCH,  EQUOUT, 


t 


1 

1 

! 

i 

j 


I 
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SOLUT  (LF,  NDEBUG,  NC,  TYPE,  SUBTLE,  ND,  JNODE,  MS,  ITYPE,  INODE, 
NODMAX,  NERROR,  V,  B,  E,  NEQU,  R,  KS,  NTYPE,  NAME,  NODE,  IP, 
NPOW,  XTIOM,  JCOMP,  II) 

SOLUT  is  called  by  the  main  program,  cnee  for  each  solution  SNAP  II  requires, 
;ifter  subroutine  ASSIGN  has  been  called.  SOLUT  tests  that  the  value  computed  in 
ASSIGN  for  a resistor  is  real.  If  a resistor  was  assigned  an  imaginary  part  through 
an  equation  in  EQUIN,  NERROR  is  set  to  1 in  ASSIGN  and  on  testing  in  SOLUT,  a 
diagnostic  is  printed,  the  circuit  is  printed,  and  the  run  stops. 

If  NERROR  = 0,  and  the  debug  print  option  is  1,  the  equivalent  circuit  is  printed 
with  the  values  for  the  circuit  elements  assigned  by  ASSIGN. 

Next,  SOLUT  tests  for  voltage  sources.  Any  voltage  source  that  has  not  already 
been  converted  to  current  (in  which  case  NSCR  (I)  = -20)  must  be  converted  to  current. 
To  do  this,  SOLUT  calls  subr'^utine  CURENT  once  for  each  voltage  source.  CURENT 
must  be  called  even  if  this  is  an  equivalent  circuit  used  many  times  before.  Since 
the  values  assigned  by  ASSIGN  are  different  for  each  solution,  the  calculation  of 
current  will  be  different,  even  if  the  topology  of  the  circuit  changes  the  first  time 
through  CURENT. 

Finally,  the  equivalent  circuit  is  again  printed  out  if  NDEBUG  = 1 and  if  this 
is  a DC  circuit  or  the  first  frequency  of  an  AC  circuit. 


B 

DUM  (201) 

E 

H 

II 

ICOMP  (400) 
INODE  (1) 

IP 

ITYPE  (1) 
IVALUE  (2,  1) 

JCOMP 
JNODE  (2,  1) 


Voltage-dependent  current  source. 

A place  holder  for  common  storage. 

AC  voltage  circuit  element. 

Voltage-dependent  voltage  source. 

A node  index  minus  1 used  to  make  the  resulting  number  equal 
to  the  correct  node  name,  since  we  use  zero  nodes,  but  cannot 
use  zero  as  an  index. 

Circuit  element  numbers  in  the  equivalent  circuit. 

One  of  the  nodes  between  which  a circuit  element  in  the  equivalent 
circuit  lies. 

Index  used  by  the  NPOW  array. 

Circuit  element  type  in  the  equivalent  circuit,  either  R,  L.  C, 

V,  E,  D,  I,  A,  Z,  B,  H,  F,  or  G. 

Value  of  the  circuit  element  in  the  equivalent  circuit.  This  is 
changed  by  CURENT  when  the  circuit  element  is  a voltage 
source. 

Number  of  circuit  elements  in  the  read-in  circuit. 

Array  which  serves  as  a table  of  contents  for  the  arrays  describ- 
ing the  equivalent  circuit.  See  appendix  on  Network  Storage. 
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Index  of  an  array  in  CURENT  used  to  store  the  resistor  values 
needed  to  convert  voltage  to  current. 

Index  of  the  frequency  we  are  using  for  this  solution. 


I 

I 


KS 

LF 

MS 

NAME  (1) 

NB 

NC 

ND 

NDl 

NDEBUG 

NEQU  (1) 

NERROR 

NF 

NODE  (2,  1) 

NODM.VX 

iDMXl 

!■«  »w  (2.  n 


0 or  1 depending  on  whether  or  not,  respectively,  we  have  a new 
equivalent  circuit.  MS  is  set  to  1 in  main  after  returning  from 
subroutine  SOLUT. 

Up  to  4 alphanumeric  characters  used  to  describe  the  circuit 
element  besides  its  type. 

The  beginning  range  of  a DO  loop  which  searches  through  one 
node's  entries  in  the  equivalent  circuit's  arrays. 

Type  of  equivalent  circuit  we  are  solving.  NC  = AC  or  DC. 

Number  of  nodes  in  the  equivalent  circuit. 

Number  of  nodes  in  the  equivalent  circuit  plus  one. 

1 or  0 depending,  respectively,  on  whether  or  not  we  want 
intermediate  debugging  prints. 

An  array  of  node  name  equivalences  between  the  read-in  circuit 
(contained  in  the  index  + 1 of  NEQU)  and  the  node  number  of  the 
equivalent  circuit.  Thus  NEQU  (5)  = 3 means  that  node  4 of  the 
read-in  circuit  is  numbered  node  3 of  the  equivalent  circuit. 

0 or  1 depending  on  whether  ASSIGN  found  that  formula-computed 
resistor  elements  were  real  or  had  imaginary  parts,  respectively. 

The  end  of  a range  of  a DO  loop.  See  NB. 

Node  numbers  between  which  this  circuit  element  lies. 

NODE  (1,  J)  is  the  primary  node,  NODE  (2,  J)  is  the  secondary 
node. 

Number  of  nodes  in  the  read-in  circuit. 

Number  of  nodes  in  the  read-in  circuit  + 1 . 

The  array  holding  resistor  codes  for  use  in  converting  voltage 
to  current.  NPOW  (1,  I)  holds  the  circuit  element  number. 

An  array  used  by  SOLUT  and  CURENT  to  flag  those  voltage 
entries  in  the  equivalent  circuit  which  have  already  been  con- 
verted to  current. 

The  circuit  element  type,  either  R,  C,  L,  A,  Z,  E,  V,  D,  1, 

H.  H.  F or  G. 
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SOLUT  (Continued) 

R Resistor-type  circuit  element. 


SUBTLE  (1)  A subtitle  labeling  the  solution  request  card  (Type  Card)  and 

used  by  SOLUT  to  label  the  printout  when  NDEBUG  = 1. 

TYPE  The  type  of  solution  we  are  calculating.  This  variable  is  used 

for  printing  a title  when  NDEBUG  = 1. 

V Voltage-source  circuit  element. 

XNOM  (1)  Nominal  values  of  the  circuit  elements. 

St)LUT  calls  subroutine  CURENT. 

r 

I 

I 

I 


I 

I 

i 
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SOLVE  (NAD,  A.,  NDIM,  N,  NER) 


SOLVE  is  called  by  subroutine  COMPUT  once  for  each  circuit  solution  required. 

SOLVE  solves  sets  of  simultaneous  linear  equations  in  double  precision  for  either 
complex  or  linear  systems.  * 

Since  double  precision  operations  on  machines  without  double  precision  hard- 
ware are  expensive,  and  since  most  circuit  admittance  matrices  have  many  zero 
elements,  SOLVE  tests  for  zero  elements  before  performing  such  operations.  The 
complex  arithmetic  is  done  in  the  program,  since  FORTRAN  IV  compilers  will  not 
do  both  double  precision  and  complex  arithmetic. 

To  reduce  round-off,  each  row  is  divided  by  its  biggest  element  so  that  all  the 
elements  are  between  -1  and  1.  This  is  a compromise  since  full  column  scaling  and 
pivoting  on  the  largest  element  would  further  reduce  roimd-off,  but  would  also  increase 
machine  time  significantly. 

A Impedance  matrix  plus  the  right  hand  sides  as  an  additional  column. 

DIV  Reciprocal  of  the  largest  element  in  the  row  If  the  matrix  is  complex, 

this  is  the  reciprocal  of  the  magnitude,  squared. 

DIVI  Imaginary  part  of  the  reciprocal  of  a complex  element  aji. 

DIVR  Real  part  of  the  reciprocal  of  a complex  element  an.  i 

II  Index  used  to  obtain  the  imaginary  part  of  the  column  elements. 

11  (I-l)  is  a DO  loop  which  computes  the  column  vector  sum  below  the 

diagonal. 

12  Index  of  the  column  entry  for  the  DC  case.  See  J2. 

IR  Index  used  to  obtain  the  real  part  of  the  column  elements. 

J2  Index  of  the  column  entry  for  the  DC  case.  We  must  skip  the  imaginary 

entries  to  save  computer  time,  so  J2  = 2*J-1  and  will  index  real  parts 
only. 

JI  Inaex  used  to  obtain  the  imaginary  part  of  the  column  elements. 

f 

JR  Index  used  to  obtain  the  real  part  of  the  column  elements . • 

K Index  used  to  go  from  the  second  to  last  row  up  to  the  first  row.  ; 

i 

K1  K + 1. 

K2  Index  of  the  column  entry  for  the  DC  case.  See  J2.  { 


♦The  method  used  is  the  Crout  method,  described  in  a paper  in  the  AIEE  i 

Transactions,  1941,  Vol.  60,  page  1235.  Further  references  are  Hildebrand’s  .1 

"Introduction  to  Numerical  Analysis"  from  McGraw-Hill,  1956  page  429.  The 
Crout  method  is  basically  the  back  solution  of  a Gauss  elimination  scheme. 

1 
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KI 

Kli 

N 

NAD 

N1 

N2 

NDIM 

NMl 

PI 

PIl 

PI2 

PR 

PRl 

PR2 

SUM 

SUMI 

SUMR 

XMAX 

XMIN 

SOLVE 


Index  used  to  obtain  the  imaginary  part  of  the  column  elements. 

Index  used  to  obtain  the  real  part  of  the  column  elements . 

Number  of  rows  in  matrix  A. 

0 if  this  is  a real-valued  matrix  (i.  e. , if  we  are  on  a DC  circuit).  Other- 
wise, if  NAD  = 1,  the  circuit  is  AC  and  the  matrix  is  complex. 

Number  of  columns  in  the  augmented  matrix. 

Twice  the  number  of  rows  in  A,  equal  to  the  number  of  real  and  imaginary 
rows. 

The  first  dimension  of  the  A array,  needed  by  the  compiler  and  sent 
from  the  COMPUT  subroutine. 

The  number  of  rows  minus  1,  used  when  we  are  computing  from  the 
second  to  last  row  to  the  first  row  to  obtain  a solution. 

Imaginary  part  of  the  product  of  aik  akj- 

Imaginary  part  of  a sum  used  to  perform  complex  arithmetic.  We  are 
computing  (x  + iy)(u  + iv)  where  (x  + iy)  = “ %j- 

Imaginary  part  of  a sum  used  to  perform  complex  arithmetic.  We  are 
computing  (x  + iy)(u  + iv)  where  (x  + iy)  = a^k  and  (u  + iv)  = akj 

Intermediate  sum  used  in  computing  the  new  elements  of  a row.  In  the 
complex  case,  it  is  the  real  part. 

Real  part  of  a sum  used  to  perform  complex  arithmetic.  We  are 
computing  (x  + iy)(u  + iv)  where  (x  + iy)  = ajk  and  (u  + iv)  = akj. 

Real  part  of  a sum  used  to  perform  complex  arithmetic.  We  are 
computing  (x  + iy)(u  + iv)  where  (x  + iy)  = a^k  and  (u  + iv)  = akj. 

Storage  variable  used  to  hold  the  sum  over  K of  ajk  akj. 

Sum  of  the  imaginary  parts  of  the  sum  over  k of  ajk  akj.  See  SUMR. 

Sum  of  the  real  parts  of  the  sum  over  k of  aik  akj,  used  in  the  case  of  an 
AC  circuit  when  we  are  accumulating  lower  half  column  sums. 

Maximum  element  for  each  row  used  to  scale  all  the  elements  of  that  row. 

Minimum  element  for  each  row  used  to  scale  all  the  elements  of  that  row. 


calls  no  subroutines. 
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SPCE  (RR,  RI,  DIST,  VR,  N,  VI,  MD)  ( 

SPCE  is  called  by  the  main  program  when  a Monte  Carlo  solution  is  required 

and  if  at  least  one  circuit  element  of  the  equivalent  circuit  must  be  varied  according 

to  a special  distribution.  SPCE  is  called  for  each  such  circuit  element. 

A Admittance-type  circuit  element. 

DIF  A variable  internal  to  SPCE,  used  to  store  the  difference  between  the 

distribution  RI  or  RR  is  close  to.  i 

DIST  (2,  1)  Cumulative  distribution  provided  by  main.  The  abscissa  is  in 
NDIST  (1,  I),  the  ordinate  in  NDIST  (2,  I). 

lOVER  Variable  used  only  by  SPCE.  It  is  0 when  computing  the  real-part  i 

specially  distributed  random  variable,  and  is  1 when  computing  the 
imaginary  part. 

MD  Number  of  pairs  of  points  in  DIST  (supplied  by  the  main  program). 

i 

N Circuit  element  type  from  the  main  program. 

RI  Uniformly  distributed  random  variable  supplied  by  main  for 

imaginary  part  of  circuit  element  if  the  element  is  A or  Z. 

RR  Uniformly  distributed  random  variable  supplied  by  main  for  real 

part  of  circuit  element. 

RV  Uniformly  distributed  random  variable  equal  to  RR  or  RI.  It  is  used 

internally  in  SPCE. 

VAL  The  specially  distributed  random  variable  which  will  equal  either  VR 

or  VI.  It  is  used  internally  in  SPCE. 

VI  The  special  distributed  random  variable  to  be  used  by  the  main  pro- 

gram for  the  imaginary  part  of  the  circuit  element  if  the  element 
is  an  A or  Z type. 

VR  Special  distributed  random  variable  to  be  used  by  the  main  program 

for  the  real  part  of  the  circuit  element. 

XM  The  slope  of  the  interpolating  points.  Used  only  internally  in  SPCE.  I 

Z Impedance-type  circuit  element.  i 

SPCE  calls  no  subroutines. 
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SPECIN(ISP,  NHOLD,  JCOMP,  TITLE,  NC,  TYPE,  SUBTLE,  XNOM,  XLOW, 
XHIGH,  NCUR,  NL,  SPEC,  NPEC,  SV,  KV,  NV,  IS,  NSPEC,  NAME, 
NTYPE,  XMCHI,  XMCLO,  NOUT,  LEQU,  JS,  NODMAX,  LF,  DC,  AC,  ALL) 

SPECIN  is  called  by  the  main  program  whenever  we  have  a special  solution 
request.  If  more  than  one  set  of  special  values  is  to  be  solved,  SPECIN  is  called  for 
each  such  set.  This  subroutine  reads  in  the  frequencies  (if  this  is  an  AC  circuit  and 
the  first  time  SPECIN  has  been  called  for  this  Type  Card),  reads  in  node  voltage, 
branch  current,  and  special  equation  output  requests  (if  this  is  the  first  time  SPECIN 
has  been  entered  for  this  Type  Card)  and  reads  in  one  set  of  special  value  designator 
cards.  SPECIN  then  prints  the  values  assigned  to  the  variables  in  the  equivalent 
circuit  and  returns  to  the  main  program. 

.AC  AC  circuit. 

.ALL  Indicates  that  all  the  node  voltages  or  branch  currents  are 

requested. 


COMP  (201) 


A dummy  array  in  SPECIN. 


DC  DC  circuit. 

ICOMP  (400)  An  array  of  all  the  component  numbers  in  the  equivalent  circuit. 

I PR  Index  of  SV,  KV,  NV  arrays.  When  IPR  = 8,  a line  is  printed. 

IS  Index  used  by  the  NPEC  and  SPEC  array. 


IVALUE 


JCOMP 


Index  of  the  particular  group  of  special  values  we  are  processing. 
The  first  time  control  enters  SPECIN  for  this  Type  Card, 

ISP  = 1. 

Values  for  the  circuit  elements  in  the  equivalent  circuit. 

The  number  of  circuit  elements  in  the  read-in  circuit. 


The  number  of  entries  in  LEQU  (i.  e.  , the  number  of  functional 
outputs  designated  on  the  Type  Card). 

Used  to  find  an  equivalent  circuit  index  corresponding  to  a 
particvilar  circuit  element. 


KV  (8,  1) 


The  other  equivalent  circuit  index,  different  from  Kl. 
Fixed-point  equivalent  array  of  SV. 


LEQU  (3,  1)  An  array  of  special  outputs. 

LF  Frequency  index. 

M A temporary  storage  for  the  NHOLD  entry. 

NAME  (201)  Up  to  4 alphanumeric  characters  used  to  describe  the  circuit 

element  besides  its  type. 


i 


SPE CIN  (Continued) 


NBR  (19) 


A temporary  storage  variable  for  a value  of  NBR. 

A holding  array  for  reading  in  individual  branch  current  requests . 


NCUR  (1) 


NHOLD  (1) 


AC  or  DC. 

An  array  of  0 or  1 indicating  branch  current  requests.  If 
NCUR  (6)  = 1,  for  example,  the  current  through  circuit  element  6 
is  desired  as  output. 

An  array  of  integers  from  0 to  9 which  directs  what  value  is  to  go 
in  a particular  component.  If  NHOLD  (8)  = 2,  then  circuit 
element  8 should  take  on  the  high  value  found  in  columns  51-60  of 
the  Circuit  Element  Card.  The  numbers  and  their  values  are 
as  follows: 

NHOLD  (I)  = 0,  nominal  value  (col.  31-40). 

NHOLD  (I)  = 1,  low  value  (col.  41-50). 

NHOLD  (I)  = 2,  high  value  (col.  51-60). 

NHOLD  (I)  = 3,  Monte  Carlo  low  (col.  61-70). 

NHOLD  (I)  = 4,  Monte  Carlo  high  (col.  71-80). 

NHOLD  (I)  = 5,  2nd  card,  col.  31-40. 

NHOLD  (1)  = 6,  2nd  card,  col.  41-50. 

NHOLD  (I)  = 7,  2nd  card,  col.  51-60. 

NHOLD  (I)  = 8,  2nd  card,  col.  61-70. 

NHOLD  (I)  = 9,  2nd  card,  col.  71-80. 

NHOLD  is  read  in  by  SPE  CIN  on  up  to  3 cards  (depending  on  the 
value  of  JCOMP).  Each  column  of  the  cards  is  an  entry  in  NHOLD. 

1 or  0 depending,  respectively,  on  whether  or  not  we  have  a 
nonlinear  circuit. 


NOD  MAX 


Number  of  nodes  in  the  read-in  circuit. 


NOUT  (1) 


NPEC  (6,  1) 


NSPEC  (201) 


NSUT  (25) 
NTYPE  (201) 


NV  (1) 


An  array  of  0 or  1 giving  the  node  voltage  requested.  If 
NOUT  (4)  = 1,  then  we  want  the  voltage  off  of  node  4. 

The  SPEC  array  in  fixed  format,  used  to  hold  the  circuit 
element's  number  in  SPEC  (1,  I). 

An  array  which  tells  the  program  if  this  circuit  element  has 
special  values  in  addition  to  those  in  columns  41-80  of  the  Circuit 
Element  Card. 

A temporary  storage  variable  for  a value  of  NSUT. 

A holding  array  for  reading  in  individual  node  voltage  requests. 

The  circuit  element  type,  either  R,  C,  L,  A,  Z,  E,  V,  D,  I, 

B,  H,  F,  or  G. 

The  code  number  from  the  NHOLD  array.  This  is  printed  out 
with  the  name  and  values  of  the  circuit  element. 
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SPKCIN  ( Continued) 


SPEC  (G,  1) 

SUBTLE  (1) 
SV  (8.  1) 

TITLE  (1) 
TYPE 
XHIGH  (1) 

XLOW  (1) 

XMCHI  (201) 

XMCLO  (201) 

XNOM  (1) 
SPECIN  calls 


The  NPEC  array  holding  special  values  for  some  of  the  circuit 
elements. 

A title  read  in  with  the  Type  Card,  used  to  title  this  solution. 

Used  to  accumulate  a line  of  circuit  element  values  for  printii^. 
Only  8 circuit  element's  values  can  be  printed  in  a line.  SV  (I,  3) 
contains  the  numerical  value,  while  SV  (I,  1)  = the  type  of  circuit 
element  and  SV  (I,  2)  = its  name.  The  arguments  SV  (I,  1)  and 
SV  (I,  2)  are  addressed  as  KV. 

Title  of  the  circuit,  used  for  labeling. 

The  letters  "SPCL",  used  for  labeling. 

High  value  of  the  circuit  element  for  sensitivity  analysis,  or  a 
special  value. 

Low  value  of  the  circuit  element  for  sensitivity  analysis,  or  a 
special  value. 

Standard  deviation  of  the  Monte  Carlo  distribution  (if  a normal 
or  lognormal  distribution  is  required),  or  a special  value. 

Mean  of  the  Monte  Carlo  distribution  (if  a normal  or  lognormal 
distribution  is  required),  or  a special  value. 

Nominal  values  of  the  circuit  elements. 

no  subroutines. 
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STAT  (ISTAT,  JS,  IW,  HIHIST,  TITLE,  NC,  TYPE,  SUBTLE,  LEQU,  IR,  NORV, 

NSTR,  OUTPUT,  XLfflST,  FREQ)  

STAT  is  called  by  the  main  program,  once  for  each  solution  required,  when  a 
Monte  Carlo  solution  is  needed.  Thus  if  NORV  = 500,  subroutine  STAT  is  called 
500  times. 

The  first  time  STAT  is  called,  ISTAT  is  zero,  which  flags  the  program  to  set 
arrays  to  zero  and  compute  the  bin  size.  Each  output  can  be  thought  of  falling  into 
a particular  box,  and  each  box  size  is  dependent  on  the  size  of  BIN.  After 
initialization,  ISTAT  is  set  to  1 and  STAT  begins  the  part  of  the  program  common 
to  all  entries  in  STAT. 

The  output  is  computed  for  up  to  five  outputs  and  the  resulting  number  is 
compared  to  the  limits  read  in  or  computed  in  the  main  program.  If  the  output  falls 
outside  these  limits,  it  is  printed  as  an  out  of  limit  output  and  a counter  is  increased 
by  one.  If  any  output's  counter  exceeds  20  percent  of  NORV  (the  number  of  solutions 
for  this  Monte  Carlo  request),  a diagnostic  and  the  incomplete  histogram  of  all  the 
outputs  are  printed. 

If  the  outputs  are  within  the  allowable  range,  STAT  adds  their  values  into  the 
sums  used  for  computing  the  mean  and  variance,  and  then  finds  what  box  they  belong 
in  for  the  histogram  plot  routine. 

For  example,  if  an  output's  value  is  26.  3 and  the  bin  size  (computed  from  the 
range  divided  by  100)  is  5.1,  then  26.  3 would  fall  in  the  fourth  box,  if  the  boxes 
were  though  of  as  ranging  as  follows; 


Box  1; 

10.  2 to  15.3 

Box  2: 

15.  3 to  20.4 

Box  3; 

20.4  to  25.5 

Box  4: 

25.5  to  30.6 

etc. 

Thus  the  entry  in  NBOX  (4,  I)  would  be  increased  by  one  if  this  were  the  Ith  output. 

When  IR  = NORV  (we  have  done  the  last  solution  required),  then  we  compute 
the  mean,  variance  standard  deviation,  and  the  values  between  which  the  median 
lies.  We  also  rescale  the  bin  size  so  it  is  either  one-half  a standard  deviation  (if 
NORV  is  less  than  300),  or  one-fourth  a standard  deviation  (if  NORV  is  greater 
than  300). 

Finally,  the  labeling  used  by  subroutine  HIST  is  constructed  and  subroutine 
HIST  is  called.  On  return  from  HIST,  the  mean,  median,  variance,  and  standard 
deviation  are  printed  and  the  intervals  of  the  12  (if  NORV  is  less  than  300)  or  24 
(if  NORV  is  greater  than  300)  boxes  are  printed.  Control  returns  to  the  main  program. 

Temporary  storage. 

Size  of  the  division  of  the  boxes  for  the  histogram  before  it  is 
scaled. 

Temporary  storage. 


BEG 
BIN  (5) 

DIF 


OBSC. 
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ST  AT  (Continued) 


FINISH 


FREQ 


Incremented  value  used  to  print  the  intervals  in  JBOX  after 
calling  HIST. 

Frequency  this  equivalent  circuit  is  being  solved  for. 


FS 


Temporary  storage  used  to  compute  increments  of  standard 
deviation  fractions. 


HIHIST  (5)  Upper  limit  of  the  outputs.  This  is  read  into  main  or  computed 

as  20  percent  of  the  nominal  value  of  the  output  if  the  engineer 
didn't  supply  limits. 

IK  Index  of  JBOX. 


IM 

m 

1ST  AT 

IT 

IW  (5) 

J3 

JBOX  (50) 

JS 

JT 

K 

KB 

KM 

KS 

LABEL  (25) 
LEQU  (3,  5) 

LIG3 

M 


Index  of  the  SINGO  array. 

Index  of  the  solution  being  sought.  IR  = NORV  when  we  are 
finished. 

0 the  first  time  STAT  is  entered  for  this  Type  Card.  ISTAT  = 1 
thereafter.  This  flags  the  program  to  initialize  the  arrays. 

Used  to  compute  entries  in  JBOX. 

Array  of  1 or  2's  for  AC  outputs.  If  IW  (I)  = 1,  the  Ith  output  will 
be  in  magnitude.  If  IW  (I)  = 2 the  Ith  output  will  be  in  phase. 

Counter  in  printing  the  intervals  in  JBOX  after  callii^  HIST. 

The  NBOX  array  is  compressed  into  the  JBOX  array  for  use  by 
HIST.  JBOX  is  divided  along  fractions  of  standard  deviations. 

Number  of  outputs  requested. 

Temporary  counter. 

Counter  in  printing  the  intervals  in  JBOX  after  calling  HIST. 
Starting  value  in  a DO  loop  to  compute  NBETWH. 

Temporary  storage. 

Used  to  scale  NBOX  to  JBOX. 

Labeling  array  needed  to  label  the  bottom  of  the  histogram  plot. 

The  functional  output  equation  number  (LEQU  (1,  I))  and  titles  for 
labeling  purposes  (LEQU  (2,  I)  and  LEQU  (3,  1)). 

2*NF  used  to  compute  JBOX  entries. 

Used  to  print  the  list  of  the  number  in  each  JBOX. 
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STAT  (Continued) 

1 

M2 

Temporary  storage. 

NB 

Used  to  compute  entries  in  JBOX. 

• 

NBETWL  (5) 

Number  of  outputs  above  the  lower  limit  but  below  the  minus-three 
sigma  point. 

t 

NBETWH 

Number  of  outputs  below  the  upper  limit  but  above  the  plus-three 
sigma  point. 

r 

i 

\ 

t 

NBOX  (100,  5) 

An  array  of  divisions  into  which  we  add  one  if  the  output  falls  into 
its  range. 

i 

NC 

Either  AC  or  DC,  depending  on  the  type  of  equivalent  circuit. 

NF 

The  part  of  the  calculation  of  the  variable  used  for  the  upper  limit 
of  a DO  loop. 

i 

NFl 

The  upper  range  of  a DO  loop. 

I 

! 

NMl 

Temporary  storage. 

1 

NM2 

Temporary  storage. 

.1 

NORV 

The  total  munber  of  solutions  we  will  do  for  this  Type  Card. 

i 

NOUND  (5) 

A counter  of  out-of-range  outputs. 

1 

NOUNHI  (5) 

Number  of  outputs  above  the  upper  limit. 

1 

NOUNLO  (5) 

Number  of  outputs  below  the  lower  limit. 

1 

1 

NOW 

Temporary  counter. 

NR  (5) 

Number  of  outputs  within  the  limits  used  in  histogram. 

NRM 

NR/2  used  in  computing  the  median. 

NRT 

Temporary  storage  for  NR. 

1 

NSUM 

Temporary  storage  when  computing  the  median. 

NT 

Temporary  coimter. 

.! 

NTl 

2*NF  used  to  compute  JBOX  entries. 

'1 

OUTPUT 

Node  voltage  output. 

1 

RANGE 

Expected  range  of  the  output,  computed  from  the  limits. 

ii 

RN 

Floating-point  value  of  NR. 

4 


STAT  (Continued) 


RNl 

SIG 

SIGNO  (7) 

SQUAR 
STAN  (5) 

STR 

SUBTLE  (35) 
THIS 

TITLE  (40) 
TYPE 
VARI  (5) 
XINT 

XLHIST  (5) 
XMEAN  (5) 

XMED  (2,  5) 
XP  (5) 


RN-1. 

Fractional  part  of  the  standard  deviation,  either  one-half  or 
one-fourth  of  it. 

Labeling  array  used  by  HIST  which  contains  the  mean  and  multiples 
of  the  standard  deviation. 

Temporary  storage  while  we  are  computing  the  variance. 

Standard  deviation  of  the  outputs. 

The  starting  value  for  the  random  number  generator,  printed  in 
STAT. 

Subtitle  on  the  solution  request  card  of  the  main  program. 

A sum  when  we  are  finding  the  right  NBOX  entry. 

Title  of  this  circuit  read  into  the  main  program. 

The  type  of  solution  requested.  In  this  case  a "MONT"  one. 

Variance  of  the  outputs. 

Sum  when  scaling  NBOX. 

Lower  bounds  for  the  outputs.  See  HIHIST. 

The  sum  of  the  outputs  as  they  are  computed;  used  to  hold  the 
means  of  the  outputs  when  IR  = NORV. 

The  two  box  limits  on  the  medians  of  the  outputs. 

The  phase  of  the  output. 


STAT  calls  subroutines  EQUOUT,  OUT,  HIST. 


3-69 


4.  LOGICAL  FLOW  DIAGRAMS 


Logical  flow  diagrams  for  the  various  subroutines  appear  in  the  following 
sequence: 


ACEQ 

EQUIN 

READFQ 

ASSIGN 

EQUOUT 

RECTAN 

BRANCH 

EXCH 

SENSPR 

CARD 

HIST 

SENSPl 

CHECK 

INOUT 

SOLUT 

COMPUT 

NOUN 

SOLVE 

CONECT 

NORM 

SPCE 

CURENT 

OUT 

SPEC  IN 

CURR 

PLOTF 

STAT 

DCEQ 

POLAR 

I 

I 

( 

I 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963  A 


ACEQ-2 


Construct  INODE, 
INODE  arrays  from 
NSCR. 


Fill  ITYPE  array. 


Count  number  of  nodes 
now  In  equivalent 
circuit  = NDR. 


1=1,  NODMXl 


( 


Yes 


Print  error 
codes . 


COMPUT-  2 


COMPUT-5 


CONECT-I 


Set  indices  of  two  arrays 
containing  separated 
nodes  to  zero.  Compute 
number  of  non-ground 
nodes,  set  NOW  = 0. 


^ First  time  > 

entered  CONECT  for  this 
equivalent  circuit, 

V i.e.,N  = l ^ 


Find  first  node  in 
equivalent  circuit. 


Assign  value  of  node  N to 
NFIRST  (N  sent  from 
SOLUT  as  node  we  want 
to  start  with) . 


Put  this  node  in 
NFIRST. 


Increment  NS  by  one  and 
put  NFIRST  in  NS  array, 
set  NOW  = NFIRST  to  get 
going  on  tracing  links  out 
of  NFIRST. 


r ^ 

Page 


^ V 


1 

I 


start 


CURENT-I 


Find  two  entries  in 
IVALUE  array,  com- 
pute current  from 
resistor  in  NSAVE. 


CURENT  - 2 


Find  NPOW  entry  for 
this  voltage.  Compute 
resistor  and  current. 


Page  2-II 


r node  we  > 
are  looking 
at  (to  add  resistor  to) 
the  last  one  of  the 
N.  equivalent  y 
circuit 


Increase  JC  by  1;  compute 
NB,  the  position  inIVALUE 
array  where  the  resistor 
entry  will  go;  increase  the 
link  count  of  that  node  by 


Move  INODE,  ITYPE, 
ICOMP,  IVALUE 
arrays  down  to  make 
room  for  new 
resistor. 


Also  move  INODE 
arrays  and  incre- 
ment the  link 
counter  array. 

JC  = JC  + 1. 


yT  Have  we 
“^one  this  for  botlr 
nodes  that  the  volt 
v age  lies  . 
between 


i-aBe 

I 


DCEQ- 4 


I 


No 


EXCH-2 


start  J 

V J 


HIST -I 


■ 


Define  data  words 
which  make  symbols 
and  bars  on 
Histogram  plot. 


Save  array  to 
be  plotted  so  we 
can  scale  it. 


Find  maximum  entry 
in  saved  array  so  we 
know  height  of  plot. 


Divide  max.  by  2. 
Keep  track  of 
this  in  K2. 


Divide  all  elements 
of  KCT  array 
by  2. 


Construct  ordinate 
axis  labels;  scale 
values  according  to 
how  KCT  was  scaled. 


Write 

title. 


1 


HIST- 2 


IP  ^ 1,50 
(50  lines  of  print) 


Page  2- 1 1 


Page  3 


I = 1,  12 


Should  we  print 
entry  for  this  bar 


Completed 
a line 


12-bar 

Histogram 


Same  as  for  12-bar 
Histogram  except 
I 1, 24  and 
spaces  smaller. 


Did  we  print 
bar  to  the  left 
of  this  one 


Fill  print 
space. 

i No 

Is  this 
top  of  bar 


Page  3 


Put  different 
edge  on 
previous  bar. 


Yes 


Yes 


No 


Print  24-bar 
bottom  line. 


1 


Print  2 lines,  labeling 
and  giving  values  of 
sigma  and  mean  points. 


Print  12-bar 
bottom  line. 


I 


Return 


start 


INOUT 


Print  C.  K.  's  values  as 
they  ap|X!ar  on  C.  K. 
card. 


Read  type  cards  until 
NKXT  or  FINI  comes 
up. 


•Arrange  C.  K.  's  in 
t'.K.  number  order. 
Call  KXCH  for  each 
C.  K. 


^ CK  ^ 

one  t . K. 

number  used  more 
than  once 


Print  diagnostic. 


TNo 


Put  I's  in  all 
NOl'T  arrays. 


Huad  in  s[)Ocial 
iunction  formula 
numixTs  ami  titic.s. 


Read  BRANCH 
request  card. 

current 

Same  procedu 
NCUR  array  a 
NOUT  array, 

re  for 
s for 
above. 

^ Any  > 
special  function 
outputs 


\ . 

— 1 

Fill  IVALUK  array 

with  nominal  values. 

■ 

IU';id  plot  mode,  and 
other  eontrol  variables. 


Put  phase  limits  in 
arrays  correctly 
instead  of  as  read. 


I 


V 


CXitput 

less  than  upper 
. limit 


Compute  output  from 

dif.  between  low  and 
guess,  times  lO'Jf. 

Output  X, 
greater  than  lower 
S.  limit  > 


Compute  output  from 

high  limits  times  10^. 

I Set  NCLOS  = 2 


Output 
less  than 
guess 


Change  NCLOS 
sign  to  -2. 


Compute  DEL  from  dif. 
between  guess  and 
output. 


Have  we  X 
examined  all  the  out 
Ik.  outs  that  are  non- 


I = 1,  NONLIN. 


Page 

4-1 


NOLlN-4 


1 NCLOS  =: 

fc2,  ±3 

Still 

No 

Now 

vNo  / 

1 depending  on  how  far 

^ above  upper  y 

> 

below  lower  > 

we  slippec 

1. 

limit 

^S^limit 

\ 

Yes 


Page 

5-II 


Yes 


Set  NCLOS  = -3; 

DELT  = -1/2  previous 
DELT. 


Page  5-III 


DELT  = -1/2  previous 
DELT. 


■f 

I 


NORM- I 


l)al:i  statonu'iU  del'ines 
uniform  dislrilnilion 
|X)ints. 


I.ST  0 
first  onlry 
in  run) 


SelNT()l>=l  ^ 

Compute  normal 

1ST  ^ 1 . 

distribution  points. 

N i x;  - 0 


ItV 

greater  than 
1/2 


I'ind  uniformly 
distributed  ixjints  that 
R\'  lies  between. 


Set  V.M,  - largest  nor- 
nuilly  distributed  HV  we 
have  plus  NTOP*0  1 . 


ItV  'v 
between 

first  and  sceond 
points 


ComiTute  \'AI,  using 
up|x-r  |K)inls  formula. 


P:if;e 

2-1 


NORM- 2 


Increase  NTOP  by  one 
unless  NTOP  = 19. 
Then  reset  NTOP  to  1. 


I 


St:irt 


PLOTF-I 


1 


1 


r 


1 


i 


I 

1 

I 


Yes 


Print  magnitude 
and  phase,  X 
and  V. 


Set  I'H  I^Q  iirruy 

to  0. 


Head  frequency 
designating  card. 


READFQ 


^ This  N 
a frequency 
card  ^ 


Write  diagnostic 
and  stop  run. 


Read 

frequencies. 


3-way  > 
branch  on 
V IT  123  > 


Linear  ^ 
or  logarithmic 
generation  of 
\ frequency  > 


Generate 

frequencies. 


Compute  frequencies  to 
Ije  generated.  Store  in 
DlMIi. 


Compute  number  of 
frequencies  to  be 
generated. 


' Linear  ^ 
or  logarithmic 
generation  of 
\ freq.  ^ 


Generate  frequencies 
at  decade  intervals  in 
range  desired. 


Logarithmic 


Save  NLF. 


IT  1 23  - 2 


^ 3-way  ^ 
branch  on 
V LTYI’l;  > 


Compute  NLF  and  save 
generated  frequencies 
in  DUMB. 


Read  in  remaining 
frequencies  into  DUMB 


Ihit  frequency  in 
middle  of  de- 
cades at  3 » 
initial  decade 


Sort  DI;mb  array  and  puti 
ordered  frequencies  in  f 
FRKy  array.  J 


i 


Page 

2 


SENSPR-3 


start 


SOLUT  - 1 


Set  NSCR  array 
to  zero. 


Yes 


Print  equivalent  circuit 
before  calling  CUR  ENT  and 
converting  voltage  to  cur- 
rent. Also  print  JNODE, 


NEQU. 


Print  diagnostic. 


SOLUT-2 


SOLUT-3 


I 


Sl:irl 


SOLVE  - 1 


Compute  N1  = number  of 
columns  in  augmented 
matrix.  Compute  N2, 
number  of  rows,  RK  :md 
IM. 


DC  circuit 


J = 1,  N2. 

(look  at  each  row). 


Find  max  ele- 
ment in  that 


Normalize  elements  in 
row  by  dividing  by  max 
element. 


Compare  scaled  ele- 
ments to  matrix  mini- 
mum save  smallest. 


Examined  ^s^No 
all  rows 


Set  min  = min 
element*10"'^ 


Divide  first  row  by  ajj 
to  initialize  process. 


SPCE-I 


Set  lOVKU  0. 

HV  - ro:il  nincloni 
vuriiible. 


Ih 


Find  ordinate  value  just 
above  RV  in  special 
distribution . 


None 

bigger 


Set  VAL  to 
biggest. 


Ordinate  = 
second  one 


Closer  x 
to  lower 

point  than  biggeiv 
V point 


Use  two  lower  and 
one  upper  points  in 
interpolation. 


At 

upper  end  of 
distribution 


Use  one  lower  and 
two  upper  points  in 
distribution. 


! 

i 


I 

I 

1 

! 


I 

I 

I 

I 

I 


I 

I 


Kill  Nf'llll  array  with  0 
anti  1 as  wc  did  NOUT. 


Head  in  special  value 
eard(s)  IPH  ■ 0. 


SPECIN-2 


Write  title. 


1 = 1 to  JfoMP. 


Find  equivalent  circuit 
indices  for  this  C.  K.  , 
store  in  Kl,  K2. 


i [ 

.1 


Did  s|)ecial 

'^alue  cards  indicat^ 
nominal  value  to  be  , 
used  in  this  y' 
\ C F / 


Special  value 
required,  i.e.,  M 
greater  than  , 
\ 4 / 


5 of  four  C.  F.  I 
eld's  in  ('.  F.  's  I 

F J 


Find  entries  in 
NPFC,  SPlCC  array 
and  put  correct  one 
in  C.  F.  's  IVAI.FF 


Print  line  of  C.  F. 
names,  values  and 
codes. 


-A 


STAT-4 


